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Abstract 

We consider doubly-periodic travelling waves at the surface of an in- 
finitely deep perfect fluid, only subjected to gravity g and resulting from 
the nonlinear interaction of two simply periodic travelling waves making 
an angle 2Q between them. 

Denoting by /x = gL/c? the dimensionless bifurcation parameter ( L is the 
wave length along the direction of the travelling wave and c is the velocity 
of the wave), bifurcation occurs for /i = cosd. For non-resonant cases, 
we first give a large family of formal three-dimensional gravity travelling 
waves, in the form of an expansion in powers of the amplitudes of two ba- 
sic travelling waves. "Diamond waves" are a particular case of such waves, 
when they are symmetric with respect to the direction of propagation. 
The main object of the paper is the proof of existence of such symmet- 
ric waves having the above mentioned asymptotic expansion. Due to the 
occurence of small divisors, the main difficulty is the inversion of the 
linearized operator at a non trivial point, for applying the Nash Moser 
theorem. This operator is the sum of a second order differentiation along 
a certain direction, and an integro-differential operator of first order, both 
depending periodically of coordinates. It is shown that for almost all an- 
gles 6, the 3-dimensional travelling waves bifurcate for a set of "good" 
values of the bifurcation parameter having asymptotically a full measure 
near the bifurcation curve in the parameter plane (0,^?). 
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1 Introduction 



1.1 Presentation and history of the problem 

We consider small-amplitude three-dimensional doubly periodic travelling grav- 
ity waves on the free surface of a perfect fluid. These unforced waves appear 
in literature as steady 3-dimensional water waves, since they are steady in a 
suitable moving frame. The fluid layer is supposed to be infinitely deep, and 
the flow is irrotational only subjected to gravity. The bifurcation parameter is 
the horizontal phase velocity, the infinite depth case being not essentially dif- 
ferent from the finite depth case, except for very degenerate situations that we 
do not consider here. The essential difficulty here, with respect to the existing 
literature is that we assume the absence of surface tension. Indeed the surface 
tension plays a major role in all existing proofs for three-dimensional travelling 
gravity-capillary waves, and when the surface tension is very small, which is the 
case in many usual situations, this implies a reduced domain of validity of these 
results. 

In 1847 Stokes 001 gave a nonlinear theory of two-dimensional travelling 
gravity waves, computing the flow up to the cubic order of the amplitude of 
the waves, and the first mathematical proofs for such periodic two-dimensional 
waves are due to Nekrasov |20], Levi-Civita and Struik 01] about 80 years 
ago. Mathematical progresses on the study of three-dimensional doubly periodic 
water waves came much later. In particular, to our knowledge, first formal 
expansions in powers of the amplitude of three-dimensional travelling waves 
can be found in papers ^Hl and (35]. One can find many references and results 
of researches on this subject in the review paper of Dias and Kharif ^3] (see 
section 6). The work of Reeder and Shinbrot (1981)23 represents a big step 
forward. These authors consider symmetric diamond patterns, resulting from 
(horizontal) wave vectors belonging to a lattice F' (dual to the spatial lattice 
r of the doubly periodic pattern) spanned by two wave vectors Ki and K2 
with the same length, the velocity of the wave being in the direction of the 
bissectrix of these two wave vectors, taken as the Xi horizontal axis. We give 
at Figure m two examples of patterns for these waves (see the detailed comment 
about these pictures at the end of subsection I2.4|l . These waves also appear 
in litterature as "short crested waves" (see Roberts and Schwartz |SZ|, Bridges, 
Dias, Menasce for an extensive discussion on various situations and numerical 
computations). If we denote by 9 the angle between Ki and the xi— axis, 
Reeder and Shinbrot proved that bifurcation to diamond waves occurs provided 
the angle 9 is not too close to or to 7r/2, and provided that the surface tension 
is not too small. In addition their result is only valid outside a "bad" set in the 
parameter space, corresponding to resonances, a quite small set indeed. This 
means that if one considers the dispersion relation A{K, c) = 0, where K and 
c G are respectively a wave vector and the velocity of the travelling wave, 
then there is no resonance if for the critical value of the velocity Cq there are 
only the four solutions zLKi,zLK2 of the dispersion equation, for K G T' (i.e. 
for K being any integer linear combination of Ki and K2). The fact that the 
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surface tension is supposed not to be too small is essential for being able to use 
Lyapunov-Schmidt technique, and the authors mention a small divisor problem 
if there is no surface tension, as computed for example in |37) . Notice that the 
existence of spatially bi-periodic gravity water waves was proved by Plotnikov in 
\62\ . PJj in the case of finite depth and for fixed rational values of gL/c^ta,n.9, 
where g,L,c are respectively the acceleration of gravity, the wave length in the 
direction of propagation, and the velocity of the wave. Indeed, such a special 
choice of parameters avoids resonances and the small divisor problem, because 
the pseudo-inverse of the linearized operator is bounded. 




Figure 1: 3-dim travelling wave, the elevation jy^ is computed with formula 
(frrnil . Top: = 11.3°, r = l/5,£ = 0.8^c; bottom: = 26.5", t = 1/2, e = 
0.6/Xc. The dashed line is the direction of propagation of the waves. Crests are 
dark and troughs are grey. 

Craig and NichoUs (2000) !P used the hamiltonian formulation introduced by 
Zakharov |44|. in coupling the Lyapunov-Schmidt technique with a variational 
method on the bifurcation equation. Still in the presence of surface tension, they 
could suppress the restriction of Reeder and Shinbrot on the "bad" resonance 
set in parameter space, but they pay this complementary result in losing the 
smoothness of the solutions. Among other results, the other paper by Craig 
and NichoUs (2002) [S] gives the principal parts of "simple" doubly periodic 
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waves (i.e. in the non resonant cases), expanded in Taylor series, taking into 
account the two-dimensions of the parameter c. They emphasize the fact that 
this expansion is only formal in the absence of surface tension. 

Mathematical results of another type are obtained in using "spatial dynam- 
ics", in which one of the horizontal coordinates (the distinguished direction) 
plays the role of a time variable, as was initiated by Kirchgassner and ex- 
tensively applied to two-dimensional water wave problems (see a review in 
The advantage of this method is that one does not choose the behavior of the 
solutions in the direction of the distinguished coordinate, and solutions periodic 
in this coordinate are a particular CRSC, clS well as quasi-periodic or localized so- 
lutions (solitary waves). In this framework one may a priori assume periodicity 
in a direction transverse to the distinguished direction, and a periodic solu- 
tion in the distinguished direction is automatically doubly periodic. The first 
mathematical results obtained by this method, containing 3-dimensional doubly 
periodic travelling waves, start with Haragus, Kirchgassner, Groves and Mielke 
(2001) Hi, HZ], EH, generalized by Groves and Haragus (2003) [Tl|. They use 
a hamiltonian formulation and center manifold reduction. This is essentially 
based on the fact that the spectrum of the linearized operator is discrete and 
has only a finite number of eigenvalues on the imaginary axis. These eigenvalues 
are related with the dispersion relation mentioned above. Here, one component 
(or multiples of such a component) of the wave vector K is imposed in a di- 
rection transverse to the distinguished one, and there is no restriction for the 
component of K in the distinguished direction, which, in solving the dispersion 
relation, gives the eigenvalues of the linearized operator on the imaginary axis. 
The resonant situations, in the terminology of Craig and NichoUs correspond 
here to more than one pair of eigenvalues on the imaginary axis, (in addition to 
the origin). In all cases it is known that the largest eigenvalue on the imaginary 
axis leads to a family of periodic solutions, via the Lyapunov center theorem 
(hamiltonian case), so, here again, there is no restriction on the resonant set 
in the parameter space at a fixed finite depth. The only restriction with this 
formulation is that it is necessary to assume that the depth of the fluid layer is 
finite. This ensures that the spectrum of the linearized operator has a spectral 
gap near the imaginary axis, which allows to use the center manifold reduction 
method. In fact if we restrict the study to periodic solutions as here, the center 
manifold reduction is not necessary, and the infinite depth case might be con- 
sidered in using an extension of the proof of Lyapunov-Devaney center theorem 
in the spirit of j23| . in this case where belongs to the continuous spectrum. 
However, it appears that the number of imaginary eigenvalues becomes infinite 
when the surface tension cancels, which prevents the use of center manifold re- 
duction in the limiting case we are considering in the present paper, not only 
because of the infinite depth. 

1.2 Formulation of the problem 

Since we are looking for waves travelling with velocity c, let us consider the 
system in the moving frame where the waves look steady. Let us denote by ip 
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the potential defined by 

if = (j) — c ■ 

wliere (j) is the usual velocity potential, X = {xi^X2) is the 2-dim horizontal 
coordinate, is the vertical coordinate, and the fluid region is 

n = {(X,X3) : -oo <X3< r?(X)}, 

which is bounded by the free surface E defined by 

i: = {{X,x^):x^=r){X)}. 

We also make a scaling in choosing |c| for the velocity scale, and L for a length 
scale (to be chosen later), and we still denote by (X, X3) the new coordinates, 
and by ip, rj the unknown functions. Now defining the parameter n — ^ (the 
Proude number is -^^) where g denotes the acceleration of gravity, and u the 
unit vector in the direction of c, the system reads 

infi, (1.1) 
on S, (1.2) 

on E, (1.3) 
as X3 —00. 

Hilbert spaces of periodic functions. Wc specialize our study to spa- 
tially periodic 3- dimensional travelling waves, i.e. the solutions rj and (p are 
bi-periodic in X. This means that there are two independent wave vectors 
Ki,K2 S generating a lattice 

T' = {K = niKi + n2K2 : nj e Z}, 

and a dual lattice F of periods in such that 

r = {A = miXi + m2\2 ■ iTij € Z, Xj ■ Ki — 2Ti8ji}. 

The Fourier expansions of r) and <p are in terms of e*^'^, where K ^T' and 
K ■\ = 2nn, n G Z, for A € F. The situation we consider in the further analysis, 
is with a lattice F' generated by the symmetric wave vectors Ki = (1,t), K2 = 
(1,— r). In such a case the functions on M^/F arc 2tt— periodic in xi, 27r/T— 
periodic in X2, and invariant under the shift (xi, ^2) {xi +Tr,X2+ tt/t) (and 
conversely). We define the Fourier coeflacients of a bi- periodic function u on 
such lattice F^ by 

u{k) = ^ / u{X) exp{-ik ■ X)dX. 

2"" J[0,27r]x[0,27r/T] 



Vxr?- (u + Vx</') - 1^ 

ax-i 



u • Vxf + 



(V^)^ 
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For m > we denote by H"^{M.^/r) the Sobolev space of bi-periodic functions of 
X G M^/r which are square integrable on a period, with their partial derivatives 
up to order m, and we can choose the norm as 



\ker' 



\ 1/2 

\k\r^\u{k)\A . 



Operator equations. Now, we reduce the above system for {(p, rf) to a 
system of two scalar equations in choosing the new unknown function 

and we define the Dirichlct-Neumann operator Q^^ by 



Q^^ = Vl + (Vx^)^^U=r,(X) (1.4) 

dtp , „ „ 

where n is normal to S, exterior to fi, and ip is the solution of the rj— dependent 
Dirichlet problem 

Aip ^ 0, X3< r]{X) 

ip = ip, r]{X), 
Viy9 —^ as a;3 — > —oo. 

Notice that this definition of Qj^ follows j27j and insures the selfadjointness and 
positivity of this linear operator in L^(IR^/r) (see Appendix lA.l|l . Our definition 
differs from another usual way of defining the Dirichlet - Neumann operator 
without the square root in factor in H1.4|l . Now we have the identitv l] 1.4(1 and 
the system to solve reads 

T{U,^l,u)^0, T^iT,,T2), (1-5) 

where U — {'ip,r/), and 

J'iiU,^i,u) - ■.g^i^)-u-VxV, (1-6) 

Let us define the 2-componcnts function space 

H'"(MVr) = ij™(RVr) X i7'"(MVr) 

We denote the norm of U in H™(RVr) by 
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where means functions with average, and U = {ip,r]). The average 
condition comes from the fact that the value ^ of the potential is defined up to 
an additive constant (easily checked in equations H1.6|l . (|1.7|) '). Moreover, the 
average of the right hand side of is as it can be easily checked (this is 
proved for instance in (S]). We have the following 

Lemma 1.1 For any fixed to > 3, the mapping 

{U,n,u)^ T{U,^i,u) IS C°° : H"(MVr) X M X §1 ^ H™-i(RVr) 

in the neighborhood of {0} x R x §i. Moreover J-{-,fj,,u) is equivariant under 
translations of the plane: 

where 

%U{X) = U{X + w). 

In addition, there is M3 > 0, such that for ||C/||3 < M3 and \fi\ < M^, T satisfies 
for any to > 3 the "tame" estimate 

||^([/,/i,u)||„,_i < c„AM3)\\U\U (1.8) 

where Cm only depends on m and M3. 

Proof. The C°° smoothness of (tA,??) ^ ^^(tA) : H"(MVr)_^_ff"-i(R7r) 
comes from the study of the Dirichlet-Neumann operator, see ljA.llA.21) . and the 
properties of elliptic operators. This result is proved in particular by Craig and 
NichoUs in 9 , and by D.Lannes in [57|. Notice that iJ''(R^/r) is an algebra for 
s > 1. Notice that it is proved by Craig et al ^U] that the mapping {ip,r]) 1-^ 
g^{^) : 7J™(RVr)xC""(RVr) ^ i7"-i(RVr) is analytic and the authors give 
the explicit Taylor expansion near 0, with the same type of "tame" estimates 
that we shall use in the following sections. We choose here to stay with {ip,r]) G 
IHI'"(R^/r) and we just use the C°° smoothness of the mapping, in addition to 
the tame estimates (see |77]'). 

The equivariance of J-^ under translations of the plane is obvious. 

We refer to |57] for the proof of the following "tame" estimate, valid for any 
fc > 1 (here simpler than in | 27| since we have periodic functions and since there 
is no bottom wall) 

\\GM\k < Cfc(||r/||3){||r/||fc+i||Vl|3 + llV-IU+i}, (1.9) 
necessary to get estimate (|1.8|l . ■ 

1.3 Results 

We are now in a position to formulate the main result of this paper on the 
existence of non- linear diamond waves satisfying operator equation H1.5|) . We 
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find an explicit solution to 1)1. 5|l in the vicinity of an approximate solution Ue^^ 
which existence is stated in the following lemma restricted to "diamond waves", 
i.e to solutions belonging to the important subspace (still with F' generated by 
(1,±^)) 

^fs) ^ = (^i'?) G IHl'^(IR^/r) : -0 odd in xi, even in X2, rj even in xi and in 2:2}. 

For these solutions the unit vector Uq = (1,0) is fixed (see a more general 
statement at Theorem 12. 31 with non necessarily symmetric formal solutions). 

Lemma 1.2 Let N > 3 be an arbitrary positive number and the critical value of 
parameter ^c{t) = + is such that the dispersion equation n^ + r'^m^ = 

IJ,~^n'^ has only the solution {n,m) — (1, 1) in the circle m? + < iV^. Then 
approximate 3- dimensional diamond waves are given by 

C/i^) = (0, ^)f > = J2 ^'^^'^ ^ for any k, (1.10) 

1<P<N 

i/^"'^'' = (sinxi cosrx2, — cosxi cosra;2), /Lt^^-* = /Ltc + /i, /i = /Lti£^ + 0(£^), 
where 

^ (J: 1 3 a*, 9 . 

^Ay^l 2y.l +2 4(2-/i,)^' 



and where for any k, 



£ 



N+1, 



uniformly bounded in H^g^, with respect to £. There is one critical value Tc 
of T such that /ii(Tc) = 0, and fJ,i{T) < for t < t^, /ii(T) > for r > Tc- 

Proof. The lemma is a particular case of the general Theorem 12.31 in the 
symmetric case. ■ 

The following theorem on existence of SZJ-diamond waves is the main result 
of the paper (notice that r = tan 0) 

Theorem 1.3 Let us choose arbitrary integers I > 23, > 3 and a real number 
6 < 1. Assume that 



€{6,1/6), Mc = (l + r2) 



-1/2 



Then there is a set OT of full measure in (0,1) with the following property. If 
/Zc e 91 and r 7^ Tc, then there exists a positive £0 — ^oilJ'c, N,l,6) and a set 
E — S{^c, N, I, 6) so that 

2 

lim -t; I s ds = 1, 



e^Q £ 



2 



£n(0,e) 
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and for every fi — fii^^ with e G equation (11.51) has a "diamond wave" type 



solution U = t/i^^ + e'^W, with G 



Moreover, W : £ 



IS a 



Lipschitz function cancelling at e — 0, and for t < Tc (resp. t > Tc), and when 
e varies in £ , the parameter /i — ^i^^ runs over a measurable set of the interval 
(fii^\ ^c) (resp. (/ic,/ii^') of asymptotically full measure near fic- 

We can roughly express our result in considering the two-dimensional pa- 
rameter plane (r, fi) where r = tanfl, 29 being the angle between the two basic 
wave vectors of same length generating the two-dimensional lattice T' dual of the 
lattice r of periods for the waves. The critical value Hc{t) of /i {— gL/c?), where 
/Xc(t) = (1-|-T^)~^/^ = cos6', corresponds to the solutions of the dispersion rela- 
tion we consider here (in particular 3-dimensional diamond waves propagating 
in the direction of the bisectrix of the wave vectors). We show that for r < 
(Ri 2.48) the bifurcating (diamond) waves of size 0(|/i — ^Jl^ciJ)\^/'^) occur for 
jjL < Hcij), while for r > Tc it occurs for fi > IjLc{t). We prove that bifurcation 
of these 3-dimensional waves occurs on half lines r = const of the plane, with 
their origin on the critical curve, for "good" values of r (which appear to be 
nearly all values of r). Moreover, we prove that on each half line, these waves 
exist for "good" values of /i, this set of "good" values being asymptotically of 
full measure at the bifurcation point /i = ij-c{t) (see Figure 

M. 




Figure 2: Small sectors where 3-dimensional waves bifurcate. Their vertices lie 
on the critical curve ^ = ^c(t). The good set of points is asymptotically of full 
measure at the vertex on each half line (see the detail above). In the paper we 
only give the proof for each half line r = const (dashed line on the figure) 

Another way to describe our result is in terms of a bifurcation from a non 
isolated eigenvalue in the spectrum of the linearized operator at the origin. 
Indeed, for our critical values (r, /ic(T)) of the parameter, the differential at 
the origin is a selfadjoint operator with in general a non isolated eigenvalue 
(see Theorem I4.1|l . Our result means that from each point (r, /ic(T)) where r 
is chosen in a full measure set of (0,oo), a branch of solutions bifurcates in the 
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following sense. Every half line r = const with origin at the point {t, fic{T)) 
and on the good side of this curve, contains a measurable set of points where 
the bi-periodic gravity waves exist, with an amplitude 0{\^i — /xd^/^), this set 
being asymptotically of full measure near Hc{t). 

In fact, we can improve our result in replacing the half lines mentioned 
above, by small sectors centered on these half lines. Each half line in each 
sector, with origin at the vertex of the sector, contains a measurable set where 
the bi-periodic gravity waves exist, with an amplitude 0(|/i — /ic|^^^), this set 
being asymptotically of full measure near /j.c(t). The proof of such a result 
introduces many technicalities, which are not essential for the understanding of 
the paper. This complication is mainly due to the fact that we then need to 
work with a lattice F now depending on e. We just mention in various places 
what is really needed for such an extension of the result proved here. 

1.4 Mathematical background 

There are some aspects of our method which deserve brief mention. First we use 
the Nash-Moser method, which is now an integral part of nonlinear analysis |12| . 
for proving Theorem 11.31 The crucial point for the Nash-Moser method is to 
obtain a priori bounds on an approximate right-inverse of the partial derivative 
du^{U, ^,u.). As it is shown in Section 2, this problem is equivalent to the 
problem of invertibility of a second-order selfadjoint pseudodifferential operator 
with multiple characteristics. Second, we use the Moser theory of foliation 
on a torus [29) and the invariant parametric representation of the Dirichlet- 
Neumann operator to reduce the linearized equation to a canonical form with 
constant coefficients in the principal part. We employ a modification of the Weil 
Theorem 021 on uniform distribution of numbers {um^} modulo 1 to deduce the 
effective estimates of small divisors, and hence to prove the invertibility of the 
principal part of the linearized operator. The most essential ingredient of our 
approach is the algebraic descent method j21|, [SSI which allows to reduce the 
canonical pseudodifferential equation on 2-dimensional torus to a Fredholm-type 
equation. 

1.5 Structure of the paper 

Now we can explain the organization of the paper. In section 2, we prove 
Theorem 12.31 which establishes the existence of approximate solutions under 
the form of power series of the amplitudes of the two incident mono-periodic 
travelling waves, corresponding to symmetric basic wave vectors. The parameter 
is two-dimensional here, due to the freedom in the direction of propagation of 
the three-dimensional wave. To show this result, we use a formal Lyapunov- 
Schmidt technique, assuming that the angle 20 between the two basic wave 
vectors satisfies that r — tan^ is such that the equation for positive integers 
(n, m) 

2 I 2 2 4/1 I 2\ 
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has the unique solution {n,m) = (1,1) (non resonance property). In playing 
with scales and parameters, this condition is not restrictive among non resonant 
situations, which indeed represent the general case. In such a case, the kernel 
of the linearized operator at rest state (taken as the origin) is four-dimensional, 
and in using extensively the symmetries of the system ()1.5|) . we obtain, for a 
fixed value of the bifurcation parameter (/z, u), doubly-periodic formal travelling 
gravity waves propagating in the direction u. Limiting cases are the mono- 
periodic travelling waves corresponding to one of the basic wave vectors. The 
Lemma ll.2l is a particular case of the above theorem. 

From now on, we restrict the study to solutions called "diamond waves", 
which are symmetric with respect to the direction of propagation, here the xi — 
axis. In section 3 we consider the linear operator C{U, /j,) corresponding to the 
differential of H1.5() at a non zero point in H^^^ , which we need to invert for using 
the Nash-Moser theorem. The principal part of this operator is the symmetric 
sum 

a 

of a second order derivative in the direction of a periodic vector field V{X), 
and of the Dirichlet-Neumann operator which is integro-differential of first or- 
der, both parts depending periodically on coordinates. More precisely, V = 
(G-'^{X){ua+VtpiX)), where G(X) dX ■ dX is the first fundamental form of the 
free surface. Recall that G is a covariant tensor field on E, and for the standard 
parametrization X3 = r]{X), it is given by G(X) = 1-1- Wrj (g) V77. It follows 
from the kinematic condition H1.2() that integral curves of the vector field V{X) 
coincide with trajectories of liquid particles moving along S and submitted to 
the vertical gravity /x. 

Section 3 is concerned with the first step of the long way towards the in- 
version of £, which consists in finding a diffeomorphism of the torus for which 
the highest order terms of the operator C become constants (depending on the 
linearization point). We begin fLemma I3.6|l with the construction of a dif- 
feomorphism which takes integral curves of the vector field V onto straight 
lines parallel to the abscissa axis. Being endowed with the Jacobi metric 
ds^ = (1/2 — i.iri{X))G{X) dX ■ dX the free surface becomes a Riemannian man- 
ifold on which the integral curves of V coincide with geodesies (see Appendix 
nj. Hence, by Lemma [3.61 they form a geodesic foliation on S. Moreover, since 
the distance between each of these curves and the abscissa axis is finite, the 
foliation has a zero rotation number. It is at this point where the restriction 
to symmetric solutions (diamond waves) is necessary, since we don't know yet 
how to manage such a diffeomorphism in the non symmetric case, see j29j for 
discussion. Recall that the Moser Theorem 29 guarantees the existence of at 
least one geodesic for any given rotation number. 

The second result of Section 3 is Theorem 13.51 which gives the parametric 
representation of the Dirichlet-Neumann operator in arbitrary coordinates Y on 
E so that a mapping X — X{Y) is a diffeomorphism of a torus. It follows from 
this theorem that for any smooth periodic function u{Y) and ii{X) = u{Y{X), 
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the Dirichlct-Ncumann operator has the decomposition 

Gr^u = GiU + QoU + G-lU, 

in which Gi are pseudodifFerential operators of order i. We give an exphcit ex- 
pression for their symbols in terms of the first fundamental form and the princi- 
pal curvatures of the free surface. In particular, we show that, up to a positive 
invariant multiplier, the symbol of the operator Gi is equal to ^yG{Y)^^k ■ k, 
and the real part of the symbol of Go coincides with the difference between the 
sum of the principal curvatures and the normal curvature of E in the direction 
of G~^k. This leads to the interesting conclusion: the manifold E is defined by 
its Dirichlet- Neumann operator up to translation and rotation of the embedding 
space. 

Combining Lemma |3 . 61 and Theorem 13.51 gives the main result of Section 3 
- Theorem 13.41 This theorem ensures the existence of a diffeomorphism X = 
X{Y) of the 2-dimensional torus, which brings the linearized operator to the 
canonical form 

2. + f) = Z + my,+'& + £.-i, 

where the remainder £_i is of order —1, 21 and 05 are zero-order pseudodifFer- 
ential operators, and the principal part 

£ = < + (-A)V^ ^ = dl+r'dl 

is a selfadjoint pseudodifferential operator. Here the parameter v depends on 
the point of linearization, with vifi) = vq = ^c{t)~^ (Lemma 13. 7|) . 

In Section 4 we study the operator £ in many details, and give estimates on 
its resolvent in Sobolev spaces of bi-periodic functions which are odd in yi , and 
even in 1/2 ■ We begin with the observation that for vq = ^c{t)~^ and almost 
every positive r, zero is a simple eigenvalue of the operator £0 = i^oS^j A)^/^ 
and 

||£o^u|U < c{T,a)\\u\\s+(l+a)/2 

for all u orthogonal to the kernel of £0 and a > 0. Next we study the pertur- 
bation of its resolvent assuming that v = vq — e^vi + 0{e^) and with a spectral 
parameter k — 0(£^), both being Lipschitz functions of a small parameter e. 
Here we have a small divisor problem, and we meet the necessity to restrict 
the parameter values to "good ones", for being able to find suitable estimates. 
Calculations fLemma I4.5|l show that the resolvent of £ satisfies the estimate 

>i)-^u\\s<c\\u\\,+i, ue(ker£o)^, (1.11) 

if parameters v and k satisfy the quadratic Diophantine inequalities 

|ijjn^ — m — C| > crt^^ for all positive integers n,m, (1-12) 

where w = i/r~^ and C = {2vt)^^ — kt^^ . Note that there is a difference 
between linear and polynomial Diophantine approximations: in classic theory 
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of linear Diophantine forms, see [H] for general theory, the integers for which 
"small divisors" are really small, form a sparse set in the integral lattice. This 
property was used in pioneering works of Siegel and in the Arnold proof of 
the Kolmogorov Theorem [2]. In contrast to the linear case, the couples (m, n), 
for which the left hand side of inequality (|1.12|) is small, can form clusters 
in Z^, and the problem of obtaining small divisors estimates becomes more 
complicated. It turns out that the validity of inequalities ()1.12() with a constant 
c independent of the small parameter is a consequence of estimate 

card {n : ujon^ modulo 1 < e and 1 < n < N} < ce for aU N > e^^. 

(1.13) 

Recall that, by the Weil Theorem 03], [H], for each fixed e, the left hand side 
tends to e as A'^ — !■ cx). Hence the inequality holds true for some c depending on 
e. In Appendix^ we make this result more precise and prove the existence of 
absolute constant c such that inequality H1.13|l is fulfilled for all A > 78 and all 
intervals of length e. This leads to the main result of this section - Theorem 
14.21 which shows that with a suitable choice of the parameters, the resolvent 
operator provides a loss of one in the degree of differentiability. Moreover, 
estimate Hl.ll|l holds true for all in an asymptotically full measure set on 
every half line 9 ~ const of the parameter plane, the origin of which being 
chosen arbitrarily in a full measure set, on the bifurcation curve fi — cos6. 

In Section 5 we take into account all remaining terms of the linear operator 
£ + i3 and prove its invertibility with a loss of differentiability. The main diffi- 
culty is that the operator £ + 9j involves the principal part £, which inverse is 
unbounded, and arbitrary operators 21, 05 with "variable coefhcients" . 

Most, if not all, existing results related to such problems were obtained by 
use of the Frdhlich-Spencer method proposed in JS], cf and developed 

by Craig and Wayne ^2 Ej and Bourgain W . The basic idea of the method 
is a representation of operators in the form of infinite matrices with elements 
labelled by some lattice and block decompositions of this lattice. Let us use the 
operator £ -I- to illustrate the main features of this method. First we have to 
replace a periodic function u by the sequence of its Fourier coefficients 
k € r', and the operator £ by the diagonal matrix with the elements L{k). 
Then we have to split the lattice F' into a " regular" part which consists of all k 
with "large" L{k), and an "irregular" part which includes all k corresponding 
to "small" values of L{k). Using the contraction mapping principle we can 
eliminate the " regular" component and reduce the inversion of £ + ^ to the 
inversion of an infinite matrix on the "irregular" subspace. The existence of an 
inverse to this matrix is established by using a special iteration process which 
is the core of the method. Note that the Frohlich-Spencer method is working in 
our case only if 2t = 0. 

Our approach is based on the descent method which was proposed in |35[I24| 
and dates back to the classic Floquct-Lyapunov theory. The descent method is 
a pure algebraic procedure which brings the canonical operator to an operator 
with constant coefficients and does not depend on the structure and spectral 
properties of the principal part £. The heart of the method is the following 
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identity fTheorem l5.2|l 

(£ + 2is)i + »)(i + e:)ii== {1 + ~ >c)u + ^u, 

which holds true for all functions u e i/^(M^/r) odd in yi. Here £ and £ are 
bounded operators in the Sobolev spaces of periodic functions H^{R'^/r); the 
remainder S' is a bounded operator : ^'^-^(K^/r) i-^ H''(M?/T); the Floquet 
exponent >c has an explicit expression in terms of operators 2t and 05. More- 
over, if ||2l||, ||*8|| - £, then ||e:||, ||£||, Hg-H - e and ^< = 0{e^). The proof of 
these results constitutes Section 15.11 and Appendix The technique used is 
more general than the one used in 24 , since we use here general properties 
of pseudodifferential operators, however taking into account of the symmetry 
properties of £. 

The descent method of algebraic character presented here, might be eas- 
ily used for example on the one-dimensional KDV and Schrodinger equations, 
avoiding the heavy technicalities of the Frohlich-Spencer method. 

Thus we reduce the problem of the inversion of the canonical operator £ + 
to the problem of the inversion of operator £— >f+£_i, where £_i is a smoothing 
remainder. It is then possible to use the result of section 4 for inverting the full 
operator and to prove Theorem 15.11 - the main result on the existence and 
estimates of (£ -I- In particular, this theorem implies that if 21 and S 

are Lipschitz operator-valued functions of a small parameter e, which vanish 
for £ = and satisfy symmetry and metric conditions (Section 5), and if £ — >f 
meets all requirements of Theorem l4.2l then for all e"^ taken in an asymptotically 
full measure set, the resolvent has the representation 

(£ + i5)-\£) = -<oo{e)+^i{e), 
c 

in which operators S)i{e) : iJ'*(R^/r) i-^ iJ''^-^(M^/r) are uniformly bounded 
in £, and SjQ{e) are bounded operators of rank 1, the coefhcient c — £^@ -I- 
0{e^) being given by (|5.9|l . We show at the end of the section (see Theorem 
I5.9|l . that the results apply to the linear operator £([/, /i) corresponding to 
the differential of 1)1. 5|l at a non zero point in H^^^. In particular, we give 
the sufficient conditions which provide the existence of the bounded inverse 
£(f/,/i)-i : iJ^+3(R/r) ^ H'iR/T). 

Section 6 applies extensively the result proved in |22 concerning the Nash- 
Moser theorem with parameters in a Cantor set. The main result, which is the 
main result of the paper is Theorem 11.31 establishing the existence of smooth 
bi-periodic travelling gravity waves symmetric with respect to the direction of 
propagation, in the region of the parameter space mentioned above. Notice 
that, a part from the last section, which heavily rests upon the self contained 
Appendix N of the rest of the paper is self contained, with some details of 
computations and basics on pseudodifferential operators put in Appendix, for 
providing an easy reading. 
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2 Formal solutions 

2.1 Differential of Q^i 

In this subsection wc study the structure of the operator C/^, and we give useful 
formulas and estimates. 

The following regularity property holds 

Lemma 2.1 The differential h i-^ d,qQri[h] of Grj satisfies for rj, ip, h smooth 
enough bi-periodic functions 

dr,gn[h]^P = -gr,{hO+yx-{iCyxV-^x4')h}, (2.1) 
^ = i + (v^^)2 {g.V' + Vx7?-V^^}. (2.2) 

Moreover, despite of the apparent loss of derivatives for (ip, rj) in 1^2. 112. ^) we 
have for ||?7||3 < M^, the following tame estimate 

WdriGrimWk < Cfc(Af3){||/^|U+l + ||C/||fc+l||/l||2}. 

Proof. We refer to Appendix lA. II for the formula H2.1I2.2|1 already proved for 
instance in |27], and we also refer to (27] for the tame estimate. ■ 

From the formulas of the above Lemma 12.11 we are now able to compute 
successive derivatives of JF. Observe that in H2.1I2.2|I there is a loss of two 
derivatives for {ip^rf). In fact there is a compensation cancelling the dependence 
into the second order derivatives and we have the following Lemma which com- 
pletes Lemma Fl. II 

Lemma 2.2 For ||?7||3 < M^, and < M3 the following tame estimates hold 
( and analogous ones for higher order derivatives ) 

\\^uT{U,^i,^x)[SU]\\k < Cfe(M3){||<5C/||fc+i + ||C/||fe+i||,5[/||3}, 
\\dlruHU,t^,n)[SUi,6U2]\\k < Cfc(M3){||(5t/i||fc+i||5t/2||3 + 

+ ||5C/2||fc+l||<5C/l||3 + ||C/||fc+l||<5f/l||3||5C/2||3}. 

2.2 Linearized equations at the origin and dispersion re- 
lation 

The linearization at the origin of system H1.6|l . (|1.7|l leads to 

g(o)(V)-u. Vx?? = 0, (2.3) 
u-Vx^ + HV = 0, (2.4) 

where the following operator 

= (-A)i/2 
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is defined more precisely in Appendix lA. 21 Now expanding in Fourier series, we 
have 

Kev iter' 
and (ESI, (E3|) give for any K (^V 

\K\'iIjk-i{K-m)iik = 0, 
i{K ■ u)ipK + mK = 0. 

Hence, the dispersion relation reads 

^{K,^,,u)^M ^\K\-{K-nf=i). (2.5) 

The point now to discuss is the number of solutions K e T' of (|2.5|l . for 
a fixed vector u g §i, and a fixed parameter /i. We restrict our analysis to a 
lattice r' generated by two vectors Ki and K2 symmetric with respect to the 
Xi— axis, taken in the direction of u, which is the situation if one is looking for 
short crested waves: 

A'i-(l,r), X2 = (1,-t) 

where r is positive. When r is small, the lattice T of periods is formed with 
diamonds elongated in the X2 direction (see Figure . Taking 1 for the first 
component of Ki implies that we choose the length scale L as the wave length 
in the xi— direction divided by 27r. 

We consider in what follows, the cases when the direction Uq of the travelling 
waves at criticality is the xi— axis, and the critical parameter /ic = (1 + t^)~^/^ 
is such that the equation for (mi,m2) G 

IJ'c\Jm\ + r^mj — m\ = Q (2.6) 

has only the solution 

{mi, 1^2) = (1, 1). 

In case we have a solution (mi, m2) ^ (1, 1), one can make the change (/ic, t) i— > 
(mT' ^mf ) recover the case we study here. Moreover, changing /i into /i/mi 
corresponds to changing the length scale L into L/mi which indeed corresponds 
to the new wave length in the xi direction. So, it is clear that we do not restrict 
the generality in choosing the case of a solution (mi, m2) = (1,1)- 

Notice that for any integer Z, when t — I ot: l/l, there is an infinite number 
of solutions (toi,TO2) of (|2.6|l . hence we need to avoid such choices for t. 

Remark. We notice here the fundamental difference between the present 
type of study and the works using spatial dynamics for finding travelling waves, 
as for instance Groves and Haragus in 18 . Their study only consider cases 
with surface tension, and cannot work without surface tension, since this would 
lead to an infinite set of imaginary eigenvalues ±zmi (hence preventing the 
use of center manifold reduction), with no restriction for mi to be an integer, 
while m2 G N (this corresponds to fixing the length scale with the period in X2 , 
transverse to the direction of the travelling waves. 
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2.3 Formal computation of 3-dimensional waves in the 
simple case 

In this subsection we make a formal bifurcation analysis for the simple case. 
We denote by fic the critical value of /i, and we denote by Uq = (1,0) the 
critical direction for the waves (this direction of propagation may be changed for 
bifurcating travelling waves). The lattice F is generated by the two symmetric 
wave vectors Ki = (l,r), K2 — (1, — r), where 

We notice that we have the following Fourier series for U — {tp, rj) : 

n={ni ,n2)GZ2 

and we notice that 

niKi ■ X + n2K2 ■ X = miXi + Tm2X2, 

mi = ni+n2, m,2 = ni - n2, 

which gives functions which are 27r— periodic in xi and 27r/r— periodic in X2. 

We already noticed the equivariance of system (|l.t)|l . (|1.7|l with respect to 
translations of the plane, represented by the linear operator v being any 
vector of the plane. Let us complete the symmetry properties of our system 
by the symmetries Sq and Si defined by the representations of respectively the 
symmetry with respect to 0, and the symmetry with respect to Xi axis 

SoU = J2 (^(:^n)e-^("^^-^+"^^-^'\ 5C/„ = (-V'„,r?„),(2.7) 

SlU = J2 t/„e*("i^^-^+"^-^i-^^ (2.8) 

n={ni ,n2)&^ 

The system (|1.6|) . (|1.7f) is equivariant, under the symmetry Sq in all cases, while 
it is equivariant under Si only if 

U • if 1 = U • 7^2- 

In particular, in such a case we have 

where we denote by £0 the symmetric linearized operator for ^ — fj,c and u = Uq 



uo • V He 

Notice that the commutation property with the linear operator Sq is not trivial 
since the choice of writing our system in the moving frame selects the direction 
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u which breaks a reflection symmetry. Indeed the symmetry property results 
from the galilean invariance of the Euler equations. 

In the following Lemma, we use the two parameters fl = ii—j-ic and u) = u— uq 
and we notice that, since u is unitary, we have 

uj = (a;i,W2), 

W2 = • {Ki - K2), 

It 

-1 = -^+0(0.1). 
In the following subsections we prove the following 

Theorem 2.3 Assume we are in the simple case, i.e. for r and the critical 
value of the parameter ndj) = (1 + r^)"-^/^ < 1 such that the equation v? + 
T^TT? = Hc^n'^ has only the solution (n, to) = (1, 1) in N^. Then, for any N >1, 
and any v e K^, approximate 3-dimensional waves are given by %rUe^^ (T"^— 
torus family of solutions) where Ki = (l,r), K2 = (1, — r) are the wave vectors, 
and 

[/W = (V,,?7)W= J2 el^el^U^'P^'P^^ &m^, for anyk, 

f/(i.o) = = (sin(ii'i • X), — cos{Kr ■ X)), 
i7(°'^) = 6 = («in(A'2 • X), — cos(if2 • X)), 

C/^''°^ = (^^sin(2i^i.X),^cos(2i^i.X)^, 
= (^sM2i^2-X),^cos(2i^2-X)), 

+ ("o, — cos((i^i - i^2) • X) 

A = -^(ao+/3o)(£?+e^) + 0{(c?+ci)2}, 
u-{K,-K2) = {el-el)(^^{ao-po) + 0{{el + el)} 

with 



4/3 2 3„„ 9 
Po-ao = [ ? + ^H S-^IJ-c 



K V Mc 2 - /ic 
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and where for any N and k 

^(C/i^),Mc + A, uo + c^) = [el + sD^Q,, 

Qe uniformly bounded in H'^, with respect to 61,62- There are critical values 
Tc and t'^ of t such that (qq + I3q){tc) — 0, and (ap — Po){Tc) = 0, and ap + Po 
is positive for r g (0, Tj,), negative for r > Tc, Tc ~ 2.48, while — /3o is 
negative for r G (0, r^), positive for r > r^, w 0.504. Moreover, for ei = £2 
we /lawe "diamond waves" where the direction of propagation u is along the xi — 
axis; and for £2 = (resp. £1 = 0) we obtain 2-dimensional travelling waves 
of wave vector Ki (resp. K2)- All solutions are invariant under the shift %,o '■ 
X i-^ X + (tt, tt/t) and ^pi^^ is odd in X, while r]i^^ is even in X. 

Remark 1: In this lemma we assume that equation (12.6(1 for (mi, m2) G 
has only the four solutions (mi, 7712) = (±1,±1), corresponding to the four 
wave vectors K = ±Ki and ±i^2- The corresponding pattern of the waves for 
£1 = £2 is in diamond form, and for r close to 0, the diamonds are flattened in 
the xi direction, and elongated in X2 direction, looking like flattened hexagons 
or flattened rectangles, because of the elongated shape of crests and troughs. 
This last case is indeed observed experimentally for deep fluid layers. Nearly 
all (in the measure sense) values of r are indeed such that we are in the simple 
case. 

Remark 2: The above Lemma is stated differently in Theorem 4.1 of |H|; 
indeed we prove here that the manifold of solutions has a simple formulation in 
terms of the two parameters. 

Remark 3: Since ^ = gL/c^, the result of Theorem 12.31 about the sign 
of ao + Po shows that for t < Tc the bifurcation of "diamond waves" (i.e. 
a; = 0) occurs for a velocity c of the waves larger than the critical velocity cq 
corresponding to fic, while for r > Tc the bifurcation occurs for c < cq. This 
is in accordance with the numerical results of Bridges et al |1] (see p. 166-167 
with Al = Al real. Tip = 1,t = (no surface tension)). Notice that « 6.15, 
i.e. this corresponds to a critical angle Q between the wave number K\ and the 
direction of the travelling wave, such that d ~ 68*^, which is very large, and not 
easy to reach experimentally. 

Remark 4: Notice that for r near we still have "diamond waves' (even 
for T = T^) and it may exist other bifurcating 3-dimensional waves, as noticed 
in However to confirm this, we need to compute at least coefficients of order 
4 in |R5|) . 

Proof : the proof is made in Appendix ^ 
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2.4 Geometric pattern of diamond waves 

Diamond waves are obtained for Si = £2, they propagate along the xi— axis, 
and possess the symmetry Si. In making £1 = ejl in Lemma 12.31 we obtain 



f/J^) = (V',?7)W= ^ e H^ for any /c, (2.10) 

nGN, n<N 

U^^'' = ^0 = (sinxi cosTa;9, cosa;i COSTX2), 

{/(2) ^ 1(^(2,0)^^(0,2)^^^(1,1))^ 

where V'e^'' G Hj^ ^, rji^^ G iJg g meaning that V'e^'' is odd in xi, even in a;2, and 
rji^"^ is even in both coordinates. Moreover, we have 



with 



= (l + r2)-i/2, 



A^i = -y|(«o + /3o) (2.11) 
1 1 3 ^ 9 



2^il 2 4(2 -;.,)■ 

For r close to (corresponds to some of the experiments shown in |2()jV one has 

77 = ecosxi cosrx2 + — cos2xi(l + cos 2x0:2) + 0{e^T^ + \e\^) 

with 

e = -e//ic ~ -£■ 

One can assume that £ > since £ < would correspond to a shift by tt of xi. 
Then the above formula shows that crests (maxima) and troughs (minima) are 
elongated in the X2 direction, with crests sharper than the troughs, and there 
are "nodal" lines {rj k, 0) (as noticed in experiments ^D]) at 

X2 — 7r/2T + WK /t^ n G Z, 

where 77 is of order 0{e^T'^ + So, the pattern roughly looks asymptotically 
like rectangles elongated in xi direction, narrow around the crests, wide around 
the troughs, organized in staggered rows. Notice that when r — > 0, we have 
/zi ~ —3/4, hence ^ (4/3)(/ic — m), i.e. 

Co 3 2-2 7^ Co 



CO 8 ' " VgL' 
where Fq is the Froude number built with the short wave length L. 
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For larger, still small, values of r, the nodal lines disappear and the pattern 
looks like hexagons where two sides of crests, parallel to the X2 axis, are con- 
nected to the nearest tip of two analogue crests, shifted by half of the wave 
length in xi and X2 directions. 

For values of t near 1, the pattern of the surface looks like juxtaposition of 
squares. 

For large r, i.e. in particular t > Tc, we have for the free surface truncated 
at order 

77 = ecosxi cosTa;2 H — — cos 2x2:2(1 + cos2a;i) + O(e^) 

where 

£ = —e/ He ^ —ST. 

The above formula shows that crests (maxima) and troughs (minima) are elon- 
gated in the Xi direction, with crests sharper than the troughs, and there are 
"nodal" lines (?/ w 0) at 

Xi — 7r/2 -|- nTT, n ^ "L. 

Moreover, in the above formula, we see that r is in factor of e^, which means that 
the second order term in the expansion influences much sooner the shape of the 
surface as e increases. In particular for small values of e there are local maxima 
between two minima in the troughs. This phenomenon is seen in experiments 
(see 1201) in the case when t is small, however the values of e allowing such 
a phenomenon for r small are 0(1) and cannot be justified mathematically. 
So, when r is large, the pattern roughly looks asymptotically like rectangles 
elongated in xi direction, narrow around the crests, wide around the troughs, 
organized in staggered rows, and where local maxima in the middle of the troughs 
may occur for a large enough amplitude. Notice in addition that when t — > 00, 
then fii ^ ^t''/^, hence e^r^ ~ 7^(m ~ Mc), i-e- 



Co 



CO 8' ' V^' 

where we observe that in this last formula L is the physical wave length along xi, 
which is in this case the long wave length {— tLi, if Li denotes the short one). 
We plot at Figure n the elevation r/f^ a) for r = 1/5 {0 - 11.3°), e = O.S^c, 
b) for T — 1/2 {9 ^ 26.5°), e = 0.6/ic. These cases correspond to r very small 
or moderately small, currently observed in experiments (see [20)1. Observe 
however that in both cases we need to consider t not exactly 1/5 or 1/2 since 
both cases are not "simple cases" as required at Theorem 12.31 Indeed, if we 
consider solutions (n,m) G different from (1,1) for the critical dispersion 
relation 

+T^m'^ = {1 + T^)n'^, (2.12) 

then the smallest values are {n,m) — (61,18971) for t = 1/5, and (13,377) for 
r = 1/2. This means that the computation fails for coefhcients of 
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first case, and e^^'' in tiie second case. Notice that the formal computations 
made by Roberts and Schwartz 137' (1983) correspond to diamond waves with 
T — 1 {9 = 45°) and r = \/3 {9 — 60°). They observed numerically the flattening 
of troughs and sharpening of crests. However, both cases are not "simple cases" 
in the sense of Theorem 12.31 and these formal computation should break at 
order e^^ for t — 1, and at order e^^^ for r = -^3 (due to solutions of H2.12|l 
(n, to) = (5, 35) for t — 1, and (n, m) = (13, 195) for r = VS). 

3 Linearized operator 

In this section we study the linearized problem at a non zero U = (iptV) for the 
system (II. 6|) . (|1.7|l . We restrict our study to diamond waves, i.e. the direction 
of the waves we are looking for is u = uq = (1, 0) and the system possesses the 
symmetries So and Si (see ()2.7|l . H2.8|l V The purpose is to invert the hnearized 
operator, for being able to use the Newton method, as it is required in Nash- 
Moser theorem. 

3.1 Linearized system in {ip,r]) ^ 

Let us write the nonlinear system 11. 6() . H1.7|l under the form (|1.5(l 



and we omit the argument Uq since it is now fixed. Then, for any given (/, g) 
the linear system 



Wm) -0, 



where 



f/ = (V',r7), 



duT(lJ,ii)\m\=F, F:^{f,g) 



can be written as follows 



V ■ V((5Vj) + fiSri + (b^V?/ - b(VV' + uo)) 



/, 

V(<5?7) -5 



where 



y = VV^ + Uo - bV?/, b 



1 



{V?? 



(uo + VV')}. 



(3.1) 



1 + |V77|2 




and after using (|2.1I2.2II . we obtain the new system 



£{U,n)[S(l),dT]] =F + 'R{T, U)[5Ul 



(3.2) 



where the linear symmetric operator C{U, /i) is defined by 




(3.3) 
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J:^V-V{-), a = V-Vh + ii. (3.4) 
The rest TZ has the form 

n{T,U)[5U] = {Ri{T,U)[5U],Q), 

and cancels when [/ is a solution of J-{U, /i) = 0. We also notice that for U = 

(^Z^, ri) e HJ^) where 

M"g^^ = |[/ e H"'(M^/r) : odd in xi, even in 0:2, r] even in xi and in 2:2}: 
then 

V = (14,^2) e H,"ri(RVr) X i/™-i(RVr), 

all these functions being invariant under the shift 

%o : {xi,X2) i-> (xi +7r, X2 +tt/t). 
Moreover we have the following "tame" estimates 

Lemma 3.1 Let U E U"^), m > 3. Then, there exists A/3 > such that for 
\\U\\3 < M3, one has 

||y-u|U_i + ||a-/x|U_2 < c„(Af3)||C/|U, 

\\Ri{J^,U)[6U]\U^2<CsiM,){\\J^,UM\^^^^^ 

Proof. The tame estimates on — u and a ~ fi result directly from their 
definitions, from the following inequality, valid for any /, 5 G iJ'"(]R^/r), to > 2 

||/5lU<c™{||/||2||5|U + ||/|U||5l|2}, (3.5) 

and from interpolation estimates like 

||/IUo+(i~A)^<c||/||^||/||^-\ (3.6) 

which leads to 

||/IUJ|/IU<c||/|UJ|/IU 

when 

< ai < Q!2 < /32 < ai + /?! = Q!2 + (32- 

The tame estimate on Ri{J-, U)[SU] follows from the tame estimates (|3.5ll . (I1.9|l 
and from the following interpolation estimate deduced from (|3.6() : 

||/||3||5lk-+l<Cfe{||/||2||.9|k.+2 + ||/|k.+ l||5l|3}. 
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The main problem in using the Nash Moser theorem, is to invert the ap- 
proximate hnearized system, i.e. invert the hnear system 



£[J0,(5r,] = (/,5), 
which leads to the scalar equation 

Q^{5(l>)-J*{-J{5(ly))^h (3.7) 
a 

with 

a 

and where we look for 64) in some H^^J'i^^ /T). 

3.2 Pseudodifferential operators and diffeomorphism of 
the torus 

In this section we use a diffeomorphism: ^ , such that the principal part of 
the symbol of the linear operator occurring in (|3.7|l has a simplified structure. 
Its new structure will allow us to use further a suitable descent method for 
obtaining, at the end of the process, a pseudodifferential operator equation 
with constant coefhcients, plus a perturbation operator of "small" order. 

Let us denote the change of coordinates hy X = X{Y), where X{-) is a not 
yet determined diffeomorphism of such that 

X{Y) = TY + V{Y), TY = {v^,V2/t), (3.8) 

We assume that V{Y) and fj{Y) — ri{X{Y)) satisfy the following 

Condition 3.2 Functions fj and V are doubly 27: -periodic, rj is even in yi and 
2/2; Vi is odd in yi and even in y2, V2 is odd in 2/2 and even in yi, and 

Vivi +7r,t/2 +7r) = r~i{yi,y2), V{yi + tt, y2 + tt) = V{yi,y2). 

In particular, X{Y) takes diffeomorphically R^/Fi onto M.'^/Tr- The lattices 
of periods Fi and F^- are respectively the dual of lattices F'^ and F^ generated 
by the wave vectors (1, ±1) and (1, ±t). In new coordinates the free surface has 
the parametric representation 

x = riY) := {X,iY),X2iY),?jiY)y, 

with 2x2 matrix Gr{Y) of the first fundamental form of the free surface defined 
by dij — dy-T ■ dy-v. We denote by J the determinant of the Jacobian matrix 

B(y) = VyX{y). 

Our aim is to simplify the structure of the operators involved in the basic 
equation (|3.7(l by choosing an appropriate change of coordinates. The most 
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suitable tool for organizing such a process is the theory of pseudodifferential 
operators, and we begin with recalling the definition of a pseudodifferential 
operator . We consider the class of integro-differential operators on a two- 
dimensional torus having the representation 

^"(^) = ^ E e"'''A{Y,k)u{k), m = i- / e-"^''u{Y)dY, 

which properties are completely characterised by the function A : x i-^ C 
named the symbol of 21. We say that 21 is a pseudodifferential operator, if its 
symbol satisfy the condition 

Condition 3.3 There are integers / > 0, m > and a real r named the order 
of the operator 2t so that 

|2t|L.i = P(-,0)||c'+ sup sup |fc|l"l-'-||a,"A(.,fc)||c< <oo. 

fceZ2\{0} \a\<m 

Pseudodifferential operators enjoy many remarkable properties including ex- 
plicit formulae for compositions and commutators (see Appendix^for references 
and more details). Important examples of such operators are the first-order 
pseudodifferential operator Q'^^^ — (— A)^/^ with the symbol |T^^fc| and the 
second-order pseudodifferential operator 

£ = j^D? + (-A)i/2, where Di = a^,,, (3.9) 

with the symbol 

L{k)^-vkf + \T-'^k\. (3.10) 
On the other hand, integro-differential operators defined by 



are not pseudodifferential for j < 0, which easy follows from the formulae 

Diu = dl u for j > 0, and D^u — Hiu — u / m(s, y2)ds. 

27r J 

Further we will consider also the special class of zero-order pseudodifferential 
operators 21 with symbols having the form of composition A{Y^ C(fc))j where the 
vector field ^(fc) = , £,2{k)) is defined by 

^(fc) = T-^k/\T-^k\ for fc 7^ 0, ^(0) = 0. 

The metric properties of such operators are characterized by the norm 

\%U,i= sup sup ||a|'A(-,e)||c' <oo, (3.12) 

\a\<m |?|<1 
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which is equivalent to the norm |2t|^ By abuse of notation, further we will 
write simply f instead of ^(fc) and use both the notations A{Y, ^1,^2) and A{Y, ^) 
for A{Y,^). 

We are now in a position to formulate the main result of this section. 

Theorem 3.4 For any integers p,m> 14, real t £ (6,6^^), and U — {ip,'n) € 
H|^^j(M2/r^), there exists Eq > so that for 

\\U\\p < e with e e [0,eo] 

there are a diffeomorphism of the torus of the form (|3.8|l . satisfying Condition 
zero-order pseudodifferential operators 21, *B and an integro- differential op- 
erator £_i of order —1 such that: 



(i) The identity 



a J 



£u + 2lDiM + S-it + £_iM 



(3.13) 



holds true for any u G Hl^^{B?/Ti) and u{X) = u{Y{X)). 
(a) The operators 2t, *B and £_i have the bounds 

|2lkm-6 + 1*814 ,m— 6 <c,„(||C/||6)|lC/|l 
||£_iu||r < c£||M||r-i, for I < r < p — 13, 
< c(£||u||s_i + ||C/||,5+i3||m||o). 



(Hi) Operators %, 05 and £_i are invariant with respect to the symmetries 
Y ±Y* , Y* ~ (~yi,?/2) which is equivalent to the equivariant property 



2lDiu(±r*) = 2lDiu*(±r), Q3u(±y*) = ^u*(±Y), 



(3.14) 



£_iu(±y*) ^ £-iu*{±Y), u*{Y) = u{Y*); 
they are also invariant with respect to transform 1" — > 1" + (tt, tt). 
(iv) Diffeomorphism (|3.8|l of the torus can be inverted asY = T^^(X — V{X)), 

yi = Xi +d{xi,X2), y2 = TX2 +Te{xi,X2)- 



Functions d £ CJ! 



'/Tr), e e C™-4(MVr,), K,J e Cl^-H^V^i) 



and parameter v satisfy the inequalities 

llrfllc— 4 + llellc"-* < c™(||;7||4)||C/||™, \v - 1/mI < c(||{7||4)||;7|l4, 

(3.15) 



\k~ 1||c,„-5 + II J - 1/t||c..-5 < c™(||;7||5)||C/|| 



(3.16) 



27 



The proof is based on two propositions, the first of which gives the represen- 
tation of the Dirichlet-Neumann operator in the invariant parametric form, and 
the second shows that trajectories of hquid particles on the free surface forms 
a fohation of the two-dimensional torus. 

In order to formulate them it is convenient to introduce the notations 

Gi {Y, k) = \/(G-ifc • fc, div q(r) = ^ divy (Vdet Gq(r)) . 

Vdct (G 

Recall that Gfc • k is the first fundamental form of the surface S. 

Theorem 3.5 Suppose that functions f) and V satisfy Condition \cl.<l\ and there 
are integers p, I such that 

M\c^ + \\V\\cp < e, 10<p</, 
||^?||c= + ||V||c= < Eu s<l. 

Then there exists Eq > depending on p and I only such that for < e < eo o.nd 
2TT-periodic sufficiently smooth function u , the operator Q^, has the representation 

g,^uo{T + v) = g^u + gnu + g^iu, u{x) = u{y{x)). (3.17) 

Here C/i is a first order pseudodifferential operator with symbol 



Gi(r,fc) = ^-M^Gi(r,fc), YeM^, kel?, (3.18) 
C/Q is a zero order pseudodifferential operator with symbol 

Go = Re Go + ilm Go , 
Re Go(y, k) = [ ^ Q(r, k) + div(G"iVyr?) ] , (3.19) 

Im Go(r, k) = -^^^^div(VfeGi). (3.20) 

Here the quadratic form Q{Y, .) is given by 

g(y, k) ^ ivy(G-ifc • k) ■ {G-^\/yv) - G-^k ■ Vy(G-iVy?? • fc), (3.21) 

and the operator C/o satisfies the estimates 

\gou\lp_^ < ce, \gou\li^^<cEi, (3.22) 

while the linear operator satisfies 

\\g-iu\\r < C£||u||,.-i, for l<r<p — 9, 

< c{e\\u\\s^i+Ei\\u\\o), for s<l-9. (3.23) 

Moreover, operators CJi, t/o and satisfy the symmetry properties 

g^u{±Y*)^g,u*{±Y), J -1,0,-1, u*{Y)^u{Y*). 
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Proof. The proof is given in Appendix [HJ ■ 

The formula for the principal term Gi is a classic result of the theory of 
pseudodifFerential operators |221- The expression for the second term in local 
Riemann coordinates was given in pp . It seems that the general formulae (|3.19|) , 
H3.20|l are new. Note that the ratio det G/J^ = 1 + |Vx??(-'^)P is a scalar 
invariant and the real part of Gq can be rewritten in the invariant form 

J 1 LG-2MF + NF 1 L^ + 2M^^^^ + Ng^ 



Vd^ 2 EG-F^ 2 Eei+ 2F^iC2 + Gf| 

Here we use the standard notations for the second fundamental form L^^f + 
^M^i^2 + N(j and the first fundamental form E(l + 2F^i(2+G£,j of the surface 
0:3 = Tj{X), the vector ^ is connected with the covector k by the relation k = G^. 
The right hand side of H3.24|l is the difference between the mean curvature of 
E and half of the normal curvature of S in the direction ^. Note also that the 
conclusion of Theorem 13.51 holds true without assumption on the smallness of e, 
but the proof becomes more complicated and goes far beyond the scope of the 
paper. 

Lemma 3.6 For m > 4, and U G M^^~^{M.'^ /Tr) with \ \U\\4 small enough, there 
exists a unique function Z € C™^''(]R^) such that 

— = niZ = Z2, (3.25) 

where Hi denotes the average over a period in zi , and Vi are the components of 
the vector field V defined by (|3.1(l . Moreover, Z is even in zi, odd in Z2, 

Z{Z) = Z{Z + (27r, 0)) = Z{Z + (0, 2-k/t)) - 2tt/t = Z{Z + (tt, tt/t)) - tt/t, 

and the shifted function TsZ = Z(- + 6, ■) is solution of the system ^3.25\} where 
V = V{- + (5, •). Moreover, the mapping 

Xl^Zi, X2 ^ Z2 ~ di{zi,Z2) = Z{zi,Z2), (3.26) 

with his inverse 

Zi = Xi, Z2 = X2 + di{xi,X2) 

define automorphisms of the torus M./r^: X i-^ Z = lAi{X) , X = U^^{Z). The 
functions di and di have symmetry (e, o), as above and we have the following 
tame estimates 

||dl||c™-3 + ||di||c..-3 <C„,(||C/||4)||C/|U. 

The automorphism Ui takes integral curves of the vector field V, which coincide 
with the bicharacteristics of the operator Qn — J*{or^J-), onto straight lines 
{z2 — const.}. In other words, bicharacteristics form a foliation of the torus 
with a rotation number equal to 0. 
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Proof. The proof is made in Appendix El ■ 

Let us turn to the proof of Theorem 13. 41 We look for the desired diffeomor- 
phism Y ^ X in the form of the composition {U2 ° ^i)^^, 

where the diffeomorphism Z/^i : M^/r^. R^/F^ is completely defined by Lemma 
13.61 and the diffeomorphism U2 ■ R^/F^ i— > R-^/Fi is unknown. We look for it 
in the form Y = U2{Z) with 

yi^ zi+d2{Z), ?;2 t(z2 + 62(22)), (3.27) 

where functions d2 and 62 will be specified below. Our first task is to make the 
formal change of variables in the left hand side of H3.13|l . We begin with the 
consideration of the second-order differential operator J*{a''^J-). It follows 
from the equality J = V ■ Vx that 

where 

V^V{X{Y)), a^a{X{Y)), B(y) = VyX(r), j(y) = det B(y). 
In particular, we have 

B-i(y(Z)) = VzY{Z) [VzX{Z)Y\ 
On the other hand, Lemma l3 . 61 and formulae H3.27|l imply 

[VzX{Z)Y^V{X{Z)) -Fi(X(Z))ei, VzY{Z)e^^{l + dMZ))eu 
where ei = (1,0). Thus we get 

From this and parametric representation (|3.17l) of the Dirichlet-Neumann oper- 
ators we conclude that the left hand side of the desired identity H3.13|l is equal 
to 

p{Z{Y))dlu{Y)+5{Y))dyMY) + Qiu{Y) + gou{Y) + G^MY), 

where 

P{Z) = ^^—{\ + dMZ)fv^{X{Z))\ (3.28) 
a{X{Z)) 

B{Y)^^^dy,{j{Y)^{Z{Y))]. 
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Note that is a second order operator with respect to the variable yi and 

a first- order operator with respect to variable 2/2- Its principal part is the 
pseudodifferential operator pdyf + Gi which symbol reads 

-pfc? + J"^(adjGfc•fc)^^^ 
which we write in the form 

1/2 



Hence operator H3.13|l can be reduced to the operator with constant coefficients 
at principal derivatives if for some constant u , 

u-^p{Z{Y)) = (rJ(r))~'V5l^- (3.29) 

This equality can be regarded as first order differential equation for functions 
d2{Z), 62(2:2) and a constant v. It becomes clear if we write the right hand 
side as a function of the variable Z. To this end note that since dy-^Z2{Y) = 
d,,Xi{Z) = Q, 

311 - \VzXdy,Z\^ + {Vxri ■ VzX ■ dy.zf = 
[dy,z^{Y)f [1 + (a,,Z)2 + (a^.ry + d,,T^d,,Zf] , 

which along with H3.25|l gives 
On the other hand, we have 

J{YiZ)y^ = det VzY{Z){dct VzX{Z)y^ = T{l+d,,d2){l+d,,e2){d,,Z)-\ 

Substituting these equalities into (|3.29() wc obtain the differential equation for 
d2 and 62, 

dMZ)^[MZ)il + e2{z2))]'^' -1. (3.31) 

where q is equal to 

yi(X(Z))-3(|V^(X(Z))p + (T/(X(Z)).Vxr7(X(Z)))')'^'a(X(Z))9,,(Z(Z))-\ 



(3.32) 

has symmetry (e,e), and is close to /i. Note that it is completely defined by 
equalities H3.1|l . (|3.4|) . and Lemma [3. 61 This equation can be solved as follows. 
We first note that it gives a unique d2 with symmetry (o, e) provided that 

+ 4(z2))-i/2 r qi/2(Z)dzi. (3.33) 
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This now determines a unique periodic odd function 62 provided that 

' — tt/ t 



v = 2t:t I ^ ([ q^''^{Z)dz^ dz2, (3.34) 



which is just when (-0, r\) = 0. We observe that the invariance of q(Z) under 
the shift Tv(,, leads to the invariance of the right hand side of (|3.33|) by the 
change zi > Z2+ t^/t. This means that 62 is indeed tt/t— periodic. It foUows 
from Appendix G (see Lemma G.3) of |2J and Lemma r3.6l that 

||q-/.||c™-4 < c™(||(7||4)||;7|U, 
li^-l^l < C||C/||4, 
IM2||c"- + l|e2||c™-4 < c™(||[/||4)||C/|U. 

From this, the identities 

d{X) = d2{xi,X2+di{X)), e{X)^dx{X)+e2{x2+dx{X))), 

and from Appendix G (see Lemma G.l) of 23' we deduce estimates (|3.15|) for 
the functions d, e, and also estimate (|3.16|) for the Jacobian J. Hence the 
diffeomorphism Y X \s well defined and meets all requirements of assertion 
{iv) of Theorem l3.4l Next set 



k{Y) = (Tj(r))-i VmM, HY) = -Tgn{Y)-'gi2{Y), l+2a = T^gniYy'g,2 

It follows from H3.15|l that 

\\gn - l\\c^-^ + \\g22 - l/r^Wc^-. + Hffiallc-^ < c„(||f/||5)||C/|U, (3.35) 

which leads to estimate (|3.16|) for k and the following estimates for the functions 
a and b 

||a||c™-5 + \\b\\c^-s < c„,(||C/||5)||t/||™. (3.36) 

We are now in a position to define the operators 2t, 03 and £_i defined by 
the left hand side of identity H3.13|l . It follows from H3.28|l and H3.29|l . that its 
coefficients satisfies the inequalities. 

p — i^K, p^^s = dy^{lnKj) 
From this we conclude that 

GriU - J* i^Ju) o {U2oUiy^ = K £ + giu + 1/9^1 (In Kj)Diu + Q5u + i2-iw 
where 

^^-Go, 2-1 ^-G-i, Gi^-Gi~{-Ay/^. (3.37) 

K K K 

By construction, the operator Gi has the following symbol 

{(1 + 2a{Y))kf + 2b{Y)Tkik2 + r'klY^^ - {kf + T^klY'^ , 
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which imphes that Qi = 2loS)i where the zero order pseudodifFerential operator 
2lo has the symbol 

-2i{a{Y)ki + Tb{Y)k2) 

{(1 + 2a{Yj)kf + 2b{Y)rkik2 + T^kl^'^ + {kl + T^^klf'^ ' 

which can be written in the standard form 

1 + (1 + 2aiY))i\ + 2b{Y)iA + if) ' 

Hence the basic identity H3.13|l holds with the operators S, £-i and the operator 
21 — 2to + i^9j/i (In gii) which yields (i). Estimates {ii) follow from formula (|3.38(l . 
estimates (|XT^ and Theorem IX^ 

Note that the functions k, J, a are even and the function h is odd both in 
2/1 and y2. Moreover, they are invariant under the shift Y ^ Y + (tTjTt/t). 
Hence the operator 21 satisfies the symmetry conditions [in). Formulae H3.37|l 
along with Theorem 13.51 imolv that the operators *B and £_i also satisfy the 
conditions (|3.14|l . which completes the proof of Theorem 13.41 

3.3 Main orders of the diffeomorphism and coefficient v 

In using later the Nash-Moser theorem, we need to set 

where C/i^-* is an approximate solution, up to order e^of the system (11.611 . H1.7|l . 
and W is the unknown perturbation. At Theorem 12. 31 for ei = £2 — e/2 (Dia- 
mond waves), we showed how to compute explicitly any order of approximation 
Ue^\ with fjL — ^c{t) = e^/xi(r) + O(e^), ^^l{T) 7^ for r 7^ Tc. Then we need 
to know in particular the principal part of the coefficient v in H3.34|l and H3.9|) . 
which implies the knowledge of the diffeomorphisms Ui and U2- We show the 
following two Lemmas: 

Lemma 3.7 For N > 3, we have 

i^{e,T,W) ^ fi-^ - e^i^iir) + 0{e^), (3.39) 

where fj,c — {I + r^)"-'^/^, and 

^i(t) 3 5 



Mr) = 



1117 1 9 



4/.5 2fii 2/.3 Sfil 4^, 16(2-/1,)' 
which is positive for any r > 0, (j^i(t) > 5/16/ic). 



(3.40) 
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Lemma 3.8 For N > 3, the diffeomorphism of the torus, which allows to 

obtain the form 1^3.13}) for the linear operator Q,^ — Sf*{-Sf{-)), and the principal 
part of f]{Y), are given by 

x = (xi,x2), r = (2/i,y2), x(r)-(Xi(r),X2(r)), 

-'^i(^) = ?yi + |sinyiC0S?/2 + £^(Ciisin2?;i +^i2sin2yiCOs2y2) + O(e^), 

X2{Y) = i?;2 + e-r cos yi sin y2+e^(C2i sin 27/2 +62 cos 2yi sin 27/2) +0(£^), 
r 

fj{Y) = cosyi cosy2 + e^{??oo + ??oi cos2j/i + 7702 cos2y2(l + cos2yi)} + 0(e^ 

Mc 

where e > is defined by the main order e^Q of [/i^"* and 

fi = ^c + e'Mi(T) + 0(e4), Mc = (1 + t2)-i/2^ 



and where 

61 



1(3 3 5 3 1 



1 ^ _ T 17 

62 - 62 



11 1 
'700 - -7-^ - a ; ^702 — Z 

1 1 1 



'701 



4/^, 4^.2 4/^3 8(2 -^e)- 



Proof: see Appendix IdI 



4 Small divisors. Estimate of £— resolvent 

It follows from Thcorcm l3.4l that the pseudodifFerential operator £ = Q^'^'' +^^yi 
with the symbol 

L{k) = -iykl + {kl + T'^k^y/'^, k€l?. (4.1) 

forms the principal part of the linear problem. In this section we describe the 
structure of the spectrum of £ and investigate in details the dependence of its 
resolvent on parameters v and r. The first result in this direction is the following 
theorem which constitutes generic properties of £. 

Theorem 4.1 For every positive v and t, the operator £ is self adjoint in 
FL^i^jV) and has the natural domain of definition 

£»(£) = {u e i7°(MVr) : ^ L{kf\u{k)\^ < 00}. 
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The spectrum of £ coincides with the closure of the discrete spectrum which 
consists of all numbers Afe = L{k), k € N^; the corresponding eigenfunctions are 
defined by 

(27r)-^ exp(±ifc • Y), (27r)-i exp(±fc* • Y) where k* = (-fci, fca). 

For each real >c =^ \k, k & I? , the resolvent (£ — is a self adjoint unbounded 
operator defined in terms of the Fourier transform by the formula 

(£ ^^-^u{k) = {L{k) - x)-^u{k), k G Z^. (4.2) 

Proof. The operator £ is unitary equivalent to a multiplication operator in 
the Hilbert space I2, which can be represented by the infinite diagonal matrix 
with diagonal elements L{k), k G I?. It remains to note that for such a matrix 
the spectrum coincides with the closure of diagonal, and the discrete spectrum 
coincides with the diagonal. ■ 

Note that zero is nontrivial eigenvalue of £ if and only if the dispersion 
equation L{k) = has a nontrivial solution k = k^. The number of such 
solutions depends on the arithmetic nature of parameters v, r. We restrict our 
attention to the simplest case when fco = (1, 1). It is easy to see that in this 
case the point {u,t) = {vojt) belongs to the positive branch of the hyperbola 

Our aim is to investigate in details the dependence of £-resolvent on param- 
eter ly. With further applications to the Nash-Moser theory in mind wc take the 
perturbed values of parameter ly, and a spectral parameter x in the form 

i'j{s) = vq — s'^vi + e^Pj{e), Xj{e) = e^kj{e), (4.3) 

where = ^0(1") and z/i = z^i(t). Here functions i)j and kj, j > 1, are defined 
on the segment [0, ro] and satisfy the inequalities 

\i>jie)\ + \^j{e)\<R, 
\i>j{e') - i>,{e")\ + - ^(e")! < R\e' - e"\, (4.4) 

\i>j+i-i>j\ + \kj+i- k\<R{2-^). 

Denoting by £0 the operator £ for ly = lyoir), the next theorem establishes 
the basic estimates for a resolvent {£, — on the orthogonal complement 

to ker £0, which are stable with respect to perturbations of parameters from a 
suitable Cantor set. 

Theorem 4.2 (a) For each a e (0, 1], there is a set of full measure Ot(, C 
(1,00) so that whenever uq g 91^ and v = vq{t) = (1 + r^)^/^, zero is a 
non-trivial eigenvalue o/£o and 

ker £0 = span |l,exp(±2fco • Y), exp(±fcQ • i^)}, 
||£o 'u|U-(i+a)/2 < c||m|1„ when u& {ker ilo)^nH'{Ryr). (4.5) 
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(h) Suppose that vq G with a G (0, 1/78), then there is a set £ C [0,ro) so 
that 

-J / ede ^ 1 as r ^ 0, (4.6) 



[0,r]n£ 

and for all e G £, j, s > 1 and v — Vj{s), k — Xj{s), 

\\{2.- Hy^u\\^_^<c\\u\\s when ue{ker £o)^r\H''{R^/r), (4.7) 

where the positive constant c depends on a, and R only. 

Remark 4.3 The fact proved at Lemma \3. ?| that i'i{t) > 0, allows the freedom 
to take T as a function of e as t — tq + e'^Ti. Then, for |ti| small enough, we 
have 

LOi := uj(){vQ^yi - T^^Ti) ^ 
and TheoremTT^ still holds true. 



It follows from formula (|4.2I) and the Parseval identity that Theorem l4.2l will 
be proved once we prove the following 

Theorem 4.4 (a) For any a G (0,1] there is a set of full measure OIq. C 
(1, oo) so that for all G D^q and v = vt^ir), 

|L(fc)| > c|fc|-(i+")/2 when k^^,±kQ,±kl. (4.8) 

(h) Suppose that vq G with a G (0, 1/78), then, there is a set £ C [0, tq) 
satisfying (|4.6|l so that for all s <E £ , j > 1 and v = Vjie), h = >Cj[e), 

\{L{k)~ >i)-^\>c\k\-^ for all fc 7^ 0, ±fco, ±fco: (4-9) 

where the positive constant c depends on vq, a and R only. 

The proof naturally falls into four steps. 

First step. We begin with proving two auxiliary lemmas which establish a 
connection between the symbol L and the Diophantinc function defined by the 
formula 

m C 

{m,n)^uj (m,n) gN^, (4.10) 

where constants uj, C are connected with parameters v, r and x by the relations 

lo = v/t, C =1/{2vt) - x/t. 

Lemma 4.5 Let for some g > 2 and a G [0, 1], parameters v , t and a vector 
k satisfy the inequalities 

Q<Q-^ <V,T<Q, \X\<Q, |/Cl|>2^^^ (4.11) 

\^kl ~ \k2\ ~C\> 5{g' + ei")|A;iri-", (4.12) 
Then \L{k)-x\ > |/c|-(i+")/2. 
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Proof. Since L{k) is even in k, it suffices to prove the lemma for 

k = (n, m) witli integers n > 2g^ , m > 0. 
We begin with the observation that for n'^ > 2g^m, 

^L{k) + >c > g^^fi^ — \/ Q^irfl + n? ^ g > 



.g 2^2 ' g,4/ y ^ - 3 = 
Hence it suffices to prove the lemma for 

< 2g'^m. (4.13) 
Since for all non-negative cr, < 1 + a/2 — (1 + aY^"^ < (7^/4, we have 

-T-^L{k) =un^ -m- {2T'^m)-^n'^ + {ir^m^y^n'^ Oi, Oi € [0, 1], 
which along with the identity 

(2T^my^n'^ = (2i^t)"^ + (2i/Tm)"^ (wn^ - m) 

yields 

-T-^e{k){L{k)- k) ^ ojr? - m - C + 02- (4.14) 

Here 

02 = C(l - e{k)) + e(fc)(4r''m3)-in4oi, e(fc) = (l - {2Tvm)-^y^ . 

Since 

(2:/Tm)"^ < 2"^(?2m"^ < p^n"^ < 2-\ 

we have 

1 < e{k) < 2, e(fc) - 1 < 2/n-^ 
which along with the inequality \C\ < 3g^ /2 leads to the estimate 



\02\<hg' + ^^)^<5g'"n 
V 2 m'^ / n-^ 

From this, (|4.12l) and H4.14|l we conclude that 

2g\L{k) - x\> \ujr? -~ m - C\ - 2e^"n"2 > 

It remains to note that due H4.13|l the right hand side is larger than |fc|"'^^+"^/2 
and the lemma follows. ■ 

Lemma 4.6 Suppose that parameters vq, tq and v, t, >f satisfy the following con- 
ditions. 



10^-2 
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(i) There is g >2 so that 

< g^^ < V, T, < g, \x\ < g. 

(a) The dispersion equation 

vokl - ^kf+T^k^ = (4.15) 

has the unique positive solution k = ko. 
(Hi) There are positive N and q so that 

\u}n^ — m — C\ > qn^^^"^ for all integers n > N, m > 

Then there are positive S and 70, depending on N , q and vo only, so that 
|L(fc)->f| >7o|fcr(^+")/2 /or an fc 7^ 0, i/co, ±fc* (4.16) 

when 

\v -va\ + \x\<5. (4.17) 

Proof. It follows from [Hi) that v, r and k meet al requirements of Lemma 
01 which implies that \L{k) - h\> foj. \k^\ > On the other 

hand, since 

lim |i(fc) — >^| > lim \j Q^'^k\ ^ k\~ Qk\ — g ^ 00, 

there is N* depending on g and N only so that \L{k) — >c\> for all 

\k\ > TV*. Since the number of wave vectors k in the circle |fc| < iV is finite, the 
existence 6 and 70 follows from continuity of L{k) as a function of parameters 
z/, T. ■ 

Second step. Next we prove that condition (Hi) from the previous lemma is 
the generic property of function H4.1()(l which leads to assertion (aj of Theorem 
14.41 In order to formulate the results let us introduce the important function 
d : I— > M defined by the formula 

TTi C 

d(m,n)^uja -, where loq — I'a/T^Ca ^ (2i^ot)^^ . (4.18) 

Lemma 4.7 For each a G (0, 1] and q > there is a set of a full measure DJla 
in (1,00) so that for all vq G OJIq, there exists N > 0, depending on vq, a and 
q only, such that 

\d{m, n)\ > qn^^~°' when n> N,m>0. (4.19) 
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Proof. Noting that Co = {clJo ^ '^)/2, we rewrite inequahty (|4.19|l in the 
equivalent form 

IdoiuJoTiTiTn)] > qn~^~°' where do{uJo,m,n) = ujq j(m. — 2~^ujo + 2~^ujQ^y 

Without loss of generality we can assume also that loq G [1, a], where a is an ar- 
bitrary positive number. The set of points loq for which |c?Q(a;o, ti, n)| < qn~^~" 
can be covered by the system of the intervals t(rn, n) labelled by integers m and 

n. It is easy to see that their extremities tf^(m, n) satisfy ^n^ + i^w^(m, n) > 

1^ — . On the other hand, since 

d^do{ujQ,m, n) = l + (2n^)-i + (2t^2^2)-i ^ 2], 

the length of each intervals is less than 2qn^^^°' . Hence for fixed n, the number 
of intervals L(m,n) having nonempty intersections with a segment [I, a] is less 
than a(rt^ + 1/2) + q/-n}^°' . From this we conclude that 

^ ^ meas ^(to, n) < 2q{-a + g) ^ n^^"" < cq{a + q)N-°'. 

n>Af m:t(m,n)n[l,a]^0 n>N 

Hence 

meas {wo G [1, a] : |do(w, m, n)| > qn^'^^" for all m>0,n>iV}> 
a — C(7(a + (7)iV^" ^ a as — > 00. 

It remains to note that the mapping i^q luq takes diffeomorphically the interval 
[1, 00) onto itself and the lemma follows. ■ 

Next lemma shows that almost each point of Tla satisfies the following 

Condition M. There are positive constant cq, ci such that for all integers 
n > 0, m > 0, 

\d{m,n)\>con-^-°', (4.20) 
le^^'"""- iP^ < cm^ (4.21) 



Lemma 4.8 For each a e (0, 1) there is a set C 97la such that meas (dJIaX 
*yict) = and Condition M holds for each vq £ OIq.. 

Proof. We begin with the observation that inequality H4.21(l holds for almost 
every loq with a constant ci depending on loq only. Hence it suffices to show that 
inequality H4.20|l holds true for each G 9TIq,. Fix q = 1, £ Tla and note that 
by lemma (|4.7|l . |(i(m, n)| > n"'^^" for all n > N. Since lim d{m,n) — — cxd, 

m — ^CJO 

we can choose N such that this inequality is fulfilled for all (m, n) outside of 
the simplex n > 0, m > 0, n + m < A'^. It remains to note that this simplex 
contains only finite number of integer points and d{m, n) for irrational wq 
different from the sum of a rational and the square root of a rational. ■ 
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Third step We intend now to study the robustness of estimate (|4.2U|I when 
one adds a small perturbation to d(m, n). In order to formulate the correspond- 
ing result we introduce some notations. For each e € [0,ro] and j > I set 

Cu,{e) = v,{e)It, C,{e) = {2v,{e)TY^ - ^,{e)T-^ 
It follows from (|4.3|) that they have the representation 

tJj(e) = t^o - e^t^i + e^i^j(e),C'j(e) = Co - e^'^j{e) e e [0,ro], 

in which u\ = t^^vi. Our task is to obtain the estimates for the function 
Dj : [0, ro] x ^ K defined by 

b{s, m,n,)= cJ, {e)-^-^+e'^. (4.22) 

For technical reason it is convenient to formulate the problem in terms of a new 
small parameter A. Set 

Wj(A) ^Luo- Xuji + A3/217j(A), a e [0,po], (4.23) 
where po = Tq. It follows from l|4.4() that for suitable choice R and ro, 

|%| + |^,| < R, + < R{2-^), m.(A)| + |^;.(A)| < 

(4.24) 

Set 



D,{\,m,n,) -.^ D^{VX,m,n) ^ ojjiX) - ^ - ^ + X^^. (4.25) 

Definition 4.9 For positive N, q, r denote hy'Hj{N,q,r) the set defined by 

T-Lj{N ,q,r) = |a G (0,7') : \Dj{X,m,n)\ > qn"'^ for all integers n> N, m > o|. 

(4.26) 

Theorem 4.10 Assume that G OIq, with a G (0, 1/78). Then, for each q > 0, 
there is N > such that 

- meas f^{nj{N,q,r)} ^ 1 as r 0, (4.27) 

and there exists c* such that 
1 



— 771 eas 
r 



f]{njiN,q,r)}>l-c*r'^ for < r < r* , (4.28) 



where <w < min{a/(3 - a), (78"^ - a)/ {3 + a)}. 
Proof. Subsection 14. II is devoted to the proof. ■ 
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Fourth step. We are now in a position to complete the proof of Theorem 14. 41 
First note that by Lemmas 14 . 71 and |4 . 81 for each i^q £ O^q, parameters v = i'o(t) 
and X — meet ah requirements of Lemma 14.61 with S — 0. Applying this 
lemma we obtain (|4.8|l . It remains to note that for almost all vq the dispersion 
equation has the only positive solution fcg and assertion (a) follows. In order to 
prove (6) choose an arbitrary i/q e D^q, q > and set 

£ = {e>0: e'e f]H,iN,q,r)} , 

j>i 

where N is given by Theorem 14.101 It follows from this theorem that the set 
£ satisfies density condition 1)4. 6() . On the other hand, (|4.26() implies that for 
e € £, parameters v = J^j(e), >f = ^ji^) meet all requirements of Lemma 14.61 
Since 

(e) — i^o\ + \Hj{e)\ as e — > 

uniformly with respect to j, there exists £2 > depending on i^q and R only 
such that for all e G £ n (0,e2), parameters — Vj{e) and x = >Cj{e) satisfy 
inequality H4.17II which yields (|4.16() (where a = 1), and the theorem follows. 

4.1 Proof of Theorem ICTTl 

Our approach is based on standard methods of metric theory of Diophantine 
approximations and Weil Theorem on the uniform distribution of a sequence 
{ojon^}. Without loss of generality we can assume that wi > 0, and po = Vq 
satisfies the inequality 

Po < ^(^1 + ^ '^(^) ^ ^ fo^' ^ ^ Po- (4-29) 

It follows from (|4.24|) that the sequences flj a.nd (pj converge uniformly on [0, pa] 
to functions fioo and 9300 such that 

|noo|+|(^oc| < R, \noo~nj\+\ipoo-iPj\ < 2^-^R, \nu>^)\+\ip'^ix)\ < 2-^\~^/^R. 

(4.30) 

Our task is to estimate the measure of intersection of the sets Tij. We begin 
with the investigation of their structure. 

Structure of a set 'Hj{N,q,r). The main result of this paragraph is the 
following covering lemma. Fix an arbitrary positive q and set 

(4.31) 

where cq is the constant from Condition M . 

Lemma 4.11 For N > N3{q), r < r2 and 1 < j < 00, the set [0,r]\Hj{N,q,r) 
is covered by the system of the intervals 

I-i{m,n) = {Xj{m,n),X'^{m,n)), n>N /,-(m, 71) n [0, r] ^ 0, 
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such that 

(i) X^{m,n) are solutions of the equations 

Xf{., - ^) - iXff/^niXf) = dim,n) ± ± (4.32) 

with d{m,n) > 0. They satisfy the inequalities 

< Aj(m,n) < A+(m,n) < 2r2, (4.33) 

2 2 
- — d(m,n) < Xf < —d(m,n) (4.34) 
bull t^i 



i^^<A;(™,n)-A7(.n,.)<^-L. (4.35) 



(a) For a fixed n> N, the left extremities X''_{m,n) strongly decreases in m, 



A-(Mj-„(r),n) < A7(Mj,„(r) - l,n) < ... < A" (m^- „(r), n), (4.36) 

2 

where 



X-{k-l,n)-Xj{k,n)>,^—^, (4.37) 



?Tij,„(r) ~ min{m > : Ij{m, n) n [0, r] ^0}, 
Mj_„{r) = max{m > : /j (m, n) n [0, r] ^0}. 

('iiij For each such interval with I.j{m,n) H [0,r] 7^ , 

5 / 2r \ ^/(^+"' 

n) < -wir, n > — , A+ (m, n) < 2r. (4.38) 

2 \ 5tJi / 

(iv) If intervals Ij{m,n) and loo have nonempty intersections with (0,r], then 

|A±(m,n)-A±(m,n)| <2--'"-V. (4.39) 
Proof. By abuse of notations, we suppress the index j. Proof of (i). 

-d^D{x, m, n) = uji- ^P±^ _ lx'/^n{x) - x'/^n'{x), 

2 

it follows from l|OT|) and that for A < 2r2 

27? 2_/? 
wi ^ - 2i?Ai/2 < -dxD{X, m, n) < wi + — + 2i?Ai/^ 

which along with H4.31|l implies the inequalities 

Y < ^dxD{X,m,n) < ^ for A E [0,2r2], n > N3. (4.40) 
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Hence the function — i?(A, m, n) is strongly monotone on the interval (0, 2r2). 
Therefore for r < r2,N > N3 the set (0, r) \ Ti{N,q,r) can be covered by the 
system of the intervals 

J{m,n) = (/3(m, n), 7(m, n)) , < (3{m,n) < 7(771, n) < r, 

such that 

-D(l3(m,n),m,n) = --^ if /3(m,n)>0, (4-41) 

—D{j{m,n),m,n) — +-^ if j{m,n)<r, (4.42) 

--^ < -L»(/3(m,n),m,n) < if /3(m, n) = (4.43) 

— ^ < —D{'-f{m,n),m,n) < if 7(771,71) = r, (4.44) 

77** 77 

Note that, by condition M, 

Q q Co q 

~ 5^ — < 1 "('77, 77) , 

which along with 1)4.3111 yields the inequalities 

i|d(777, 77)1 < |d(777, n) ± -^1 < ||d(777 , 77|, |d(777, Ti) | > 2-^ for n>N3. 

(4.45) 

From this and the equality D{0,m,n) = d{m,n) we conclude that case (|4.43ll 
is impossible and (3{m,n) > for all intervals J{m,n). Let us consider case 
1)4.44)1 . Since the function —D{X,m,n) increases in A on the segment [0, 27-2], 
there is a maximal 7* in (0, 2r2] such that 

{(3{m,n),r) C (/3(r(7, 77), 7*) C (/3(777, 77), 2r2), 
|-D(A,TO,77)| < g77~* in (/3(77i, 77), 7*). 

Let us show that 7* < 2r2. If the assertion is false, then r2,2r2 G (/3(7J7, 77), 7*] 
which yields 

~qn~^ < —D{r2,m,n) < —D(2r2,m,n) < qn~'^. 

Thus we get 

r2 min {—dxD{X,m,n)} < D{r2,m,n) — D{2r2,m,n) < 2qn~'^ < 2qN^^. 

[r2,2r2] 

From this and ()4.4(J)) we obtain the inequality 



!1 ^ o„Ar-4 

2 



which contradicts l)4.31)l . and the assertion follows. In particular, we have 
~D{'y* ,m,n) = qn^'^. Since the equations —D{X±,m,n) = zLqn~'^ are equiv- 
alent to ()4.32)l . in the cases ()4.41)1 . ()4.42)l we have in the cases ()4.41ll . ()4.42)l 
that 

f3{m,n) = X~ {111,71), 7(771, 77) = A^(7r7, 77), J{m,n) ~ I{m,n), 
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in the case (14 .441) we have 

l3{m,n) — X" {m,n), j* = {m,n), J{m,n) C I{m,n). 

Hence the set [0, r] \ Ti{N, q, r) is covered by the intervals /(to, n) with the ends 
satisfying (g^lll, 1233 • K follows from llOT|l and that for n > N^^q), 

A < r2, 

<x(oj,-^)- X'/'n{X) < -Lu.X. (4.46) 
4 \ IT? J 4 

In particular, the left side of equation H4.32|l is positive, which along with H4.45|l 

yields positivity of d(rn,n) in the right side of equations (|4.32|) . Combining 

gives Next equation implies 

D{X~ {m, n), TO, n) — Z?(A^(to, n), to, n) = 2qn~'*. 

From this and (|4.4U|I we conclude that 

a;i(A^(TO, n) — A^(m, n)) < 4:qn~^, 



which yields 



a;i(A+(m, n) — A {m,n))>—qn 

o 



4q I , , ^ , ^ ^ 4g 1 
J < A"^(TO,n) - A_(to, n) < j 



and the assertion follows. 

Let us turn to the proof of {ii). By Lemma [4.11^ 1. I{m,n) n [0,r] 7^ if 
and only if A~(to,7i) < r. In particular, < A^(TO„(r), 71), A~ (M„(r), 71) < r. 
On the other hand, A^ (to, 71) is a solution to the equation 

-D{X'' (to, 71), TO, 71) = -qn^'^. (4.47) 

By (|4.40() we have 

-dmD{X,m,n) = ^ - ^^>^^i^^^''^) - ~Y' ^'^^ AG[0,r]. 

Hence for fixed n > N3, the implicit function A~(to,7i) satisfying the equation 
(|4.47|l is defined and strongly decreases on the interval [mn{r) , Mn{r)] which 
yields (|4.36(l . The relation I{k, n)n [0, 7-] 7^ for each integer k e [TO„(r), Mn{r)] 
follows from (|4.36|l . Next if m„(r) < fc — 1 < fc < M„(r), then we have 

= -D{X^{k - l,n),k- l,n) + D{X- {k,n),k,n) = 
= ~D{X^{k - l,n),k~ 1,71) + L>(A"(fc,n),fc- l,7i) + 
-D{X~{k,n),k~- l,7i) + D{X~{k,n),k,n) 

= -D(X^(k - l,n),k- l,n) + D(X-{k,n),k- l,7i) - ^ 

71^ 

< max {~d\D(X, fc — 1, 7i)}(A^(fc — 1, 71) — A^ (fc, 71)) 

AG[0,r] 71^ 



<^{X-{k-l,n)-X-{k,n)) ^, 



44 



and hence 

which completes the proof of {ii). 

Next note that is a consequence of the inequahties 

2co 1 2 ^ 

and (iii) follows. 

It remains to prove (w). We have 

-Doo(A^,m,n) - Doo{\f,m,n) = Dj{\f,m,n) ~ Doo{\f ,m,n). 

On the other hand, 

\D,{Xf,m,n) - Do.{Xf,m,n)\ < - + {Xf^/^in^o - ^1 < 

2r2R2-^{N-^ + 1) < 4r2i?2"^ 
which along with H4.4()|l and (|4.31(l implies 

|A±(m,n)-A±(m,n)| < ^r22-^ < ^2-^ 
■' uJi 16 

and the lemma follows. ■ 

Cardinality of the set {Ij{m,n)}. Set 

On = {^^Qi^^ ^ C} = decimal part of ujon^ — C, 

and introduce the sequence of numbers Sn (r) defined by 

5 

Sn{r) = 1 when 9n < —ujirri^, and (5„(r) = otherwise. 
Lemma 4.12 For each n > and < r < r2, 1 < J < 00, 

card {m : Ij{m, n) n [0, r] 7^ 0} < (5„ + crn^. (4.48) 

Proof. For simplicity we omit the index j. Denote by I the totality of all inter- 
vals I{m,n) given by Lemma f4. 1 II such that I{m,n) n [0, ^ 0. By assertion 
(ii) of Lemma 14.111 there is one-to-one correspondence between the intervals 
I{m,n) S X, having nonempty intersection with [0, r], and the sequence of 
A~(fc,n) given by H4.36|l . In particular, we have 

card {m : I{m, n) n [0, r] 7^ 0} = Af„(r) - m„(r) + 1. (4.49) 

On the other hand, inequality (|4.37(l yields 

A-(TO„(r),n) - A-(M„(r),n) > (M„(r) - m„(r))-^. 
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Since mn{r) , Mn{r) G [0,r], we conclude from this that 

M„(r) — m„(r) < cri^r. 
Combining this inequahty with H4.49|l we obtain 

card {m : I{m, n) n [0, r] 7^ 0} < 1 + crn^ (4.50) 
On the other hand, for given n > N3, 

( I J /(m,n) ) n [0,r] = if min X^(m,n) > r. 

\ ^ I m:l{m,n)el 

Since, by Lemma [4. Ill < ■^d{m, n) < A_(m, n), we conchide from this that 
( I J /(to, 71)) n [0,r] = if min d(m,n) > -uJir 

V ^ / m:d(m.,n)>0 2 

m:7(m,n)tET 

Obviously 

1 



min c?(to, n) — — min {ujon — C — to) = — 

m:d{m,n)>0 Tl-^ m:uJan^ -C -m>0 Tl 



which yields 

^ 5 
n2 2 



card {to : /(to, n) n [0, r] 7^ 0} = for > -uj^r (4.51) 



Combining (|4.5U|) and H4.51|l we obtain 

1 card {m : /(to, n) n [0, r] 7^ 0} < (5„(r) 

m:/(m,n)n[O,r]^0 



and the lemma follows. ■ 

Proof of Theorem I4.10L First note that 

meas ([0,r]\P| nj{N,q,r)^ < ^ meas |J Ij{m,n)j p(r). 

j=l (m,„):„>JV3 J = l 

n [0,r] ^ 

The following lemma gives an estimate for p in terms of the sequence <5„(r). 
Lemma 4.13 For each N > N3, <r <r2, and a G (0, 1/(3 + a)), 



p(r)< y c^(ln?i + l) + cri+'^, (4.52) 

n>C3r-i/(3+Q) 

where c depends on uji, a, a and R only, — i 5^ ) > co constant 

from Condition M . 
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= —Inn- —In [—29], n>l. 



Proof. Introduce the sequence 

ln2 In2 '^^Swf 

It is easy to see that <x'i2^^ * < Aq/ (SoJin*) for all j > <^(n). It follows from this 
and l|4.39(l that for all such j, the intervals Ij{m, n) have nonempty intersections 
with Ioo{m,n) and Ujy^(^n)Ij{m,n) C /ooC?^,'^), which yields 

20q 1 
3^1 



y Ij{m,n)j < meas /oo(to, n) < ^-j, 
where Iao('m,n) is the 4q/(3ti;in'*)- neighborhood of /oo(w, n). Thus we get 



me as 

j><;{n) 



pW < H 51 X! meas/j(m,n)+ 

" > «3 j = l (m:/j(m,n)n[O,r]^0 

n [0, r] ^ J\ ■ ^ I 

^ ^ meas ( |J lj{m,n) 

" > N3 (m:l^lm.n)n\0,r]^(tl 7><?(n) 

n [0. r] ^ ' V - y L J 

which leads to 

1 > «3 i=l (m:/,(m,n)n[O,r]^0 

E 1' 

r!,>A''3 (m:/oo(m,n)n[O,r]^0 

Since, by Lemma l4.11[ the intersection Ijim, n)n[0, r] is empty for n < c3r^^/'^'^+"\ 
we have 

c ^^"^ 

^W- E ^E^^^'^{'^--^j(™'")^[O,r]^0} + 
^ — ' 

n>C3r-i/(3+a) 



which along with H4.48|l yields 



- E ^?(n)(<5„(r) +crn2) + E ^(f^^ + crn^), 

n>C3r-i/(3+a) „>c3r-i/(3+a) 



(5n(?') , , ^ , 1 + Inn 



<c E ^(l + lnn)+- E 

n>C3r-i/(3+a) „>c3r-l/(3+Q) 

n^ 



< c^(l + Inn) +cri 

n>C3r-i/(3+Q) 



47 



which completes the proof. ■ 
Introduce the sequence 

5 l-Q 

i^n{r) — 1 when 9n < —ujir^-" , and i^„(r) = otherwise. 

Lemma 4.14 Let < r < r2and < /3 < a/(3 — a). 
Ifcsr^+o >r3~°, then 

p(r) < cr^+f^. (4.53) 



p(r)<c ^ + Inn) + cri+'5, (4.54) 



//car 3+° <r3-°, t/ien 

l'n(r) 

C3r-i/(3+<>)<„<r-i/(3-a) 

where c does not depend on r. 

Proof. Since /3 < 1/(3 + a), we have from inequality (|4.52|) 

p(r) < cn(r) + y c^^(l + Inn) + cr^+f^ , (4.55) 

n>r-i/(3-Q) 

where 

n(r) = " . (1 + Inn) for C3r~ < r~ , 

^ — ' n* 

C3r-i/(3+<>)<„<r-i/(3-a) 

and n(r) = otherwise. It is easy to see that 

J2 ^(l + lnn)< ■^{l + lnn)<cr^+l' (4.56) 

n>r-i/(3-Q) „>r-i/(3-a) 

If n < r-i/(3-") thcn^,, > ftJir^i-^^/^^-a) yields 6'„ > ftJin^r. In other words, 
equality fc'„(r) = yields (5„(r) = 0. Hence 

n(r) < ^^^(l + lnn) for car^ < , (4.57) 

C3r-i/(3+o)<„<r-i/(3-a) 

and n(r) = otherwise. Substituting H4.56|l and 14.57|l into (|4.55l) and noting 
that we obtain needed inequalities (|4.53|) . 14.54|l . ■ 

Now set 



78(3+o) ^ _ ^ ^ 



i = 



3 - a 78(3 + a) ' 



Remark This complicated formulae simply express the fact that the number 
p = (^c^^r^^^^'^"^^^^'^^ satisfies the inequalities 

^wir(i"")/(3-a) < p < 1/4 for r< r^. 
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Lemma 4.15 Let , < r < rs, c3r-V(3+«) < r-V(3-«)^ and < 7 < (1/78 - 
a)/(3 + a). Then 

J2 ^(i + inn)<cri+^. (4.58) 

Proof. Let p and q be the minimal and maximal integers from the interval 

(c3r~^/^^+"\r~^/(^~")). Introduce the average 



1 " 



n 



Noting that ^'„(r) = nSn — {n— we obtain 

C3r-i/(3+a)<„<r-i/(3-c,) P<n<q 

^ /1 + lnn l + ln(n + l)\ ^^1^ 
^ L (^^4 (n + 1)^ + 



p<n<q—l 

Since 9 > c3r~-^/(^+'*^ and 

'1 + lnn l + ln(n+ 1)\ ^ c , , 

^7TI)^j"-^^' + '""^' 

we conclude from this that 

2. ^(1 + 1^") 



C3 



r-l/(3+a)<„<r-l/(3-a) 



<c{ ^ J_(i + inn) + r3/(3+")} sup 5„ < 

C3r-l/(3+c)<„<r-l/(3-o) „>e3r-l/(3+a) 

cr(3-^)/(3+«) sup S„, (4.59) 



n>C3r-i/(3+a) 

where 1? is an arbitrary positive number. Now set 

1/78 



P = 



(c3-Vl/(3+")j 



It follows from the choice of rs and the inequality r <rs that 

^^^r(i-«)/(3-a) < p < 1/4^ „ > c3r-V(3+«) ^ n > (4.60) 
It follows from this that 

n ^-^ n ^-^ 

1 < k < n, 1 < k < n, 

«n < |<.'ir(l-°)/(3-») en€[0,p] 
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Applying Proposition IE . 1 1 we obtain that for any n > c^r — p 

5„ <cp = cri/78(3+")^ 

Combining this inequahty with l|4.59|) we finaUy obtain 

y + lnn)< cri/7«(3+a)^(3-^)/(3+a) ^ ^^1+7^ 

C3r-l/(3 + a)<„<r-l/(3-<») 

and the lemma follows. ■ 

Finally, combining inequalities 1)4.53(1 . H4.54|l and H4.58|l we conclude that for 

< r < rs, 

p(r) < cr{r^ +r''') < cr^^"^ , 
which completes the proof of Theorem 14.101 

5 Descent method-Inversion of the Unearized op- 
erator 

In this section we give a general method of reduction, the descent method, 
allowing to transform the original linear operator of order 2 into the sum of a 
main operator with constant coefficients and a smoothing perturbation operator. 
Let us consider the basic operator equation 

£M + 20iu + Q5u + £_iit = /, (5.1) 

with zero-order pseudodifferential operators 21, *B and an intcgro-differential 
operator £_i of order —1. Assume that they satisfy the following conditions. 
Symmetry condition: 

(i) For all y e T2 and \^\ < 1, 



A{Y,^O^A{Y,0, (5.2) 

which means in particular that Stw is real for real- valued functions u. 

(ii) Equation 1(5. 1() is invariant with respect to the symmetry Y ^ Y* ^ 
(— 2/1, 2/2). This is equivalent to the equivariant property 



2lDiu(r*) = 2lI)iM*(r), ^u{Y*) = <Bu*(y), 
£_iu(r*) = £_iu*(r), u*iY)=uiY*), 



(5.3) 



which can be also written in the form 

A{Y*,C)=-A{Y,^), B{Y*,C) =B{Y,^). (5.4) 

(Hi) For each s e [1,^ — 3], ^{M.'^ /T) is invariant subspace of operators 
2tDi, «B and £_i. 
Metric condition: 
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(iv) There are exponents r, s, and Z, satisfying inequalities 

1 < r < s < I - 10, 

and e E (0, 1] so that 

mi^r + \'&kr<e, |2l|4,i + I^U,; < (». (5.5) 
For each s G [r,l — 10] there is a constant Cg so that 

||£_iu||r < c£||u||r-i, (5.6) 
\\2-iu\\s < c{e\\u\\s-i + Cs\\u\\o)). 

Restrictions on the spectrum and resolvent of £: 

(v) The selfadjoint operator £ : H^^iR^/T) ^ Hl^{R^/r) has a sim- 
ple eigenvalue vie"^ + 0(£^) with corresponding eigenfunction ip'^^^; the space 
iJg g(M^/r) is the sum of orthogonal subspaces 

<,(MVr)= span{^(0)}©<'i-; 

H^'^ are invariant subspaces of the operator £. Denote by Q the orthogonal 
projector of H^^^{R^/r) onto H^'^. 

(vi) For 

^ y" (^A{Y,0,lf - AiyB{Y,0,l))dY, (5.7) 



"-16. 



T2 



and all s > 1, the inverse (£ — ><) ^ : iJo^g^ ^o,e^'^ is continuous and 

||(£-><)-iu||,_i < c(s)||u|U. (5.8) 

Nondegeneracy condition: 

(vii) 

x=0{e^), J ((£ + io)i^("' ~iiQ(£-><)-iQi3?/'^°^)i^'"'rfr = @e2 + 0(e3)^ 

(5.9) 

where @ is a non-zero absolute constant, and the operator — 212) i -|- + £-i . 

Here and below we denote by c generic constants depending on r and I only, 
and use the standard notation 0(e") for quantities which absolute value does 
not exceed ce". The following theorem is the main result of this section. 

Theorem 5.1 Under the above assumptions, there is a positive constant Eq 
depending on l,r only so that for any f e ^(R'^ /T) and e < Eq , equation 
H5.1|l has a unique solution u G H^'^^{R^^ /T) satisfying the inequalities 

<e~'c||/||,, (5.10) 
\\u\\s-i < £"'c||/||,(l + Cs + |2tk,+io + l^k.^+io) + c||/|U. (5.11) 
Proof. The proof constitutes the next two subsections. ■ 
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5.1 Descent method 

In this subsection we develop an algebraic method which allows us to reduce 
H5.1|l to a Fredholm -type equation. The main result in this direction is the 
following 

Theorem 5.2 Let zero- order pseudodijferential operators 2t and S satisfy sym- 
metry conditions (i), (ii), metric condition (iv), and k is given by (|5.7() . Then 
there exist integro- differential operators (L, 2;, ^ satisfying symmetry condition 
H5.3|l with 03 replaced by £, (E, ^ so that 

(£ + 2ii)i + »)(i + e:)ni = + + (5.12) 

and for 1 < s < I — 10, 

llCulU + < ce||u||, + c(|2l|4,.+9 + |5B|4,.+9)||w||o, 

, . (o.loj 

\\du\\s < ce||u||s-i + c(|2l|4,s+io + |S|4,s+io)||w||o- 

The proof falls into three steps and is based on the following proposition, 
which gives the special decomposition of zero-order operators and plays the key 
role in our analysis. In order to formulate it we introduce the important notion 
of an elementary operator. 

Definition 5.3 Let W : T'^ C be a function of class C'(KVr)- We say 

that W is the elementary operator associated with W , if W is a zero-order 
pseudodifferential operator with the symbol WiY^ £,2) = W{Y)-\-ii^2l'<^ W(Y). 

Proposition 5.4 Let a zero-order pseudodifferential operators 21 with symbol 
A{Y,S,), and an elementary operator W, satisfy conditions H3.12|) . H5.2|l . and 
& be a pseudodifferential operator with the symbol S'(y,^) = W{Y,£^2)A{Y,£^). 
Then there exist pseudodifferential operators OTa, '^a, iis and '^a, O.a, 27s so 
that 

1 

SlHi = ^ 2ljS^^' + OTaH + *nA, (5.14) 

3=0 
2 

2tS)i = ^21jD}"^ -H^Pa^ + Oa, (5.15) 
2 

6Di =^6jD}"^+i;s£-f2Js. (5.16) 

i=o 

Here 21^-, &j are elementary pseudodifferential operators associated with the 
complex-valued functions 

io(i')-A(y,o,i), ii(y) = i[%A](r,o,i), (5.17) 

MY) = -%)A + z/m^](y,0,l), (5.18) 

Sj = AjW for j = Q,l, §2 = A2W - Lm AoImW. (5.19) 
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For any u e H^iW^/T) with 1 < s < I ~ 4, 

WMauWs + W^AuWs + Wau\\s + \\Qau\\s < (5.20) 

c(s)(|a|3,.||u||o + |2t|3.3||^i||s-l), 

\\iisu\\s-i + \\'Ssu\\s< (5.21) 

c(s)((|a|3,3|2n|3,.+3 + |a|3,.+3||2n|3,3)||u||o + |a|3,3|2IT|3,3|k||.-l) 

Proof: The proof is given in Appendix ^ Let us turn to the proof of 
Theorem O 
First step 

We begin with the calculation of a commutator of an elementary operator 
and the left side of (|5.12|) . 

Lemma 5.5 Let W be the elementary operator associated with a function W € 
C'(RVr). Then 

(£ + 212)1 + Qs)!m = 2II£+ + e^Si +T + 5n, (5.22) 

where symbols of an elementary operator Wi and zero-order pseudodifferential 
operators &, 1. are given by 

Wi^2vdy^W, S = AW, 
T = udl^ W - i^idy, W - iT^2dy, W + Ady, W + BW + (5.23) 
(i{e-diA){(2dy,~r(idy,)W, 

where = (^2, ^Ci)- The remainder has the estimate 
\\mu\\s < c(s)(|l - 21J|o,6|k||.-l + |1 - 21J|o,s+6|k||o) + 

c(s)(|2n|o,6|2l|2,3||«||.-l + (|2IJ|o,.+6|2l|2,3 + |2IJ|o,6 |2t|2,«) |1 «|1 o) + (5.24) 
c(s)(|21J|o.4|»|2,3||«||.-l + (|2n|o,.+4|»|2,3 + |21J|o,5 1»|2,.) || u|| o) - 

Proof. It is easy to see that 

i^DlW = :^W^l+Wi^i+W2, (5.25) 

where 2112 is the elementary operator associated with the functions vdyi W. Next 
since (— A)~^/^ and 2IJ are first and zero order pseudodifferential operators, 
formulae for commutators (|F.6|) . ljF.8|) from Proposition IF . 31 imply 

(_A)i/22U = 2U(-A)i/2 + [(-A)i/2^2U]i + (5.26) 

where the symbol of the operator [(— A)^/^, 21J]i is equal to 

-i{^idy,+TC2dy,)W. (5.27) 
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< 111- 



and the symbol of the operator (—A) 



Moreover, since |1 — 2U|r 
does not depend on Y, inequahty IF.lOp yields the estimate 

W^^'^'^^^ulU < c{s){\\l - W\\ce+s\\u\\o + 111 - M^llcahll.-i). 



1/2 



(5.28) 



Next applying formulae ljF.5|) , l|F.7p to the composition of the pseudodifferential 
operators SlSi and 211 we arrive at the equality 



(5.29) 



where (212)1221)1 is a zero order pseudodifferential operator which symbol is 
given by formula ljF.7|) with A replaced by ikiA and B replaced by W. Noting 
that for A: ^ 0, 

kM^^i^e, hOk^^^-re^e with e^ = (6,-ei) 
we can rewrite expression (jF.7|) for the symbol (2t!l)i2n)i in the form 



AdyiW + £,i 



(5.30) 



Recall that it vanishes for fc = 0. Since 2lSi is a first order operator with 
|2lJ)i|} ; < c|2l|i,; inequality l|F.9|l from Proposition IF .31 implies the estimate 

IID^^il^ull, < c(|2t|2,.|2n|o,6 + |2l|2,3|21J|o.6+.) Il^llo + C|21|2,3|M^|0,6||"||.-1. 

(5.31) 

Setting 

T - 22J2 + [(-A)i/2, 21J]i + (2lJ)i2n)i + (^8211)0 
we obtain the needed representation, with the remainder 

and the lemma follows. ■ 
Second step 

Now we give a formal construction of the operators £, € and g'. We take the 
operator £ in the form 



G: = ^2ir(p)D7^'-ni, 

p=0 



(5.32) 



where elementary operators W^^^ will be specified below. Applying Lemma 1^31 
and noting that £ commutes with J)j we obtain 

+ +») 211(f) Dj^P = W^p^D^P£, + w[^^DI-p + 
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Here &'^p\ m^P^ are given by Lemma ICT with 2U replaced by W^p\ Com- 
bining these identities we arrive at 

(£ + 2l2)i + <B) (Hi +(£.) = (Hi + £)£+ 

2 

p=0 

(5.33) 

Recah that S'-p^ and T^^) are zero-order pseudodifferential operators, hence by 
Proposition 15.41 they have the decomposition 

2 

3=0 
1 

3=0 

in which the symbols of elementary operators Gj are given by the formulae 

H5.19|l with W replaced by W^p\ and the symbols of elementary operators 1^-^^ 
are given by formulae (|5.17|) . H5.18|l with A replaced by T. Substituting these 



relations into (|5.33|) we obtain the identity 

(£-^2lS)i +Q3)(ni +£) = (Hi + e:)(£- x)+ 

where the elementary operators 3j and the reminders IH" are given by 



(5.35) 



















:J2 








33 


= 6« 4 


- 6f ) ^ 





p— p—O p—O 

Noting that 

S)r' = (-A)-i/2d-2£ _ i.(-A)-i/2. 



(5.36) 
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we can rewrite (|5.34ll in the form 

2 

(£ + + S) (Hi + €) = (Hi + e){£- >c)+Y^ ^pD}"^ + d (5.37) 

with 

e = € + m' + i3(-A)-i/2j)-2^ 

Now our task is to show that 3j, j <2, vanish for an appropriate choice of 
elementary operators 211^^' . Note that the operator equations Jj = are equiv- 
alent to the scalar equations Ij{Y) = 0, in which the complex- valued functions 
Ij are associated with the operators Jj. This observation along with formulae 
H5.35|l leads to the equations 

W^"^ + S^"^ = 0, 
W^'^ + S^o'^ + S["^ + fo<°) + - 0, (5.39) 

w^f ) + + s!f'> + ^« + f(") + f + ^ 

By Proposition^ we have fg^^^ = T'-p\Y,0, 1) and f^''^ = i^-i[%r(p)](y,0, 1), 
which along with equality (|5.23(l yields 

f^p^ = vd^w'-p^ - iTdy.w'-p^ + Aody.w'^p^ + BoW'^p\ 

f'f^ = -i^dy^W^P^ + 2Axdy^W'^P^ + Bidy,W'^P\ 

Substituting these identities into H5.39|l and using we obtain the 

recurrent system of ordinary differential equations for functions W'''p\ p = 0, 1, 2, 

(21/9^, +io)W^(°) =0, (5.40) 
{2ydy,+Ao)W^^^+g,=Q, j = l,2, (5.41) 

where 

gi = [vdl^W^^^^ + A^dy.W^^^ + BoW^"^^ - iTdy^W^^^) + ^W^^\ (5.42) 
52 = (iaVK^") +iiVK(i)) + (2ii5y,VF(°) +Bi9j,,VK(") -li^y.iyf")) 
+ +^o^!/i^^^^ +^oW^^^^ -iray.iy^i)) +>{W'^^\ (5.43) 

By the equivariant property, the function is odd in yi and IIi^o = Aq. 
Therefore, the general solutions of homogeneous equation (|5.4U|I has the form 

W^°\Y) =C{y2)a-{Y), where a±(r) = exp(±^£)^iio) (5.44) 
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and C(y2) is an arbitrary function. On the other hand, inhomogeneous equation 
(|5.41|l has a periodic solution if and only if 



gj{yi-,V2)a^{yi,V2)dyi^0 for aU 2/2- (5.45) 

Now we aim to show that solvability condition (|5.45l) is fulfilled for an appro- 
priate choice of C. Substituting (|5.44|) into the expression for gi and next to 
(|5.45|l we obtain the ordinary differential equation for C, 

-i27rC" + 27rC>f - iCr j dy^a^ dyx + Cu J a+dy^a^dyi + 

— TT — TT 

+C J Ana+dy.a'' dyi + C J Bq dyi = 0. 

— TT — TT 

Noting that 

TT TT 

J a+dy^a^dyi J dy^Xl^^Ao dyi = 0, 



1/ 



a+dy^a dyi + j A^a^dy^a dyi = J Aldyi, 

TT — TT — TT 



we can rewrite it in the form 

C'(y2) + %2)C(y2)=0, (5.46) 

with the coefficient 

TT 

b (y2) J {Bo- ^Al) dyi + IK. (5.47) 

It follows from formula (|5.7I) for >f, that the mean value of h over a period is 
zero. Therefore, equation 15.46|l has a periodic solution 

C = exp(-2)2"i6). (5.48) 

In this case the particular periodic solution to (|5.45|) is 



=_i-a-D-i(gia+). (5.49) 

Next note that, by the equivariant property, the functions Aq, A2, and Bi are 
odd and the functions ^1, Bq are even in yi. Hence a±, ll'(o) are even, and 
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W^^^ is odd in yi. In particular, g2 is odd in yi and automatically satisfies 
solvability condition H5.45|l . From this we conclude that the operators 3j, j <2 
vanish when 

VK'"' = exp(-D^i&)a-, W^^'> ^ ~^a-D^^ (^gja+y j = l,2, (5.50) 

which along with H5.37|l leads to identity (|5.12|l . 
Third step 

It remains to show that operators introduced above are well-defined and 
satisfy inequalities H5.13|l . We begin with the estimating of the functions W'^p\ 
Recall that for all smooth functions a,b G C"'(R^/r), 

||1 - expallc- < c(||a||co)||a||c-^, ||a6||c= < c(s)||a||co||6||c- + lla|icHl^llc"- 
It follows from this and and lfCTS|l that 

111 - W^°'>\\c^ < c(||io||co, ||Bo||co)(||io||c^ + ||Bo||c=). (5.51) 
On the other hand, (|5.42|l implies the estimate 

||<?i||c= <c(||l-iy(")|| Cs+2 + ||^o||co||W^^°''||c= + i + 
\\Ao\\c4W^'^\\c^ + WMc^) + {\\Bo\\c4W^'^\\co + ||So||co||W^(°^||c=). 
Combining it with H5.51|l and noting that 

||M^(o)||ci < c||M^(")||co +c||T^(°)||c= 

we arrive at 

lll^'llc^ <c(||io||c",||So||co)(l|io|| 

Cs+2 + \\Bo\\cs+2 J. (5.52) 

Arguing as before and using H5.43|) we obtain 

||52||c= < c(||io||co, ll^ollco, WMco, WMC, ll^i||co)(||io||c-+^ + 

Po||c=+4 + ||il||c»+3 + \\Bi\\cs + 2 + ||i2||c=), 

which along with H5.50|l leads to the estimate 

2 1 

lll^^'^llc-^ <c(||i,||co,||S,||co)(^||i,||c=+4 + 5]||S,||c.+4). (5.53) 
Noting that ||^j||c- < Ma.s, we conclude from (|F3T|) . lf^3^ . that 

111 - iy(°)||c= + IIW^^'^llc^ + WW^^^WC^ < C(|a|3,0, 1*813,0) (|2l|3,.+4 + l^|3,s+4 , 
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which leads to 

|1 - !m(")|„,,,, + |2IJ«U,, + |2»(2)|„,, < c(|2l|3,0, |»|3,0)(|2l|3.s+4 + m^.s+i) 

(5.54) 

for all ni > 0. Now we can estimate the norm of the operator £. Since = £, 
Proposition IF . II along with H5.54|) implies 

\\€u\\, < 11(1 - 21j("))u||, + WW'^'^uWs + ||2IJ(2)u||, < 

C(|a|3,0, |'8|3,0)((|2l|3,s+4 + |»|3.s+4)h||o + (|2l|3,4 + |3,4) II «ll ■ 

(5.55) 

Let us estimate the operators €, ^. First note that inequalities (|5.20() H5.21|l 
from Proposition 15.41 imply the estimate 

\\iXsi.>u\\s-l + W^s^^Mls < C(|a|3,0, |*B|3,o)((|a|3,3|21j(f)|3.s+3 + 
+ |2l|3,.+3|2B(f)|3,3)||"||0+|2l|3,3|21j(^')|3,3||^i||s-l). 

From this and ()5.54|) we conclude that 

||il5(p)u||,_i + ||5Js(p)u||, < (5.56) 

C(|2l|3,7, |»|3.7)((|2t|3,.+7 + l'B|3,.+7) hllo + {^3,7 + l^bj) II^IU-l- 

Next ()5.54|l along with inequality (|5.24|l from Lemma l5 . 51 implies the inequality 

Wm'^P^ulU < c(|a|3,io,|'B|3ao)((|2t|3,.+io + |»|3,.+ io)||w||o + 

+ (|a|3,io + |'8|3,io)||«||s-i). (5.57) 

Recalling that for arbitrary zero-order operators 2l,22J and the operator 6 with 
a symbol S = AW, 

\6\m,s < c{s,m){\mU,o\W\rn^s + mm,s\Mm,o), 

and using formula (|5.23|) for symbol 1 we obtain 

|5:^^^|3,« < |l-iOT('')|3,.+2 + (|2l|4,0 + |'8|3,o)|2Ij(P)|o.s+l + (|2l|4,. + |»|3,.)|2Ij(P)|o,l, 

which being combined with (|5.54() gives 

|2:(^^|3,. < c(|2l|4,o, l»l3,o)(|2t|4,.+6 + l^ks+e). (5.58) 
In particular, inequality H5.58|l along with 15.20|l yields the estimate 

||9JIt(p)u||. + Pt(p)w||. < (5.59) 

C(|2l|4,0, |»|4,o)((|2l|4,.+6 + |»|4,.+6|)||w||o + (|2t|4,9 + |*B|4,9) ||"|U-l) • 
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Similar arguments applying to the formula H5.23|l for the symbol of the operator 
6 give the estimate 

|6(^')|3,. < c(|a|3,ol»l3.o)(|2l|3,s+4 + |'B|3.s+4). (5.60) 

Finally estimate J^i-^y^^^- Applying f^I^ . fE^HOI) to we arrive 

at 

Psis.s < C(|2t|4,0, |'B|3,o)(|2l|4.s+6 + |»|4,.+6), 

which along with inequality ljF.l|l from Proposition lF.l| implies 



p3i-A)-'/^u\\s < C(|2l|4,0, |5B|3,0)((|2t|4,.+9 + |'B|4,.+9)h||0 + (|2l|4,9 + |»|4,9)|k||.-l) 

(5.61) 

Next note that 

\\€u\\s<\\Ms + \\M^^)-'^Ms + \\'n'u\\s, 

\\Ms<p3{-^r'^Ms + m"u\\s. 

It follows from (|CT7|) . and |[?3^ that 

2 1 
\\m'u\\s < J2 \\^Sip)u\\s + J2 W'^TiP^uWs < 

C(|2t|4,9, |*B|4,9)((|2t|4,9+. + |*B|4.9+.)||"||0 + (|2l|4,9 + 1*814,9)1^11. 
2 1 

ii5H"uiu < w'^sipMis + ii^T(p)"iu + m^^^uWs < 

p=0 p=0 
C(|2t|4,10, |*B|4.10)((|a|440+. + |*B|4,10+s)||w||o + (|2t|4,10 + l*B|4,io)||i 

Combining these results with H5.61I) and (|5.13|) we finally obtain 

WMs + W^uWs < C(|a|4,9, |'B|4.9)((|2l|4.9 + |»k9)|klU + 
(|2t|4,.+9 + |'Bks+9)h||o), 

ll^wll. < c(|2l|4ao, |5B|4,io)((|a|44o + |»|4ao)||w|U-i + 

(|2l|4,.+10 + |*B|4,.+lo)||w||o), 

which gives ()5.13(l and the theorem l5.2l follows . 
5.2 Proof of Theorem lOl 

The proof is based on the following lemma on the invertibility of the operator 
1 + € from Theorem 
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Lemma 5.6 Under the assumptions of Theorem \5.1i there is a positive £i 
depending on r and I only, so that for all e G (0,ei), the operator 1 + £ 
has the bounded inverse (1 + (E)-^ : H^^iR'^/T) ^ H^^{R'^/r) and, for all 

11(1 + < c||ii||^, (5.62) 

11(1 + < c\\u\\r{Cs + |2l|4.s+10 + |5B|4..+ io) + cMs- (5.63) 

Proof. Formally we have 

oo 

(l + £)-l = ^(-l)"e. 

n=0 

By Theorem E3 the operator £ : H';;^^{M.'^ /T) iJ„^,.(RVr) is bounded and 
its norm does not exceed ce. It follows from this and inequality H5.13|l that for 
ce < 1, 

< c£|l€"-iw|U + c(|2t|4,.+io + |'B|4,.+io)IHS""'"ll. < 
ce\\<B^-\\\s + (c£)"-i(|2l|4,.+io + |'B|4.s+io)||^i||. < 
{ce)^\\(B"-Ms + 2(ce)"-i(|2t|4,.+io + |Q5|4,.+io)hl|r < 
... (ce)"-i(||u||, + n(|2t|4,.+io + |'B|4,.+io)h|| 

which leads to the estimate 



< c(\\u\\s + (|2t|4,.+io + |*B|4,.+io)hl|r.) (l + ^ n(ce)"-i 



n=0 

It remains to note that for e < 1/c, the series in the right hand side converges 
absolutely and the lemma follows. ■ 

Let us turn to the proof of Theorem l5.1l Since 21 and 5B satisfy all hypotheses 
of Theorem l5.2l the corresponding operators £, (£, ^ are well defined and meet all 
requirements of this theorem. Moreover, condition (iv) along with inequalities 
H5.13|l yields the estimates 



+ ll^ulU < ce||u|l,, + c(|2l|4,.+io + |»|4,.+io)llu||o, 
\\Su\\s < ce||u||s-i + c(|a|4,s+io + |S|4,s+io)||w||o, 



(5.64) 



\\€u\\t + \\iBu\\t < ce\\u\\t, for t e [2,r]. (5.65) 
We look for a solution to basic equation (|5.1|l in the form 

u = X^p''^^ + {1 + €)v with veH^-^^-^. (5.66) 
Substituting this representation into H5.1|) we obtain the equation 

A(£ + i3)^(°) + (£ + io)(l + (!:)i; = /. (5.67) 
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Since Hi coincides with the identity mapping on H^^{M.'^/T), it follows from 
Theorem E21 that 

(£ + i3)(i + e:) = (! + €)(£ -^) + 3 on i/;,(RVr), (5.68) 

where 

3 = 5 + £-i(i + e:). 

Hence we can rewrite II5.67|) in the equivalent form 

A(£ + + (1 + (£)(£ - >t)v + 3w = /. 

Applying to both sides the operators Q and 1 — Q and noting that Q£'i/;'^°) = 
we obtain the system of operator equations for a scalar A and unknown function 

Q(l + €)(£- + Q3« = Q(/ - A^V^"^), (5.69) 
A(l -Q){L + + (1 - Q)(£ + + €)v = (1 - Q)/. (5.70) 

Our aim is to resolve the first equation with respect to v. Hence the task now 
is to prove the solvability of the equation 

Q{l + €)i£-M:)v + Q3v=Qg. (5.71) 

The corresponding result is given by 

Lemma 5.7 Under the above assumptions, there is a positive S2 depending on 
r,l only, so that for g G ^{M.^ /T) and < e < 62, equation H5.71|l has the 
unique solution satisfying the inequalities 

lkl|r-i < c||.9|l„ (5.72) 
\Hs-i < c\\g\\r{Cs + |ak.+io + l^ks+io) + c||.g||.. (5.73) 

Proof. We look for a solution to equation (|5.71|) in the form 

V = Q{£.- >c)-\l + <a)-^Q<f (5.74) 

with the new unknown function ip £ . Since the operator £ commutes 

with the projector Q, we have 

2(1 + (£)(£- >^)Q(£ - x)-\l + (£)-iQ = QXoS + Q (5.75) 

with 

Xo - -{<^Qf + (1 + €)q((1 + £)"^ - 1 + . 
Hence (|5.71|l is equivalent to the equation 

if + QXQif = Qg, (5.76) 
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in which the operator X is defined by 

X - Xo + 3Q(£ - + 

We use the regularization method to prove the existence and uniqueness of 
solutions to H5.76|l . and consider a family of regularising equations depending 
on small positive parameter S, 

f + QXQ(-A)- V = Qg-. S>0. (5.77) 

Let us estimate the norm of the operator X. It is easy to see that Lemma |5. 61 
implies the estimates 

llXo(/^|l. < ce^Mr, (5.78) 
W^ovWs < c|l^||,(|a|4,.+io + I^Im+io) + ce^ll,^!!,. 

On the other hand, inequality (|5.8|l along with Lemma leads to the estimate 

||Q(£->^)-i(l + (£)-iQ(^||,_i <c||^||,, (5.79) 
+ £)-iQv'll.-i < c||(^||,(|a|4,.+io + I^Im+io) + cMs- 

Next estimates (|5.64|l for € and g' along with inequalities H5.6|l for £_i imply 

||3(^||r < ce||(^||,_i, (5.80) 
113(^11, < cy\\r-i{Cs + M4.S+1Q + I^U.^+io) + ce\\(p\\s-i. 

Combining (|5.78|) - (|5.8U|) we finally arrive at 

\\Xip\\r < Ce\\ip\\r, 

IIX^IU < cMr{Cs + |2l|4,s+io + I^U.^+io) + ceMs. 

Since the operators (— A)^*^ are uniformly bounded in ^{M? /T), it follows 
that each solution to equation (|5.77|l satisfies the inequalities 

\\(p\\r{i - ce) < c\\g\\r, 
119511.(1 - ce) < cMr{Cs + |2t|4,.+10 + |»|4,.+io) + 

c\\g\\s + c||.g||r(|2l|4,s+10 + |*B|4,s+lo)- 

Hence for e < l/2c, 

\Mr<c\\g\\r, (5.81) 
M\s < c\\g\\s + c\\g\\r{Cs + Mi^s+w + |5B|4,.+io), 

In particular, the solution is unique. Since the operator 



63 



is compact, uniqueness along with the Fredholm Theorem imphes the solvability 
of equations H5.77|l for all 5 > . Hence they have a family of solutions (^5 
satisfying (|5.81|) . After passing to a subsequence we can assume that 

LPS LP, {-Ay'^ifis ^ (fi weakly in H^R^/T) as 6\0. 

Obviously cp serves as a solution to equation H5.76(l . Recalling H5.74(l we obtain 
the solvability and uniqueness result for H5.71|l . It remains to note that estimates 
(jSZZSIl, follow from formula (|57^ and inequalities JEZHl) , EHU ■ ■ 

We are now in a position to complete the proof of Theorem 15.11 Applying 
Lemma I5. 71 to equation H5.69|l we obtain the representation 

v^Xvo+Vf, (5.82) 

in which vq and Vf are solutions to equations H5.71|l with g ~ —QSjip^^^ and 
g^Qf. It follows from fH7!^ that 

ll^^/ll.-l <c||/||.r, 
\\Vf\\s-l < c\\f\\r{Cs + |a|4,. + 10 + |»|4,s + 10) + c||/|U (5.83) 

||t^o||r-i < ce, 
\\vo\\s~i < c{Cs + |2t|4,.+io + I^U.^+io). 

The functions vq and Vf are completely defined by the eigenfunction -ipQ and the 
right hand side / of equation H5.1(l . but the scalar A still remains unknown. In 
order to find it we substitute representation (|5.82|l into (|5.70(l to obtain 

A(l - Q)(£ + i3) (V^(°^ + (1 + €)vo) = (1 - Q) (/ - (£ + jo)(l + €)vf) . (5.84) 

It follows from identity H5.68|l that equation (|5.71|) for vq is equivalent to 

Q(£ - }i)QvQ + Q((£ + io)£ + f) + k)Qvo = -Q^i^^°\ 

Since the operator Q{£ — k) has the bounded inverse Q(£ — >f)^^ : H^'^ i-^ 
H^~^'-^, we have 

vo + Qi&- ^)-^Q{{£ + Sj)€ + Sj + K) Qva - - Q(£ - }<)-^ QSjyj^"'^ . (5.85) 

Next, using inequalities H5.83|l along with (|5.65(l and (|5.72(l and noting that for 
2<t <r, 

\\Sju\\t-i < ce\\u\\t 

we arrive at 

\\Qi£-^)-^QS)^^°'^\\r-2 < ce, \\Qi£~^)-^{{£ + Sj)€ + Sj + x)vo\\r-3 < cs\ 
which leads to 

vo = -Q{2- }<)-^QSj^^'''^ +v' with ||t/||._3 < ce'. 
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Substituting this expression into H5.84|l gives 

A(l-Q)(£+i5)(V'(°)-Q(£->^)-iQi3V*°)) W = (l-Q)(/-(£+i3)(l+£)«/) 

(5.86) 

where 

v" = (1- Q)(£ + io)(w' + e:wo). 

Next noting that 

\\{l-Q)^u\\s<ce\\u\\o 
and using estimate (|5.65|l we obtain 

||w"||r-4 < ce||t;'||^_3 + ce^\\vo\\r-2 < ce^- 

Multiplying both sides of (|5.84|) by V'*''^'' and integrating the result over leads 
to equality 

>^K = j (/- (£ + i3)(l + c:)^^/)V'^°^dr, (5.87) 

T2 

in which 

j (£+i3)(V'^"^-S(£->*)"^S^^^°^)V'^*'^ dY+ j w'V^"^ dY = @e'^+0{e'^). 

Hence there is a positive ei < £2 depending on r, I and i/'o only so that 2\K\ > 
for all < e < Si. From this, relation H5.87|l and inequalities (|5.83() we 
conclude that for such e the unknown A is well defined and 

eM < ^11/11.- (5.88) 

Hence for < e < £1, equation l|5.1|l has the unique solution u E Hoe- It 
remains to note that inequalities H5.83|l and (|5.88(l along with the identity 

u = X^'-°^ + A(l + €)vo + (1 + €)vf 

imply estimates H5.1U|I . (|5.11|) and this ends the proof of Theorem 15. II 

5.3 Verification of assumptions of Theorem 15.11 

In this subsection we check all abstract conditions required for solving the linear 
equation (|5.1I) . in the case when the operators involved in this equation are 
defined by Theorem 13.41 Note that the symmetry conditions (|5.3|) and (m) 
follows from assertion {Hi) of this theorem. The same conclusion can be drawn 
for the Metric conditions (|5.5|) . H5.6|l which are realized, as soon as ||J7||p < e 
and p = I + 3. 

Let us consider the restrictions on the spectrum and resolvent of £. Part (v) 
results immediately from Theorem 14. 191 restricted to the space of functions odd 
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in t/i and even in t/2, and once we notice that the eigenvalue corresponding to 
the eigenvector ip^^^ = sinj/i cos?/2 satisfies 



+ + t2 = e^iyiir) + 0{e^). (5.89) 

This results immediately from the fact that the linear operator C (see ()3.3|l 
corresponds to the differentiation of the basic system at the point 

U = t/J^) + e^W, N>3 (5.90) 

where J/i^' is the approximate solution at order e^, and from Lemma 13.71 
Moreover the coefficient vi is positive (larger than a positive number) for any 
value of r. 

Now, from Lemma |3.8I we see that the Fourier expansion of the diffeomor- 
phism of the torus has, at order e, only harmonic 1 terms in Y, and at order 
only harmonics 2 and in y = (j/i, 1/2), all terms being invariant under the shift 
: Y i—f Y + {tTjIt). Let us now consider the expression of the linear operator 
£, + 9) which is obtained after applying the above diffeomorphism on a linear 
operator which may be formally expanded in powers of e, having coefficients 
with the same property as the diffeomorphism in terms of harmonics in Y . The 
result is that the formal expansion in powers of e of has the same property 
as above i.e. the order is independent of Y, while orders e and only contain 
respectively harmonic 1 and harmonics 2 and in Y. 

Let us consider the formula 1)5. 7|l giving the coefficient >c. In the integral over 
the functions A and B are 0{e) and can be expanded in powers of e, their 
principal part containing only harmonic 1-terms in Y. It results immediately 
that 

K^O{e^). (5.91) 

It is not useful in our proof to give more precision on this coefficient, however 
the interested reader might check (after few days of computations) that 

^ = £^f^{8T^^5T^--} + Oie^). (5.92) 
16 4 

Now, from the estimate on operators 21 and 05 in theorem 13.41 and from 1)5. 7|l 
we have 

1^1 <C||C/||14, 

hence from ()5.91)l and 1)5.90)1 we deduce that 

K = £^H, < c(l + e||M^||i4). (5.93) 

It should be noticed that in formulae ()5.89)l and 1)5.92)1 the terms of order O(e^) 
depend on the unknown W of ()5.9U)) which is assumed to be bounded in Hj^^ 
(because of the smoothness required for computing operators 05 and £_i above). 
Hence, if we want to apply the result of Theorem l4.2l for obtaining the required 
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estimate (|5.8|l for {£, — >c)~^ we need to show that properties H4.4I) hold for the 
functions of e 

Vj{e) — — e^vi + e^Vj{e), x = e'^>Cj[e) 

where vq = (1 +t^)^/^, vi — vi{t), I'jie) and ><j(£) are obtained with some Wj 
in H^g^ through the iteration of the Newton method of Nash Moser theorem. 
Because of the smooth dependence of these coefficients in function of (e, U) (see 
in particular Theorem 13 . 41 for v, and (|5.93() for x), the properties (|4.4|1 for i'j{e) 
and >«j(e) are verified as soon as there exists R> such that 

\\W,is')-W,ie")\U, < R\e'-s"l (5.94) 
||l¥,+i(e)-l¥,(£)||i4 < 

holds. This property will be checked in sectional Assuming that this is true, 
we have all conditions of Theorem 14 . 21 realized . This completes the verification 
of restrictions on the spectrum and resolvent of £. 

The nondegeneracy condition H5.9|) is fundamental for obtaining the esti- 
mates (|5.10|) . (|5.11|l for the solution of the linearized system. So, let us now 
consider the coefficient @ defined in (|5.9|) . We notice that the eigenfunction 
= sin yi cos 1/2 only contains harmonic 1-terms in Y, and assuming that W 
in formula H5.90|l is bounded in Hjlj , and denoting the coefficient of order e in 

the operator by S)^-^"^ we have 

and, since £0 is invertible on finite Fourier series orthogonal to one obtains 
i3Q(£->^)-iQi3V'*°) =£2^(i)£o'^''V'"^ +0(e3). 

Now defining the coefficient of order in the operator £ + by i^'^-' we have 
@ = J (io(2)^(o) _ i3(i)£-ii5(i)^(o)) ^(0) dY. (5.95) 

Then, we prove the following 

Lemma 5.8 The coefficient @ of nondegeneracy condition is given by 

@ = -27r2^ = ^(ao+/3o), 

where /ii(r) is given by 1^2.11]) and (oiq + (3o){t) is given at Theorem \2.t\ This 
coefficient is non zero for t ^ t^- 

Proof: the proof is made at Appendix ^ 

This completes the verification of Nondegeneracy condition. 
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5.4 Inversion of £ 

Let us now consider the inversion of the operator H defined in H3.3|l 

C{U,^x)V^F, (5.96) 

where 

and where we look for 

SU — {Si/;, Srj), Si/j = 64) + bSij, 

where b is defined in In this subsection we prove the foUowing Theorem 

which coUects the results of previous sections on the operator £ 

Theorem 5.9 Consider M > , t > 14, and tt/A > S > 0, and set 

U = U^/^ + e^W, > 3, (5.97) 
^^ = A'f ' =A*c + eViW + 0(e3), 

where \\W\l < M, and (J/i^^ , /li^-*) is the approximate solution at order 
obtained at Theorem in the case of diamond waves (ei = £2 = and 
T = tanfiig, /ic = cos 6*0, 6 < Oq < tt/2 — S. Assume moreover in j5.^^7| ) that W = 
Wj{e), j € N satisfies the property ^5.94}) , then for fi^^ £ a g (0, 1/78), 
there exists eo > and a subset £ of [0,eo] such that for any F G H^^-j, s > 5, 
and for e £ E, the linear equation i|5. 96\) has a unique solution V corresponding 
to 6U e IH^g-)'^, such that the following estimates hold 

m\\r-3 < ^\\F\\r, 5<r<t-8, (5.98) 
\\SU\U-3 < ^\\FUl + e''\\W,\U^s)+cis)\\F\U. 
Moreover the following property holds for the "good" set £ : 

-^j^ meas{/^ = e £n(0,r)} -> 1, as r -> 0. (5.99) 

Proof. By construction, we have 

gn{S(b)-J*i-J{S(j)))=h, 
a 

where 

h = f-J*ig/a), 
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and which, after the change of coordinates of Theorem 13. 41 becomes 

(£ + 2l2)i +'8 + £_i)<50= -, 

which we learnt to invert at Theorem 15. II 

First, for U G Et^^-j, and \\U\\r < Mr, we have the following estimates, due 
to Lemma FH. II 

\\h\U<c{M,){\\F\U, + \\FMU\U2), 
and with Theorem 13 . 41 we obtain 

|N|.<c(M4)(||/i||. + ||C/||.+4||/l||o). 

Taking into account of the estimate on k{Y) in Theorem l3.4l we then arrive to 
11-11. <c(M5){||i^||,+i + ||F||2||{/||.+5}. (5.100) 

K 

In the same way, thanks to the Theorem 13.41 we also have 

\mu<ciM,){\MU + \\u\u,\mo}, 

and since 

a a 

we obtain 

\M\s < c(Af4){||^||s+l + ||[/|U+5||^||o} + 

+c{Mi){\\F\\s + \\U\U+2\\F\\2}. 

Hence 

\\V\\s < c(A/4){||^|U+l + ||C/||s+5||^||o} + 

+c(M4){||F|U + ||C/||.+2||F||2}, 

and thanks to Lemma l3. II 

\\5U\\s < c(M3){||y|U + ||C/||.+i||y||2}, 

< c{Mr){\M\s+i + ||C/||s+5||<50||o + ||C/||.+i||50||3} + 

+c(M4){||F||, + ||C/||,+2||F||2}. (5.101) 

It remains to use Theorem 15. II which connects S<j) and ^, taking into account 
of estimates ii) of Theorem 13.41 Assuming that properties (|4.4|) hold for za, (e) 
and %(e) we have (see Theorem 14. 2|) for e < Eq and e G a G (0, 1/78) the 
following estimates for ||C/||p < e, i > 14, 1 < r < s 

ll^ll.-i < "-^Whr, l<r<^~9 

IMU-l < ^^^11-11.(1 +||C/|U+19)+c(Mi4)||-|U. 
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From this and (|5.1UU|) we deduce that 

11^11.-2 < ^^||F||.(l + ||[/||,+i8)+c(Afi4)||F|U 

and from (|5.1U1|I we obtain the estimates for 5U (|5.98|l . The rest of Theorem 
15.91 foUows directly from the results of previous subsection and from Theorem 

6 Nonlinear problem. Proof of Theorem 11.31 

In this section we complete the proof of the main Theorem 11.31 on existence of 
diamond nonlinear waves of finite amplitude. To this end we exploit the general 
version of the Nash-Moser Implicit Function Theorem proved in Appendix N of 
1^4, . This resuh concerns the solvability of the operator equation 

$(W^,£) = (6.1) 

in scales of Banach spaces Es and Fs parametrized by s G Nq = N U {0}, and 
supplemented with the norms || • ||s and |-|s- It is supposed that they satisfy the 
following conditions. 

(Al) For t < s there exists c(<, s) such that 

\\-\\t<c{t,s)\\-U \-\t<c{t,s)\-\s. 

(A2) For A e [0, 1] with At + (1 - A)s e N, 

IMUt+(l-A). < C{t, s)\\-\\^ |-Ut+(i_A). < C{t,s)\-\t^ \-\\~^. 

(A3) There exists a family of smoothing operators Sp defined over the first scale 
such that for p > and t < s, 

\\SpW\\t < c(t, s)\\W\\,, II SpWWs < c{t, s)p'-'\\W\\u 
\\SpW-W\\t<c{t,s)p'-'\\W\\,, 

and, if e I— > p(e) is a smooth, increasing, convex function on [0,cxd) with 
p(0) = 0, then, for < ei < £2, 

II {S^^,,^ - Sp^,,-^)W\\s < c{t,s)\e^ ~ e2\p\e2)p{eiY-'-'\\W\\t. 

(Bl) Operators <I>(-,e), depend on a small parameter e g [0, Eo]; and map a 
neighborhood of in Er into Fp. Suppose that there exist 

CT < p < r - 1, cr, p, r e No, 

and, for all / E Nq, numbers c{l) > and e{l) E (0,eo] with the following 
properties for all W, U, Wi, Ui E B and e, E [0,£o]j * = 1, 2, where 
B^{W EEr-. \W\r < Ro} for some Rq > 0: 
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(B2) The operator $ : B x [0, £o] — * Fp is twice continuously differentiable, 

\^{W,E)\p+l<c{l){l + \\W\\r+l) (6.2) 

and, for W,U € Er+u e e [0,e(Z)], 

\D{W,U,e)l+i < c(0(l + ||l^||.+i + ||C/|Ui)||V^-t/||^ + 

+C{1)\\W - U\\r\\W - U\\r+l, (6.3) 

where 

D{W,U,e) = <^{W,s) - $([/,£) - <^'nr{U,e){W - U). 

Moreover, 

\D{WuUuei)-D{W2,U2,S2)\p < C (|£i - £2! + llT^l - W^2||r + - W^2||0 

•(||W^l-C/l||r + ||W^2-C/2||r). (6.4) 

(B3) There exists a family of bounded linear operators A(W, e) : — > Fp, 
depending on {W,£) G B x [0,£o], with 

\A{W,s)U\p<c{0)\\U\\r, U&Er, (6.5) 

that approximates the Prechet derivative as follows. For W G Er+i fl 
S, £e [0,£(/)] and f/ e Er+u 

\h{W,e)U-^'^{W,e)U\p+i < c{l){l + \\W\\r+m{W,e)\r\\U\\r + (6.6) 

+c{l)\^{W,e)\r+l\\U\\r + c{l)\^{W,e)\r\\U\\r+l. 

(B4) When e B n E^+u £i e [0, £(0], i = 1, 2, 

|$(W^l,£l)-*(T^2,£2)U+i < c{l){l+\\Wi\\r+l + \\W2\\r+l) (6.7) 

•(|ei - £2| + IIW^l - W2II.) + C(0||M^1 - W2\\r+U 



\{^'w{Wuei) - ^'w{W2,S2))U\p+i + |(A(W^i,£i) - K{W2,S2))U\p+i < 

<c(0(||W^l-W^2||r+i + (|£l-£2| + ||W^l-W^2||r)(||VFl||r+/ + ||W^2||r+0)l|C^I|. 

+ (kl - S2\ + \\W^ - W2\\r)\\U\\r+U (6.8) 

• A set 5 C fO, 00) is dense at if lim / £ d£ = 1. 

r\,o J 

£n[0,r] 

(B5) If a set £ C [0, e{l)] is dense at and a mapping ■!? : f — > B fl Ej.+i is 
Lipschitz in the sense that for £i, £2 S f , 

||i?(£i) - i?(£2)||r < C\e-i - £2] where C = C(i?), constant. 
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then there is a set C £, which is also dense at 0, such that, for any 
e e £{'d) and / e Fp+i, the equation K{d{e),e)8W ~ f has a unique 
solution satisfying 

||<5M^||p-.+; < e-^cmfl+i + Me)\\r+i\f\p)- (6.9) 

(B6) Suppose that -do : Sq ^ BnEr+i and mappings i!)k ■ njilg^f BHEr+i 
satisfy, for a constant C independent of fc g N sufficiently large, 

||^fc(£i) - ^k{e2)\\r < C\ei - £21, £1, 62 G n^Zo£{i9j), 
- M6)\\r < ^, ee n^^o^id,). 

Then I^j^qE {-d j) is dense at 0, where the sets are defined in (B5). 

Theorem 6.1 Suppose (Al)-(B6) hold and, for iV £ N with N > 2, equation 

H6.1|l has approximate solution W = We^^ G n^gNo^s, with, for a constant 
k{N, s), 

\\W^^^\\s<kiN,s)e, Me,W^^%<HN,s)\ef+' (6.10) 

and 

IIW^W _ w^PlU < k{N,s)\e^-e2\. (6.11) 
Then there is a set £ , which is dense at 0, and a family 

{W = i?(e) -.eeSjdEr 

of solutions to H6.1(l with ||'!9(ei) — i9(£2)||r < c\ei — £2] for some constant c. 

In order to apply Theorem 16.11 to the 3D wave problem, let us introduce 
some notations. Fix an arbitrary a € (0, 1/78) and < (5 < 1. Denote by Af 
the set of all /ic so that — ji^^ belongs to the set given by Theorem 
14.21 Since OIq. C [l,oo) is a set of full measure, 01 is the set of full measure in 
(0, 1). Choose an arbitrary /Ltc = (1 + t^)~^/^ G 01 so that r e (5, 1/^) and set 
T = tan^. It is clear that /ic and r meet all requirements of Theorem II .31 

Next we fix the lattice T such that the dual lattice is spanned by the wave 
vectors (1,±t), and set 

E, = F,^ K^M^/v) X ff|^,(RVr). 

Since the scaling mapping u ^ u o establishes an isomorphism between 
i?*(R^/r) and a closed subspace of the Sobolev space if* of doubly 2ti— periodic 
functions, the properties {Al) and {A2) are clear. A smoothing operator with 
the required properties can be defined by 

E <P\K\)uKe'''-'' , 

ZTT ^ — ' 
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where <^ : R"*" i-^ M"*" is a smooth function which equals 1 on [0, 1] and on 
[2,oo). 

Fix an arbitrary iV > 3 and define the operator (f> — ($i, $2) by the equaU- 

ties 

=£-^.F([/f^) +£^W^,/x(2^)), (6.12) 

where {Ue^'^\ /ii^^"*) is the approximate solution at order obtained at Theo- 
rem l2.3l in the case of diamond waves (ei = £2 = e/2), and r = tanfl, /ic — cos9, 
6 < T < 1/S. By construction W^i^'' — 0, hence (|6.10ll and H6.11|l are satisfied. 

It follows from Lemma fl. II that the continuous mapping {W,e) 1-^ ^{W,e) 
from Es+i X R into Fg for s > 2, is of class of C°°. Applying the same arguments 
as in section 9 of we conclude from this that for 5 < p < r — 18, the operator 
$ satisfies Conditions (Bl) and (B2 ), and inequality H6.7|) from Condition (B4). 
Let us denote by A the approximate differential, defined as 

A{W,e) = C{U^^^^ +e^W,4^^^) 

where the operator £ and coefBcient b are defined by formulae (|3.1|) and H3.3|) . 
It follows from this and H3.2|l that the linear operator 

n^A{W,e)-dw'S?iW,e), 

is defined by USW = {RiSW,0) with 

where 5W — {Sip, Srj) and $1 ~ ^i{W, e). Then, thanks to Lemma ITTl we have 
for s > 2, and \ \U\\3 < M3 

1 1 A(W^, £)w| I, < C(M3)(1 + 1 1 W^l |,+2) I |u| l.+l , 

\\nu\\s < C.(M3){||$||2(1 + ||W^||.+2)|H|2 + IHI.+ l) + ||$||.+l|H|2. 

Therefore, the operators $ and A satisfy Conditions (B3) and (B4). Now taking 
into account the result of Theorem l5.9l it appears that Condition (B5) is satisfied 
for 

5 < p < r - 18 
cr = 3, g = 2, 

Finally, since iV > 3, condition (B6) is also satisfied thanks to Theorem 14.21 so 
we can apply Theorem 16. II and Theorem 1 1.31 is proved. 
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A Analytical study of 

A.l Computation of the differential of Qr] 

Let us make the following change of coordinate for 2:3 

S = X3- T]{X) 

and denote by the function defined by 

Then the Dirichlet-Neumann operator is defined by 

Ae-2WxV^x{^)-^AxV+j-^{^xvf^0, s<0 (A.l) 
OS OS OS'' 

S79 as s ^ —00, 

and 

g^^P = (1 + {VxTlf)^\s=o - Vx?7 • VxV-- (A.2) 
OS 

Notice that the above equations ljA.l|) . ljA.2|) may be written into the form 

V • (F,,V6') = 0, s < 0, 
Gn^ - (P,V0)-e3 

where 63 is the unit vertical vector, and Pr, is the following symmetric matrix 



-(Vx?7)* l + (Vx»7)^ 



We can see easily that the operator Q^, is symmetric and non negative in 
L^(R^/r) : for 77, ipi, ■02 smooth enough bi-periodic functions 

(S^^i,V'2) = ((F„V0i)- 63,^2) 

= [ [ (V6'2 • Pr,^ei)dXds 

J-00 JR2/r 



which is symmetric. Moreover, we have 



OS OS 

hence, for 77 and ip smooth enough bi-periodic functions 
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Let us now compute formally the differential of C/^. /^From (|A.2|I . we obtain 

d^,g^[h]ip = (p^vei) ■ 63 + iQ^,[h]v0) ■ eg (a.3) 

with 9 as above, and where 9i satisfies the system 

V-(F„V6ii) = -V- (Q^[/i]V6i), s<0 
Oils^o = 0, 
V^i ~> 0, s-^-oo, 

n \h]- f ^^^'^ 

Let us notice from HA.1() . 1A.2() . that we have 



V- (p^v(/i|^)) +v-(g„[/i]v0) = o, s<o. 



Indeed, we have 



hence 



V- [P^W{h^^)]+W-iQ,,[h]W0) 



do, (, do, do d^e\ 

-) + V./..|(P,V^)-Vx-+V.,^j, 



and this cancels, thanks to V • (P,jV6') — 0, and to the definition of P^. 
It results that 

OS 

V • (P,,V6i2) = 0, s < 0, 
do, 



V02 ^ 0, s 



72|s=0 — — 't^|s=0, 



hence looking at relationship (|A.3p . we obtain 

an 

= -Gvihg-Jo) + h{{l + (Vx^)')^|o - Vxr? • Vx^lo} 
OS 
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and using V • (Pj^VO) = 0, we get 

06 09 09 

d^g^[h]i; = -g^^{h—\o) + h{-Ax^P + Vxv'^x-g^\o + Axv-g^\o} 

06 



OS OS 
as required in Lemma [2. II thanks to (|A.2|I 



A. 2 Second order Taylor expansion of in = 

Let us consider the Taylor expansion of Q^i in " powers" of rj : 
= ^(0)^- + g^^^rj}^ + ^(2){^(2)}^ + ... 

where f/^'^^ is k— Hnear symmetric with respect to 77, and hnear in ^p. Moreover, 
for any fc > 1 and to > 2, t/^'^' is bounded from 

{iJ™+i(RVr)}'= into C{H^'+\R'^/T),H^{RyT)). 

From the lemma above and (|2.1I) . 12. 2() . we obtain 

g(i){?7}?A = -g(°)('7^^°^V') - V • (TyV^A). (A.4) 

Differentiating once more (|2.2|l with respect to ry, we now obtain in 

d^C[h]^-g^°\hg^''^)-hAi,, 

and differentiating H2.1|) with respect to 77, we then obtain the (symmetric) 
second order derivative in 0, and 

^ gm{r^gm{r^gmi,)) + ^g(°Hv'^i^) + ^Alry^tjWv,). (A.5) 
Now, if we Fourier expand any bi-periodic function as 

then we have 

= \K\^K, (A.6) 
{g^'HvmK= {{K-K^)-\K\\K,\}^K,VK,, (A.7) 

\^^{\K^ + X2I + \K, + K,\ - \K\ - \K^w^K,m.m,, 
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and we check that, for to > 2, the operators are bounded from {iJ'""'"^(R^/r)}'^ 
into £(i/™+i(MVr),i?^(KVr)) since we have 

\{Ki + K2) ■ Ki) - \Ki +K2\\Ki\\< MKi\\K2\, 

and there exists a constant c such that 

\Ki + K2 + KsWlKi + X2I + \R\ + A'al - \Ki + if 2 + K^] - \Ki\\ < 

<C\K2\\K3\. 

B Formal computation of 3-dimensional waves 

In taking Uq = (1,0), the symmetric hnearized operator for fi — and u = Uq 
reads 

d 

and £0 has a four-dimensional kernel, spanned by the vectors 

Co = (l,-z/Mc)e'^-^, Co = (l,VMc)e-'^^-^, 
Ci = Ci = (l,*/Mc)e-*''^-^. 

We observe that the action of different symmetries of the system on eigenvectors 
is as follows: 

'^iCo = Ci, SiC,i — Co- 
Let us write formally the nonlinear system \i.<6\ . (|1.7|l under the form 

CqU + fiCiU + C2{uj, U)+Af2{U, U)+Af3{U, U, U) + 0{\\U\\^ + \uj\\\U\\^) = 0, 

(B.l) 

with U — (7/1, rj), fl = fi — fic, w = u — Uo 

CiU = (0,7?), 

C2{UJ,U) = (-c^- V77,C^- VV'), 

■^3(f/,c/,c/)-| -^(v,-vv^) • 
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B.l Formal Fredholm alternative 

Let us consider the formal resolution of the linear system 

CoU = F={f,g), 

with 

U = Yl t/„e*("^^-^+"=^-^), C/„ = (Vn,^n), ^0 = 0, 

n=(ni,rt2)£Z2 

F = Yl F„e^("^^-^+"=^-^), Fn = {fn,9n), /o = 0. 

n=(ni,Ti2)GZ2 

Then, we have 

\niKi + n2K2\il)n - i{{niKi + n^K^) ■ uo}??„ = /„ 

i{{niKi + n2K2) ■ Ufi^^Jn + Hctln = 9n, 

hence for 

{{niKi + 722^2) • uo}2 - iJ.c\niKi + n2K2\ ^ 
i.e. by assumption for (ni,n2) 7^ (±1,0), (0, ±1), this leads to 

{iniKi+n2K2)-uoV- ^^c\nlK,+n2K2\' 
_ i{{niKi + ^2-^2) ■ upj/n - Ini^i + 712-^2 Ifln „s 
~ {(nii^i+n2i^2)-uo}2-Mc|nii^i+ 712^2! ' ^ 



and for (ni, 722) = (0, 0) 



V'0,0 = 0, 7/0,0 = —50,0, 



while for (rii, 712) = (±1, 0), (0, ±1), we need to satisfy the compatibility condi- 
tions 

{F,Co) = {F,Co) = {F,Ci) = {F,Ci) = 

which gives 

Mc/1,0 + m,o = 0, 

Mc/-i,o - ig-1,0 = 0, 

Mc/0,1 + igo,! = 0, 

Mc/0,-1 - igo-i = 0. 

For uniqueness of the definition of the pseudo- inverse Cq^, we fix U such that 

{u,Co) = iu,Co) = {u,Ci) = {u,c,) = o, 



78 



hence this leads to 



1 » , 



B.2 Bifurcation equation 



Now coming back to (|B.ip . we use a formal Lyapunov - Schmidt method and 
decompose U as foUows 

u = w + v 

(t/,Co) = {vxo) = {v,ci) = {vrci) = Q- 



We can solve formally with respect to V the part of equ. I|B.1(I orthogonal to 
the 4-dimensional kernel of >Coj a uniquely determined formal power series in 
oj, fl, A, A, B, B, which we write as 

V = V{il,Lu,A,A,B,B). 

The uniqueness of the series and the symmetries of our system lead to the 
following identities 

%Vifi,u;,A,A,B,B) = V{jl,uj,Ae'^^-^,Ae-'^'-^,Be'^^-^,Be-'^^-^), 
SoV{fi,uj,A,A,B,B) = V{fi,u;,-A,-A,-B,-B) 
SV{fi,0,A,A,B,B) = V{fi,0,B,B,A,A). 

The principal part of V is given by 

V - -C^'{l-PoW2{W,W) + 0{i\fA + \u;\)\\W\\ + \\W\n, 

-C^'M2iW,W) = A^U2ooo + \A\''Unoo + A^Uo2oo + ABUww + ABUoiio + 

+ABU1001 + ABUqioi + B^Uoo2o + l-Bpf/ooii + B [/0002 

where we observe easily that 

VoU2{W,W) = 0. 

Using (|A.7p . and 
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wc find 



M2{W, W) = ^2^2000 + |A| Vnoo + ^Vo2oo + ABV^ow + ABV^iw + 
-\-ABViQQi + ABVqiqx + i?^Vbo2o + |-BpVboii + B V0002 



with 



V2000 




V^llOO 


= 0, Fooii 


Vioio 


= 


Vbiio 


= (0,-2t2) 



= v^oioi — i^ioio 

i(-Ki-X + K2-X) _ tT 

Now, thanks to llET2|) . ||R3)| wc obtain -£0^2(1^, VP") as follows 

TT - ( * "'^ I „2i_fs:i X _ f ~* 1 A 2iK.2 X 

Unoa = 0, C/0011 = 0, C/0200 = C^20oo, ^^0002 = C^oo20, 

_ ^ -2»(2a., - 1) -2{fil + 2A.e - 2) ^ ^^(K,■X+K2■X) rj _Jj 

t/0110 = (^0,^)e^(-^-^+^-^), [/iooi = C7oiio. 

Replacing V by V(/i, >1, A, i3, B) in the compatibility conditions, i.e. the 
components of (|B.1I) on ker£oi 

(/i£i(VF + V) + £2(^^,VK + V) +AA2(H^ + V,T4^ + V) + ...,Co> = 0, 
(/x£i(iy + V) + /:2(t^,VF + V)+AA2(VK + V,VK + V) + ...,Ci) = 0, 

noticing that the complex conjugate equations are then automatically satisfied, 
lead to two complex equations of the form 



fill,iu,A,A,B,B) = 0, 
g{fi,uj,A,A,B,B) = 0, 



for which the equivariance of the system ljB.l(l with respect to various symme- 
tries leads to the following properties 

/(A, u;, Ae'^^'' ,Ae-'^'-^ ,Be''^^-^ ,Be-'^^-^) = e''^'-'' f{il, lu, A, A, B,B) 

f{ii,Lo,-A,-A,-B,-B) = -Jifl,Lj,A,A,B,B) 
g{fi,uj,-A,-A,-B,-B) = -g{il,u;, A,A, B,B) 
fill,0,B,B,A,A) = giil,0,A,A,B,B). 
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Since Ki and K2 are linearly independent, it results that / and g take formally 
the form 



g{ti,Lu,A,A,B,B) = B0i(A,c^, |^|Mb|2), 

where functions (f)o and 0i are real valued. Moreover, for a; = we have 

Mf^,0,\A\^,\B\') = Mf^,0,\B\^,\A\^). 

It results immediately that we have the following (formal) solutions of our sys- 
tem (in addition to the trivial solution 0): 

Proof, i) B — 0, \A\ satisfying (j)o{^,uj, \A\'^,0) = 0, which is not else that the 
2-dimensional travelling wave with basic wave vector Ki, and where, with no 
loss of generality, we can choose the velocity c in the direction of Ki. 

ii) A = 0, \B\ satisfying (j)i{fl,Lu,0,\B\'^) = 0, which is not else that the 
2-dimensional travelling wave with basic wave vector K2, and where, with no 
loss of generality, we can choose the velocity c in the direction of 7^2- 

iii) \A\ and \B\ such that 

0o(/i,t^,|A|Mi3n = 0, (B.5) 
0i(/l,o.,|A|MBn = 0, 

which gives a family of 3-dimensional travelling waves. Moreover we notice that 
if = there is a family of solutions where \A\ = \B\ and 

<^o(A,o,iAiMAn^o, 

representing the "diamond waves" of The leading terms of (j)o and (f>i are 
computed by Bridges, Dias, Menasce in |^, even in cases with a finite depth 
and with surface tension. Since our case has less parameters, our computations 
may look simpler. The leading terms independent of \A\ and \B\ in (j)Q come 
from 

iilCiCo + C2{uJ, Co), Co) = _ • 

and we notice that (in using ^ = —2-^) 

Mc Co 

We then have 

U~^,io,\A\\\B\'') = - (j^-2uj-K,] +ao\A\^+(3o\B\^ + Om + \Lu\ + \A\^ + \Bm 

Mc \^J■c ) 

<^i(A,c.,|A|MBn = 1 (}L^2uj-K^ +/3o|Ap + ao|i3p + 0{(|M| + k| + |^P + |Snn, 
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with 

ao = ^(2AA2(Co,C/2ooo) + 3AA3(Co,Co,Co):Co), 

l3o - ^(2AA2(Ci,C/iooi) + 2AA2(Ci,C/ioio)+6AA3(Co,Ci,Ci),Co). 
We can formally solve the system of equations 

Mt^,^AA\^AB\-') = 0, 

Mp,u;,\A\^,\B\^) = 0, 

with respect to fi, — 27^^' (^i ^^2)- Indeed, we obtain respectively in adding 
and subtracting the two equations, a system easy to solve, in taking into account 
of = 0(w|), 

A = -^{ao+mA\' + m + 0{i\A\^ + \Bm, 

C.2 = {\A\'-\B\^)[!^iao-Po) + 0{\A\' + \Bn). 

For the computation of coefficients ao and /3o we use again HA.7|) and (|B.4|I and 
obtain 

2A/'2(Co,C^2ooo) = (^^,o)e^^-^, 
2AA2(Ci,C/iooi) = (^-^,0^e*^-^, 

^ f 4(/^g + 4A^g - 5a^c + 1) Sij^a - 1)(1 - ^g) ^ ^.K^-X 

Now, with HA.8|I and the form of A/a we have 

3A(3(Co,Co,Co) = f^'^V^'""''' 

6AA3(Co,Ci,Ci) - f^(2-l),^(2-i,))e^^-^. 



Finally we obtain 



/3o = 



4 

12 32 8 36 



Mc Mc A^c Mc Mc M?(2-Mc)' 
If ao + /3o and ao — Po are both different from 0, we can solve the system (jB.5|l 
with respect to and and obtain a formal expansion in powers of (/i, w) 
of the form 



-/ij(ao+/5o) Aic(ao-Po) 
\B\ - 57 — — 7 --^+C( (^,t^j j. 
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The indeterminacy on the phases of A and B means that we obtain in fact, 
for each fixed leading to positive expressions for and \B\'^, a torus 

of solutions, which is generated by acting the operator Tv on a particular so- 
lution, for instance with A and B pure imaginary. This two-parameter family 
of tori of 3-dimensional waves connects with the 2-dimensional travelling waves 
respectively of wave vectors Ki and K2 (choosing oj orthogonal to K2 or to Ki. 
If we choose a; = which means that we choose the direction of the waves as 
Xi axis, then we obtain the "diamond waves" as in j36j, and 0, here without 
surface tension. 

Let us study the sign of (ao -I- /3o) and (ao — Po)- We have 

ao + Po = ^ ( + + — - 8 - 2^c - 
K \ K K Mc 

4 / 3 2 3 „ „ 

Po-ao = [ T + ^H 8-2^c 



2-^1, 
9 



where we notice that 



Aic I Aic Mc /^c 2(2 - ^c) 

and it is easy to show (study of the factor as a function of l//ic) that ao + Po > 
for r G (0, Tc), and a^+Po < for r > Tc where Tc G (\/6, a/7), i.e. more precisely 
A*c,c ~ 0.374. We also notice that the function of /ic in the factor for Pq — ao is 
strictly increasing for fic G (0, 1), and cancels for fic ^ 0.893. Now, defining new 
parameters ei and £2 by 

^~ 2'^- 2 ' 

for a particular choice of a solution belonging to the torus, we get the results of 
Theorem 10 ■ 



C Proof of Lemma 13.61 

Let first consider the unique solution of the Cauchy problem (notice that ^ g 
F^-i(RVr) C Co"V^(KVr), hence the vector field is Lipschitz for m > 4) 

dX V2, 
ax I Vi 

which we denote by X{xi,y) G C""~'^. We can successively show that X{xi,y) 
is even in xi, odd in y, and such that 

X{xi,y) = X{xi + 2n,y) = X{xi,y+—) - — , 

T T 

in using the uniqueness of the solution of the Cauchy problem, in looking for 
the system satisfied by X{—xi,y), —X{xi,~y), X{xi,y + ^) — which is 
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the same as X(xi,y), and in comparing the systems satisfied by X(xi + tt, y) 
and X{xi — tt, y), using the evenness in xi. Notice that 

dX ( V2 

-(a;i,y)=exp / d^^{ — ){t,y)dt 



equals 1 for ^ = 0, i.e. if [/ = 0. Hence, for \\U\\i smaU enough, we can solve 
(implicit function theorem) with respect to y = 2/(^2) the equation 

UiXi-,y)^Z2 

where Hi represents the average over a period in xi. Then, we set 

Z{Z)=X{z,,yiz2)) 
and, thanks to ^{xi + 7r,.X2 + 7) = ^{xi,X2) we observe that the function 

, TT TT 

Z(Zi + TT, Z2 + -) 

T T 

satisfies the same differential equation as Z{zi, z-i) with the same average in zi, 
hence by uniqueness it is identical to Z(z\^Z2). It is then clear that {z\^Z2) 1— > 
Zi — Z(Z) = di{Z) satisfies the properties indicated at Lemma l3 . 61 hence lies in 
C™o"^(IR^/r). The proof of the tame estimates for di and di is identical to the 
one made in Appendix G of [23 . Notice in addition, that replacing the initial 
condition at = by an average condition, allows to keep the equivariance 
under shifts Tg parallel to xi direction. 

D Proofs of Lemmas 13.71 and 13.81 

Let us start with the expressions for U, given in (|2.10() 

with ^0 — ("011 ^i) a-iid /^i(t) given in H2.10|l and H2.11|) . and where 

{ . ^^I^^" X sin 2xi — — !^ sin 2xi cos 2tx2 
^ Vc(2-Mc) ^ -^4/^2 i ^ 

4Mg(2-Mc) + ^ ^"^^2 + 4^ COS 2xi COS 2tX2 

= (V'2,?72)- 

We successively find from and H3.4|l 

b = — sina;i cosra;2 + e^b*^^ + O(e^), 



1 1 1 



2/x2 Afil A{2-^i, 
' - — 5- sin 2x1 cos 2Ta;2 , 



sin 2xi 
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Vi = 1 + ecosxi cosrx2 +e^F/^^ + O(e^), 
V2 = — er sinzi sinTa;2 + e^Vj'^'' + 0(£^), 



1 

4^: 



.(2) 1 



F]^^ = — - — + cos 2Ta;2 + cos 2a;i cos 2Ta;2) 



1 1 3 



V? 4(2 - Mc) 



cos 2xi, 



r 

= - — ^ sin 2x1 sin 2Ta;2, 



a = + — cosxi C0STX2 + e^a^^^ + O(e^), 



-— ^(1 — cos2xi cos2rx2) H — 



{y2 ^ v^)2|l/2 ^ ^ ^ ^^^^^^^ ^^^^^^ ^ ^^^^^^ ^ ^ 



..f2i 3 5 /7 1 3 3 ^ „ 



7 1 ^ X 

+ (- + -— ^) cos2Ta;2 + (77 + tt^) cos2a;i cos2Ta;2. 
8 2/^2 



For the determination of the function Z{Z) = Z2 — di {Z) of Lemma 13.61 we have 



V2 

= --^(a^i,a;2 - di), Ilidi = 0, 



I.e. 

~ ~ (1) o ~ (2) , 

with 
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This leads to 

r(i) 

di = — TCOSZi SinTZ2, 

~(2) T 

di = — COS 2zi sin 2tz2, 
4 

hence the diffeomorphisni U^^ takes the form 

xi = zi, X2 — Z2 + ercoszi sinrz2 — cos2zi sin2TZ2 + 0{e^). 

Now we have 

with 

dj^ = — rcosxi sinra;2, 

= j{T2 + (l + r2)cos2a;i}sin2Ta;2, 



T 

4 

and finally ()3.32|l gives 

q = /Ltc(l — ecosxi cosTa;2 + £^q'^^) + 0(e'^), 



with 



(2) 



— + - H — J + q2o^ cos 2xi + q^^^ cos 2tx2 + 



(2) 

+q22 COS 2xi cos 2Ta;2 , 



(2) 

q2o 



3 5 111 



4(2 -Mc) 8 4a*c 2/i2 ^3 fit' 

"^O^ - 8 + 4^2 ' ^22 - 8 + 4^2 ■ 

Now, we obtain 

{qoU^^){Z) = /io(l - £ cos zi cos rz2 +£^q^^^) + 0{e^) 

with 

q(2) = . +(qy+ )cos2zi + -cos2tz2 + 

/ic O 4/t2 4 8 

3 

+ - cos 2zi cos 2tz2, 
8 

and 

{{qoU-')(Z)y/' = ^i/2|l_£eosz,cosrz2 + £^(|;--| + ^) 

q(2) ^2 2^ 

+ y - ^)cos22i+£2— cos2rz2 + 

5 

+£2— cos2zi cos2tz2} + 0{e^), 
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hence 



+e^^cos2TZ2} + 0{e^). 



Then with H3.34|l . we obtain 



as indicated in Lemma l3.7l Furthermore, in using H3.33|l , H3.31|) . we notice that 







. , 




62 = 




sm 2rzo 4 

32r ^ 






e . 




d2 = 




— sm zi cos TZ2 






1 




1 


d2i = 


4 


{4(2-^,) + 


16 " 



3 111 , 5 

£122 



,2 ' ,,3 ,,4 
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This, with (|3.27|) . leads to the principal part of the bi-periodic functions occur- 
ring in the diffeomorphism of the torus of Theorem 13. 41 

d{x,y) = — - sina;i cosTa;2 + sin2a;i{(i2i 3- + ('^22 + -3-) cos 2x0:2} + 0(e 



. „ r 1 2 9 r2 



e(x, y) = ~ST cos Xi sin rx2 + — sin 2Ta;2 \ — 2~s~' — 2 '^'^^ r + 



Inverting ZY2 o Z^i, and changing the coordinates in 



T]{X) = — — cosxi cosTa;2 + £^ I ^-^-^3^— — ^cos2xi + -^ — cos2ra;2 
Mc L 4^^ (2 - ^c) 'ifJ-c 

1 



cos 2xi cos 2Ta;2 r' + 0(e^) 

4a*c 

leads to the results of Lemma 

E Distribution of numbers {uon^} 

Recall that for each a; G M, 

[x] = max{7V : iV e N, N<x}, {x} = x - [x] e [0,1). 
In this section we consider the distribution of the numbers 

On = {won^ - C'}, an = {loou}, n > 1 (E.l) 
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Due the famous Weil Theorem 0^1; for each polynomial f{x) = akx'' + .. + ao 
with irrational ak, the numbers {f{n)} are uniformly distributed in [0, 1]. This 
means that for each [a, (3] C [0, 1], 

card {n:0<n<N, {/(n)} e [a, /?]} = N{f3 - a) + o{N). 

This result is the asymptotic relation in which the remainder strongly depends 
on the choice of the interval. Wc consider the simplest case 

fix) = Loox^ ~ C with {fin)} = 0„, [a, /?] = [0, p] 

and deduce the rough, uniform in p estimate which is sufficient for our needs. 
The main result is the following . 



Proposition E.l Suppose that 
1 



< Cil for all positive integers I. 



Then there is a constant c depending only on c\ such that 

i l<cp for all pe (0,1/4) and N > p-"^^ . (E.2) 

l<n<Ar,e„G[0,p] 

The proof is based on the following lemma. 

Lemma E.2 Under the assumptions of Provosition \Ej\ there exists positive c 
depending on ci only such that 

^ 1 < ce for all ee (0,1/4) and N > . (E.3) 

l<n<JV,CT„G[0,e]U[l-e,l] 

Proof. Fix an arbitrary positive e G (0,1/4) and introduce the function de- 
pending on parameter e and given by the equalities 

iPe{x) = 1 for X e [0, e] U [1 - e, 1], 
ip.ix) = for X e [2e,l - 2e], 
2e — X 

ipeix) = for X e [e, 2e], 

e 

, . X - 1 + 2e 
feix) — tor X e 1 — 2e, 1 — ej. 

£ 

We will assume that (pe is extended 1-periodically onto R. Obviously the ex- 
tended function is absolutely continuous and 



ipeix)dx = 3e, / \ip'^ix)fdx = 2e-\ (E.4) 
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It has the representation 

oo 

( \ \ ^ 







It is clear that 



Obviously 



1 < y </?e(a;)(ix = 3e. 





1 



l<n<Ar,cr„e[0,e]U[l-£,l] l<n<Ar 

Represent (ys in the form 



(E.5) 



(E.6) 



^,(x) = (p,,o+ <^.,;e'"''" + QfcW, Qfc(a:)= ^ ^e.ze^^''" (E.7) 

i<|i|<fc fc+i<|i| 

We have, by equahty HE.4(I . 



fc+l<i 



fc+i</ v ; 



*;+!<; 



\ J vfce 
\ 



(E.8) 



Here c is some absolute constant. Combining (|E.6|) - (jE.8|l we obtain 



1 

TV 



l<n<JV,cr„e[0,e]u[l-e,l] 



l<Z<fc 



TV 



l<n<JV 



ke 



(E.9) 



Since the numbers 
form a geometric progression, we have 



E 

l<r!,<7V 



< 



a'2TZiljJ(-\l 



1| 



< 2ciV 



Substituting this inequality into (|E.9|) we obtain 



1 

TV 



l<n<A',(T„e[0,e]U[l-e,l] 



c cefc'^ 

< 3£ 



l</<fc 



ke 



N Vke 



(E.IO) 
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Now set 

and note that for N > e^^, 
1 



7k- (1 - £3)1/2 N 



which along with ljE.10|) yields 



1 

N 



l<ce for N > s-^, 

l<n<JV,CT„e[0,e]U [1-6,1] 



and the lemma \EI2\ follows. ■ 
Proof of Proposition IE. n 

The proof in main part imitates the proof of the Weil Theorem. Fix p S 
(0, 1/4) and consider the function as defined above. Obviously 



l<n<JV,6(„e[0,p] l<n<Af 

Combining (jE.7|l - (jE.8|) we obtain 



(E.ll) 



i<n<N.,e„e[o,p] 



KKk 



E 

l<n<Ar 



+ 



Next set 



Wi, 



N 



l<n<N 



We have from IE. 1211 



1 E i<sp+6pj2 j^w, 



N 



i<ji<A',e„e[o,p] 



i<i<fc 



kp 



(E.12) 



(E.13) 



(E.14) 



Noting that 
we obtain 

Next we have 



<n<N 



(E.15) 



l<m.n<N 



90 



Setting r = m + n, q = n~m we arrive at 



2<r<N \q\<r 

Introduce the quantities 



Xrl 



kl<'- 



^2-KiujQlrq 



E E 

N<r<2N |g|<2Ar-r 



for 2 < r < A^, 



Xrl 



\q\<2N-7- 



for N <r <2N 



Thus we get 



Obviously 



W'n^ E 

2<r<2N 



Xrl < ^N. 

On the other hand, since Xri is a geometric progression in q, 

1 



(E.16) 
(E.17) 



Xrl 



< 



smiirujorl) 

Choose an arbitrary e E (0, 1/4) and denote by J; the set of all r such that 

2 < r < 27V, a^i = {luoH} G [0, e] U [1 - e, 1]. 
It is easy to see that for r £ [2, N] \ Ji, 



(E.18) 



1 c 

Xrl < T < -• 

sin(7re) e 



From this and (IE.16II. (IE. 1711 we conclude that 



Wl^ < 27V ^ 1 + c— 

re.h 



(E.19) 



(E.20) 



On the other hand for fixed I, r G Ji, and p = rl we have Up = {i^op} G 
[0, e] U [1 — e, 1]. Hence, since I < k, 

card Ji < card {p : 1 < p < 2kN, ap e [0, e] U [1 - e, 1]}. 

By lemma Ie!^ we have 

^ card {p:l<p< 2kN, ap G [0, e] U [1 - £, 1]} < ce for kN > e^^, 



2kN 
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which gives 

^ 1 < c2kNe for kN > e" 

re.] I 



Substituting this incquahty in (|E.20p we obtain 



W^fjv < + c— for kN > e'^, 



or 

— Wi N < cVke + for kN > e" 

N ' y/Ne 



Substituting this resuh in (|E.14|I we finaUy obtain 



cpk c 



for kN > e 



4 y l<3p + cpVfc%+ . , ^ 

l<ri<A',0„e[O,p] ^ ^ 

It follows from this 

^ V 1 < 3p + cp\/fc3e + cpfc^/ V + for fciV > e^^ (E.21) 

l<n<Ar,6l„e[0,p] ^ ^ 

Now choose 



Obviously 



fc>^>32, e<l/4, Nk>e~^, k^'^e'^ <l, k^e^l. 



From this and (|E.21|I we conclude that 

^ E 1 < CP for iV > 

l<n<Ar,6»„6[0,p] 

which completes the proof of the Proposition. 

F PseudodifFerential operators 

In this section we collect basic facts from the theory of pseudodifferential op- 
erators. We refer to the pioneering paper and monographs 05]' ESI for 
general theory. Note only that different maps from functions A{Y, k) to opera- 
tors 21 = A{Y, —idy) give rise to different theories of pseudodifferential calculus. 
In these notes we assume that F is a coordinate on the 2£'-torus and that 
the dual variable k belongs to the lattice I?. The first result constitutes the 
continuity properties of general pseudodifferential operators. 



92 



Proposition F.l Let |21|q ; < oo and 0<s<l~3,r + s>0. Then there is a 
constant c depending on s only so that for all u £ iJ''+''(R^/r), 

\\%u\U<c{\^lM\r+mUMr+s)- (F.l) 

The proof is based on the foUowing estimate of the convolution of non- 
negative sequences. Let us consider a non- negative sequences — (a-* (ft))„gz2 , 
1 < j < and V = (w(n))„gx2. Set 

|a^|, = sup(l + H)V(n), |||v|||2= + 

Lemma F.2 Under the above assumptions, the convolution w = a"'^ * ... * a™ * v 
has the hound 

IIHII. <c(s)^(n|a^|3)|a^U+3|||v|||o + (nia'l3)llMI|.. (F.2) 

Proof. We begin with proving ljF.2|) for m = 1. RecaUing the formula 
w{n) = ^ a^(fci)..a"(fc„i)w(fcm+i) 

fci + ...+A;,„ + i=Ti 

and noting that for ki + k2 = n, 

{i + \n\y<cis){{i + \k^\r + {i + \k,\r) 

we obtain for m = I, 

il + \n\yw{n) < c{s)\a%+, ^ (l + \ki\)-'vik2) + 

ki+k2=n 

+C{S)\2}\, (l + |fc2|)Mfc2)(l + |fcl|)-^ 

ki+k2=n 

From this and the classic inequality 

E( E Kfci)^'(fc2)i)'< (EKfc)i)'Ei"(^)i'- 

n k\+k2—n k k 

we obtain in the case m = 1. The general case obviously follows from the 

mathematical induction principle and the distributive property of the convolu- 
tion. ■ 

Let us turn to the proof of Proposition lF.il We have 
^{^) ^ Mn-k, k)u{k) where A{p, k) ^ ^ [ A{Y, k)e-'^P dY, 
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which yields 

\^{n)\ < cY.lM'n- k,k)\\u{k)\ < [a*v](n), (F.3) 
fc 

with 

a(n) = sup (1 + \k\y\A{n, k)\ v{k) = (1 + \k\y\u{k)\. 
It is easy to see that 

|a|. < c|2l|5,„ |||v|||, <c||u||,+,. (F.4) 
Applying Lemma FF. 21 to ljF.3|) . using inequalities (jF.4p and noting that 



2lw||J <c^(l + |n|)2^|aMi2 



we obtain IjF.ip . and the proposition follows. 

The next proposition gives the representation for the composition and com- 
mutators of pseudodifferential operators 

Proposition F.3 Let 21 and 05 be pseudodifferential operators so that for some 
r, p eM.^ and non-negative integers to, I, 

+ < oo. 

Let also 

I > \r\+5 + s, l>s + 3, m>2. 

Then the composition 215B and the commutator 2125 — 5B21 have the representa- 
tions 



2tQ5-^(2tQ5)p + sl^f , d = 0,l, (F.5) 

d 

2l*B-'82l = ^[2l,Q5]p + D|;^fl, d = 0,l, (F.6) 



in which (2tQ5)p and [21, 25]p are the pseudodifferential operators with symbols 

{AB)o{Y,k)=A{Y,k)B{Y,k), iAB)i{Y,k) = -dkA{Y,k)dYB{Y,k), (F.7) 

[A,BUY,k)^Q, [A,Bl{Y,k) = ^{dkA(Y,k)dYB{Y,k) ~ dYA{Y,k)dkB{Y,k)) 

(F.8) 

for k^O, and {AB)p{Y, 0) = [A, B]p{Y, 0) = 0. 
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The reminders have the estimates 

|2l|d+i,3l»l^+i,|.|+d+4l^*ll.'+r+P-<i-i, (F.9) 

l|S)L+fi z^lU < c(|ar,+i,,+|,|+,+4l»l2+i.H+,+4+ 

|2l|d+l,|p|+<i+4|l - ^\d+l,\r\+d+4+s) \\u\\r+p-d-l + 
|2l|rf+ljp|+rf+4|l - "Sld+Url+d+Mls+r+p-d-l, (F.IO) 

in which the constant c depends on s,r, p only. 
The proof is based on the foUowing lemma 
Lemma F.4 Let 

^ 1 ^ 

7^d+l(r/,c,fc) = ^(r^C + fc)-E E ;;i77TH^^"^('^'^)(*'^)" 

(F.ll) 



A{r], k) = ^ I e-'^'^^A{Y, k)dY. 



and TZd+i{ri, C, 0) = A{ri, C), where 
Then for all rj,Q,k E 'i? and < s < I, 

\n,+,iv,c,k)\ < c(d,r)|A|5+i,(i + + |c|)i'-i+'^+i(i + \k\y-'-\ (F.12) 

Proof. It suffices to prove ijO^ for fc 7^ only. If |C| < \k\/2, then the Taylor 
formula 

d+1 ^ 



7^rf+l(77,c,fc)= [d^A]{7j,k + tc){i-trdt}c 

implies the estimate 

1 

I7^rf+l(r?, c, fc)i < cm'd+iM + ivir^a + \c\r^' /a + ifc + tciy-^-'dt < 





cmwiM + + \c\r+\i + \k\y-'-\ (F.13) 

which obviously yields H4.32|l . If |C| > \k\/2, we have 

d 

l7^,+l(c, fc)i < c^WiA^ + \\k + cr + E(i + icififcr"''' 



< 
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Noting that for < p < d + 1 and |C| > |fc|/2, 

(1 + icinfcr^ < c(i + icD'+^i + ifcir'-' 

we obtain 

\n,MCk)\ < c{s)mWiA^ + + \c\P+'+\i + \k\y-'-\ 

and the lemma follows. ■ 

Let us turn to the proof of the proposition. Since [21,05] — —[21, (1 — 5B)], 
it suffices to prove ()4.33|l only. To this end note that, by the definition of 
pseudodifferential operator, 

2tS^t(n) = ^ l(n - p,p)B{p - k, k)u{k) = ^ l(n - fc - C, fc + C)5(C, k)u{k). 

p,k {^,k 

Applying Lemma l-b'.4l to the Fourier transform A{'r], k) of the symbol A we arrive 
at the identity 



E EE; 

fc / 



1 



[d^A]{n-k-C,k) {tCrSiCk) u{k)+ 



+ J2 7^d+l(n-fc-C,C,fc)-S(C,fc)M(fc) =2l»^i(n). 

Noting that (zC)"S(C, k) = d^iC, k), we obtain 

E E ^^mA]in-k-0[iiCrB{C,k)]uik)^i^,uin), 
fc / 

which leads to representation H4.33|l with the remainder 



,(n)= ^ Ud+iin- k - C,C,k)B{C,k)u{k). 



In particular, we have the inequality 



[sL'fi^^.l <a*b*v. 



1+1 ' 

in which the elements of the sequences a, b, v are given by 

)-'^+'^+l|7^rf+l(n,C,fc)|, 



(F.14) 



a(7i) = sup(l + |C|)-l''l-''-^(H 
HO = (1 + |C|)l''l+''+^ sup(l + |fc|)-''|S(C, A:)|,t;(fc) = (1 + |fc|)'-+''-'^-i|?^(fc)|. 

k 
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Applying Lemma rF.2l to the right side of ljF.14|l and noting that by Lemma FF. 41 

|aU < |2t|S+i,,, \Hs < m'o,s+\r\+d+V \M\s<c\\u\\p+r-d-l + s 

we obtain H4.41|l and the proposition follows. 

It is useful to reformulate the above results in terms of infinite matrices. 
To this end we introduce the Hilbert space F^^i which consists of all operators 
y : 7J'*(R2/r) 1-^ i7*(R2/r) having the representation 

pal? 

such that 

Pill., sup { ^(1 + \k\f'{j^ |3^fcpPb)l)' : + < ^- 

k p k 

(F.15) 

Corollary F.5 (i) Under the assumptions of Provosition ^FAl operator^ has 
a matrix representation with Sl^p ~ A{k — p,p) and 

E(i + i^i)''(Ei2i'opii"(p)i)'< 

fc p 

cmiif E(i + i^D'i^wi' + c(i2ii5.3)' E(i + ifci)''^+'i^(fc)i'- 

fc k 

In particular, ||2t||F,,+s . <c|2l|g;. 

(ii) Under the assumptions of Provosition VF.SI the operator Ti'^J^^ has a 
matrix representation so that 

T.(^ + \k\f'^{T.\^^f+Zp\^^^^^^^ 

k p 

i2iis+i,3i»r.+i.M+.+4+s)'E(i + i^i)'^^""''"'i"Wi'+ 

fc 

c(i2iis+i,3i»i^+i.i.i+<i+4)' E(i + \k\r+''^^'''-"'-'\m\' (F.16) 

k 

In particular, 

||S)i+f ||f.,.,,_,_,,. < c(|a|5+i,J»l2+i,|,|+,+4 + |2i|5+i,3|si:;^, 

Proof. Assertion (i) is integral part of the proof of ProDOsition lF.il In order 
to prove (ii) note that 

» i+tfcp = E - ^ - C, C, k)BiC k), 
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where TZd+i is defined by formula (|F.11|I . The needed result follows from (|F.14|I 
and Lemma FF. 21 ■ 

Proof of Proposition 15.41 

We give the proof of representations (|5.15(l and H5.1fi|l only, and begin with 
proving (ISmi). Note that 



where 
A 



Assuming fci 7^ and noting that in this case ^| 1 — we arrive at the 
identity 

A{Y,0 = MY,^i)+t^2A{Y,^,) = A(Y,0 where Ap{Y,^i) = Ap{Y,C,, 1-^?). 
From the Taylor formula we conclude that for fci 7^ 0, 

J=2 



AiY,0^J2 



1 



3=0 



dlA 



(F.17) 



7^3(^,0 = IJ [dlA]{Y,s^,,^2){l-sfds. 



Recalling symmetry property (|5.2|l we obtain 



Ar + i^2A 
d^,{Ar + i^2Ai) 



?l=0 

= i 



Re A + i^2^m A 



(Y,0A) = Ao(Y,^2), 



ei=o 



d^, (Im A - i^Rc A) (Y, 0, 1) = -lyAi (Y, 6) 



dl{Ar + t^2A^) 



(<9|i - c^cJlRe A + i^slm A) + i^alm A (Y, 0, 1) 



which being substituted into ljF.17|l leads to 

AiY.o = Ao{Y,^2) - i.Ai(y,e2Ki + i.'A2(r,6)(*6)' + 7^3(r,c)e^ (f.is) 

Noting that 

1 1 



V v\T-^(k)\ 



L{k), 



z/2|T-i(fc)| J.2|T-l(fc)|3 



{vkl + \T-Hk)\)L{k), 
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we get for ki ^ 0, 

2 

thA{Y, + Pa{Y, Om + Qa{Y, 0- (F.20) 

Here the remainders are given by 

PA(Y,k) = + J^;f^{MY,^^)- (F.21) 

Next set P^(r,fc) = QA{Y,k) = for fci = and denote by 0a the 
pseudodifFerential operators with the symbols Pa, Qa- With this notation, 
decomposition H5.15(l easy foUows from l|F.20|l . It remains to note that formulae 
l)F.17|l and 1F.21|I imply the estimate 

\^a\-,} + \qa\-,} < m3,i, 

which along with Proposition IF . II yields inequalities (|5.2l)|l for and £Ia- 

Let us turn to the proof of (|5.16() . We begin with the observation that for 
fci 7^ 0, the product of the symbols of the elementary operators 21^ and 2IJ is 
equal to 

AjW = S- + Cilm ijim W, (F.22) 

where 5*^ are the symbols of the elementary operators associated with the func- 
tions AjW. Multiplying both sides of ljF.20p by W{Y, ^2) and using the identities 

^1 iy\r-\k)\ j/|T-i(fc)|2 ^ ' 

we arrive at 

2 



ikiS{Y,k) ^y^iiki)^-^ S'j{Y,^2)-(-) lmAolmW—+L{k)Us{Y,k)+Vs{Y,k), 



where the symbols Us{Y, k) and VsiY, k) are defined by 



Us = PaW + (^)'^p4^(-fc? + (fc)|)lm iolm W- 
^|T-i(fc)|2 (l'" ^1 + (*^i)''lni i2)lm 
Vs = ^|T-i(fc)| (li" ^1 + (*^i)"'lm i2)lni W. 
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Recall that W and A are smooth function on tori which do not depend on 
k. Setting Us{Y,k) = VsiY,k) = for ki = 0, denoting by lis and QJs the 
pseudodifferential operators with the symbols Us, Vs and arguing as before we 
obtain desired identity H5.15|l . It remains to note that estimate (|5.21ll follows 
from Proposition (|F.1|I which completes the proof. 

G Dirichlet-Neuman operator 

In this section we deduce the basic decomposition for the Dirichlet-Neumann 
operator and prove Theorem 13. 51 Let us denote the change of coordinates by 

X = x{y), x=(X,X3), y^{Y,y3), 
X - X{Y), X3^y3 + fi{Y), (G.l) 

where 

77(r) = viXiY)), 

and X{-) is a diffeomorphism of the torus of the form H3.8|) satisfying the Condi- 
tion Let us notice that with the new coordinate Y = (2/1,2/2) the lattice F' 
is generated by the two wave vectors (1, ±1). We still denote the lattice of peri- 
ods by r. For a function u{x) we define u{y) by u{y) = u{x{y)). The Jacobian 
matrix of the above diffeomorphism reads 

/ dy.Xl dy,Xl \ 

^liY) = dy,X2 dy,X2 

\ dy.fl dy^f] 1 / 

and the determinant satisfies J = detB(F) — detBi(y). Now, we use the 
following identities for any scalar function u, and vector function V : 

W,-V{x{y)) = jWy{JM^'V{y)). 

With these identities, the Dirichlet-Neumann operator (|1.4|) takes the new fol- 
lowing form 

A(f = 0, 2/36 (-00,0), 
"SKp as 2/3 ^ —00, 

g^^j = i(A(F)V^(^).e3, (G.2) 

where 

= V,y • (A(y)Vy(^), 
A J(B^Bi)^"'" (symmetric matrix) 
det A = J. 
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Notice that for computing the new expression of Griip, we used the fact that 
$(x) = ^3 - f]{X) is such that $(y) = 2/3, hence = (BJ)~^e3. 

We aheady defined the 2x2 matrix Gr{Y) of the first fundamental form of 
the free surface, and we have 

/ 511 912 dy,7j \ 
B*iBi(y)= 512^ 322^ dy,?j . 

\ dyiV dy^V 1 / 

We assume that elements of the matrix A{Y) are smooth 27r-periodic functions, 
and for some p and Z, satisfying the following inequalities for 9 < p < I, 

||Ao - A||cp < c£, \\A\\ci<Ei, where Ao = T~Miag {1,T^ 1}, (G.3) 

where, by construction, the estimates of A and Aq — A would come from 

ll^llcp+i + l|V(-)llcp+i < ce, (G.4) 
lhllc=+i + l|V(-)llc=+i < c(0, s<l, 

where V is defined in (|3.8|) . By the factorization theorem, there are first order 
pseudodifferential operators so that 

A = a33{dy,+g+){dy,+g-), (G.5) 

^-^k■Yg±^^k■Y ^ |^| _^ 

It follows from (jG.2|l that 

Gjji^ ^ ji^a^jdy. - a^sG )tp. (G.6) 
i=i 

Hence the task now is to split G~ into a sum of first and zero order pseudodif- 
ferential operators. The corresponding result is given by the following 

Theorem G.l Under the above assumptions there is Eq depending on p and I 
only such that for all e € (0,eo); the operator has the representation 

g- = g^ + go +Qz^, (g.?) 

in which the pseudodifferential operators g^ ,gi , have the symbols Gq and 
defined by 

Go = g ( zVfcG+ VyGr - 26Gr + Gi) , (G.8) 

Gf{Y,k) = — (za3ifci+za32fe±£i), (G.9) 
033 
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where 

D{Y,k)^l 2^ assajmkjkm - {asiki + a32k2) | , 

l<j.m<2 

.2 ^2 (G.IO) 

C\{Y,k)^ — {dy^ajm)km, 26= — dyMsj. 

0-33 . , 033 , 

j,m=l ] = l 

The zero-order pseudodifferential operator satisfies the inequality 

|eo"lm,„-i <c||A-Ao||c", m>0, n<l, (G.ll) 

and for any u G H^^^{M.^ /T), r < p — 8 and 1 < s < I ~ 8, the rest term has 
the bound 

WGZMlr < ce\\u\\r-i, \\gi^u\\, <c{Ei\\u\\r-i+e\\u\\s-i). (G.12) 
Moreover, operators , ^ QZi verify the following symmetry properties 

g-u{±Y*) = g-u*{±Y), J -1,0,-1, u*{Y)^u{Y*). (G.13) 
Proof. First we rewrite the operator A in the form 

^ = 033773 + 2a33T^S + a33C, C = C2+Ci, B = B, + b (G.14) 
dyi dy3 

with differential operators 

2 2 
^2 = asi ^2 o-ijdyidy., Ci = 033^ ^2 i9vi°-ij)'^vj^ 



2 2 

^1 = «■« "3j^w ' 26 = 033^ dy^ asj . 



(G.15) 



Combining (jG.5|) and ljG.14|l we obtain 

g+ + g- = 2B, g+g-=c. (g.i6) 

We find the solution of ljG.16|l in the form 

g^ =g^ +g^±x, (g.i7) 

with 

gt + Qi='2Bi, g+ + g^ = 2b, (g.is) 

where the symbol of pseudodifferential operator g^ is given by formula (|G.9p , 
gQ satisfying l)G.ll(l and X being unknown. It follows from these formulae that 

g± = ±(-A)i/2 + ^± 
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and for any integers m > and n < I 

\Gi\ln,n + |5o"l™,„-i + \l3L,n-i < c||A - Ao||c". (G.19) 
Representation l|F.5p from Proposition IF .31 yields the decompositions 

g^g^ = a^/ + 7^ol, g^go - gil + ^lo- 

Here the second order pseudodifferential operator g^^'' have the symbol G^G^ , 
and satisfies the identity 

ef^-C2, (G.21) 
the first order pseudodifferential operators g[^\ ^^io\ and g^^^^ have the symbols 



(G.22) 



g[^^ =-idkG+dYG^, 

"-^10 — "-^1 "-^0 ' ^01 — "-^O "-^1 ' 

and satisfy the identity 

g['^ + g^u)^ + goi - Ci = o. (G.23) 

Now from (|G.16|) and HG.21|I we have 

Gt + G^=2Bi, G+G^=C2 

which leads to 

Gf ^Bi± {Bl - C^fl^ 
and since 033 > for e small enough, and noticing from l(niO)l that 



033 

the formula (|G.9|) follows. For obtaining Gq", we use (|G.23|I . (IG.22p which leads 
to (ITTSll 

It is then clear that G^ is a pseudodifferential operator of zero order which 
satisfies (jG.lip . Now, the symmetry properties ljG.13|l for g^ and g^ follow 
from (|G.9|) . ljG.8|l . the evenness of a, 033, an, 022 in y\ and 2/2, the oddness of 
h and a\2 in 2/1 and ?/2j and the evenness in yi, oddness in 2/2 of ai3, and the 
oddness in yi, evenness in 2/2 of 023. 

Inequality ljF.91) from Proposition lF.3l alone: with (|G.19|) implies the estimates 

ll7^ol^i|U + ||7^lo^i|U + ||7^lM|U < c||A- Ao||c6+^-||A- Ao||c6||w||o + 

+c||A- Ao||^6||u||.. (G.24) 
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Substituting identities HG.20|) . (|G.21|I and HG.23|) into HG.16|I gives the equation 
for the operator X 

X'^ + {-AY^^X + A'(-A)i/2 +UX + XV + W = 0, (G.25) 

where 

w = -{111 + + 7^ol + g^g^)- 

Our task is to prove the existence of a "small" solution X to ljC25|l . 
Introduce the Banach spaces of bounded operators 

= Fs_i^s n F<;,s+i, Y s = Fs^s_i, Zs = Fs^s 

supplemented with the norms 

IIA-Ilx. = ||A'||f.,.+, + IIA-IIf.^,., IIJ^IIy. = ||3^||f._,, ||2||z. - ||z||f.„. 

It follows from ^J^, CorollaryESl and (jU^M that for aU r < p - 8, 

\m\^r+i + I|V||y„ + ||W||z. < c£, (G.26) 

where the constant c depends on I and p only. The rest of the proof is based on 
the following lemma which is proved at the end of this section. 

Lemma G.2 Under the above assumptions, there exist eo > and c > de- 
pending on p and I only such that for e € (0,eo) equation (|G.25|) has a solution 
satisfying the inequalities 

llA'llx, <ce whenQ<r<p-8, 
\\Xu\\s <c{e\\u\\s-i + Ei\\u\\o) when l<s<l~8. (G.27) 

Moreover the operator X satisfies the following symmetry property 

Xu{±Y*) = Xu*{±Y). 

Proof. We start with the consideration of the simple linear operator equation 

{-A)i/2^ + A'(-A)i/2 = Z with ZeZs. (G.28) 

It is easy to see that for any Z e Zs satisfying ljG.28|l and having the matrix form 
with elements Zkp , operator X has the matrix representation with elements 

Xkp = jv(fc) + jV(p) ^^^^' ^^^^ " + ^'^2'^'^' 

which obviously yields the estimate 

IIA-Ilx, <c(r)||Z||z.. (G.29) 
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Hence the mapping Z ^ X defines a bounded linear operator S G £(Zs,Xs). 
Let us consider the sequence of operators defined by the equahties 



A-o = 0, A'^.+i = -S A-^f + UXr, + Xry + W 



Note that is a pseudodifferential operator which symbol U{Y, k) satisfying the 
inequalities 

\\U{;k)\\cr <ce, \Uip,q)\<csil + \p\)-''+'\q\. 

It is easy to see that for any X G X^, the operator UX has a matrix represen- 
tation with the matrix elements 

iUX)kp = U{k-q,q)X^p. 



We have 

(1 + \k\y\\{ux)kp\ < ceJ2{i + \k- qiy-p+'ix.piii + \q\y+\ 

q 

which gives 

^(l + |A.|f'-[5]|(ZYA'),,|Rp)|]'< 

k p 
k p q 

E [E(i + 1'? - ^1)'^""^' ((1 + \i\y'^' E \^qpMP)\ 



< 



k q 



< 



cs' ^ [(1 + \q\Y+' Yl \^qpMp)\\ ' < c^'II'^IIf...,, E(i + bD'I^WI' 



Thus we get 

\\UX\\z^ <c||Z^||y.^J|A'||x„. 
Repeating these arguments we obtain 

IIA-Vllz. <c||V||y„||A'||x„, 

and 

\\X'h.<c\\X\\l,<c\\Xf^^, 
inequality ljG.29|l yields the estimates 

||^«+i||x. < cedlA-^llx,, + 1) + \\Xn\\ji^, 
which holds true for all r G (1, p — 8). On the other hand, since 

'^Ji+l ^ '^ri " i'^n+l — Xn)Xn^i + Xn{Xn+l ~ Xn) 



(G.30) 
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we have 

\\X.n+i - A-^Jlx. < C\\X,, - Xn-l\WM + \\Xn\W^ + \\X,,-l\U^)- (G.31) 



Here the constant c depends on g and t only. It follows from ljG!.30|) that for all 
e G (0, £o(Pi t), the values HAViHx^ are less than ce. From this and (|G.31|) we 
conclude that for all small e the sequence Xn converges in X^. Repeating these 
arguments and using Corollary IF. 51 gives the tame estimate (|G.27|) , and the 
lemma follows once we observe that the symmetry property of X follows from 
the uniqueness of X and from the equivariance with respect to the required 
symmetry of the equation ljG.25p . ■ 

In order to complete the proof of Theorem IG.ll it remains to note that 
operator QZi := —X with X given by Lemma FG . 21 satisfies (jG.12p . ■ 

Proof of Theorem It follows from formulae (|nj|) . and ((^7)1 

that we can write ^ 

where the first order pseudodifferential operator T) has the symbol D{Y, k). Then 
we define 

yo = -—yo , y-i = -— 

and the symmetry properties follow from the evenness of J, a^^ and from Theo- 
rem lG.ll and Lemma rG.2l The zero order pseudodifferential operator Qq satisfies 
H3.22|l and from Proposition lF. II we have 

||Sow||r < ce||it||r, 0<r<p — 4 

WGouWs < c{e\\u\\, + Ei\\u\\o), 0<s</-4, (G.32) 
while the operator Q^i satisfies 

< C£||w||r-1, 1 < f < P — 8, 

< c(e||w||,_i+£;,|l^llo), l<s</-8. (G.33) 



We then deduce the estimates (|3.23|l in using I|(t.4|I . 

Now our task is to calculate the symbols Gj. It is convenient to introduce 
the scalar 1{Y) and linear form Il{Y, k) defined by 



I = J/Vdct G, U = G-^\7yV ■ k 

Recall the identity 

a33 0,jmkjkm — a^mkn\ 

l<j,m<2 l<m<2 

= det A((A-i)22fc? - 2(A-i)i2fci6 + (A-i)iife 



Noting that det A = J, «v = j juj^jui 

, 2 



aj^kjkm-[ ^3mkrn) = (B^ Bl )22 A:? "2 (1^ Bi ) 12 fci fc2 + (B^ »! ) n fc| . 



033 

l<j",m<2 l<m<2 



106 



Hence, by the definition of the metric tensor G, 



l<j",m<2 l<m<2 



which yields 

033 X! %m^j fcm-( ^ a^mkm^ = (dct G)G"^/c • /c. (G.34) 

l<j",m<2 l<m<2 

Substituting this relation into D (sec (|G.10|) ') finally gives 

D(Y,k) ^ VdetGGi(r,fc). (G.35) 

Noting that 5i = we obtain the needed formula H3.18|l . The calculation of 
Go is more delicate task. Since 033 — det G/J, formula IjG.Sp yields 

G . ^ ^ i ^ ^ 
{2IGi)Gq ^iWkG^WyG^ -y^dyMs-i^ dy^a^mkm- 

033 , 033 . , 

j = l 3,m=l 

It follows from the definition of the form 11 that 

2 2 
^ a3„ifc„ = -033!!, y^^dy.azj = - divy(a33Vfcn). (G.36) 

m—l J — 1 

From this, (|G.9|) . and (|G.35p we conclude that 

G^ = IGi - iU, Gt = IGi - ffl, 

and hence 

(2IGi)G^ =iVfc(IGi-m)Vy(-IGi-zn) — -{IGi+iIl) divy(a33Vfen)+ 

033 

^ 2 

^ ^ dyjO'jmkm- (G.37) 
y ,m— 1 

Next differentiating both sides of l|G.34|l with respect to kj we arrive to 

2 2 
2033 ^ ajmkm - '2a3j ^ a3„fc„i = (det G)dk^Gl, 

m—l m—l 

which along with ljG.36|) and the identity det 6/033 = J leads to 
2 

ajmkm = 033X19^^11 + JGi9fc^.Gi. 

m—l 
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Substituting this expression into ljG.37|l we finally obtain 

(2IGi)Go- =zVfc(IGi-m)Vy(-IGi-zn)- — (IGi+m) divy(a33Vfen)+ 

033 

— divy f a33nVfcn + JGiVfcGi) . (G.38) 
033 ^ ' 

Let us calculate the real part of Gq. It is easy to see that 

(2IG1) Re Go = IVfeGiVyn - VfenVy(IGi) - I^divy(a33Vfcn), 

033 



which along with the equality 1033 — \/dct G implies 



(2Gi) Re = VfcGi-Vyn-Vfcn-VyGi-Gi divy(Vfcn)-GiVfcn-Vy InVdet 
Noting that 

divy(Vfen) + VfcH • Vy In Vdet G ^ div(Vfcn), 

and recalling 



Go — — T"*^" ^ j2~^o ' (G.39) 



we obtain 



Re Go - |^(Vfen . VyGi - V^Gi • VyH) + ^div(V,n). 

From this and the identities 

VfeH^G-iVyT/, VfcGi = -^G-ift, VyGi = -^Vy(G-ifc-fc) 

we obtain the desired formula H3.19|l for the real part of Gq. Next (|G.38p yields 
(2IG1) Im G^ = -IVfcGi • Vy(IGi) - VfcH • VyH- 

— ndivy(a33Vfcn) + ^divy [nVfc(a33n) + jGiV^Gi] 

033 ^33 

— Gidivy(VfeGi) + Gi f — Vy J - IVyl) • VfcGi. 

a33 ^ 033 ^ 

Since J/ 033 — and 

1 . J2 1 



-Vy J - IVyl = — — VyVdet G = , VyVdet G, 

a33 (det G)3/2 Vdet G 

we have 

2 Im G^7 = I div (VfeGi). 
Recalling ljG.39|) we obtain H3.20|l and the theorem follows. 
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Invariant form of Re Go- In the rest of the section we prove the formula 
H3.24|l . Wc start with the calculation of the quadratic form 

It follows from (|3.21|l and the identity 'Gdy.Gs^^'G = —dy-G- that 

QiY, = Y.(^,^y^^)^ ■ e - ^dy^ (GC • Oq, - ■ L (g.40) 

where the vector field q = G^^Vyry. Thus we get 
^ 2 

which along with the equality ga/3 — dy^r ■ dy^r yields 

2 

Qa0 = iT.li^y^'^ - • dly^r. (G.41) 

i=i 

On the other hand, since r) = r • 63, the expression for q reads 

J- f i922dy,r - gi2dy.,r) ■ eg 
det G \ {-gudy^r + ggndy^r) ■ 63 

Now set 

a = dy^r, h = dy^r, c = a x b. (G.42) 

Noting that 

gii = a • a, gi2=a - b, 522 = b • b, 
b X c = (b • b)a — (a • b)b, a x c = (a • b)a — (a • a)b, 

we obtain 

(det G)qi (b X c) • 83, (det G)q2 = -(a x c) • 83. 
From this and the identity 

[(b X c) • 83]a - [(a x c) • 83]b - |cp83 = -(c • 83)0, 

which holds true for all a, b e and c = a x b, we conclude that 

(det G){qidy^r + q2dy^v) - |cpe3 = -(c • 83)0. 

Noting that jcp = det G and c • 83 = J we arrive at 

2 J 
qjdy^r - 83 = , n, 
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where n — c/|c| is the unit normal vector to S. Substituting this identity into 
HG.41|) gives Qap — — (n • S^^y^r) J/Vdet G which leads to 

QiY, e) = --=i=(Lc? + + ml). 

Vdet G 

Since 

G?(y, GO - GC • ^= + 2^^ei6 + 

we finally obtain 



'^^''^ ' -QjY, gqU - ^ + + ^1 . (G.43) 



2J2 iG2(y,G?)"^ ' ^'l J 2(i?e? + 2i^ei6 + Ge|)' 

Our next task is to express divq via the geometric characteristics of S. First 
we do this in the standard coordinates Y = X with fj = i]. In this case 



G-^ = ' /^^^ , dct G = l+|V77|^ J = 1, 

and q = (1 + |V77p)^^V77, which leads to the formula 

det G ,. Vdet Gl / V77 \ Vdet G h + h 

-divq= -div -j= , 1 = — , (G.44) 



2J2 ^ J 2 \^TTWW^ J 

where 6^ are the principal curvatures of E at the point r. Next note that V77 is 
a covariant vector field on E, hence G~^V?7 is a vector field on S. Since div is 
an invariant operator on the space of vector fields on S, the left side of ljG.44|l 
does not depend on the choice of coordinates, hence 



det G VdctGJi+ta VdctG LG - 2MF + A^'E^ ,^ 

^div q = = — — — . (G.45) 

2J2 ^ J 2 J 2{EG-F^) ^ ' 

Combining (|(t.43|I and (|G.45(I gives the desired identity H3.24II . If we define 
by n(0 the normal curvature of E in the direction ^ at a point r, then H3.24|l 
becomes 

^=1= Re Go(y,Ge) = + £2 - n(e)), 

which leads to 

Corollary G.3 Assume that the manifold E has a parametric representation 
r{Y) ^ (TF + V{X),fj{Y)), F e so that V and ij, which are not defined yet, 
satisfy all hypotheses of Theorem \3.5l Assume also that the parametric form 

&u:=^d=g,,uo{T + V), u{X) = u{Y{Xj), 
vdet G 

of the normal derivative operator is given for any hi- periodic smooth function 
u{Y). Then the manifold E is defined by the operator & up to a translation and 
a rotation of the embedding space. 
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Proof. Note that for all fc G Z^, 



1 _oVti..V iVtt-.V J 



"-'OO VdetG 
lim Re |e-™fe-^0e"'=^-^ - Gi{Y,k) \ ^ , ^ Re GQiY.k). 

Since Go is a homogeneous function of fc, it follows from this that the right hand 
sides of these equalities are defined by the operator © for all fc e and, in 
particular, for k = with an arbitrary ^ S R^. Hence the first fundamental 
form • f and the difference 61 + ^2 — "(0 ^"^^ completely defined by for 
all directions ^ at each point of E. Hence the principle curvatures of S are also 
defined by the operator ©. It remains to note that, by the Bonnet Theorem, the 
first fundamental form and the principal curvatures define E up to a translation 
and a rotation of the embedding space. ■ 



H Proof of Lemma 15.81 

To be able to compute all terms in H5.95|l . let us rewrite the system (|1.6|) . 
H1.7|l formally as a scalar equation for -0 and express the orders e and of 
the differential computed at > 3. Then the operator £ + ^ is up to 

order closely linked with the new form of this operator after applying the 
diffeomorphism computed at Lemma 13.81 

Indeed we can formally solve H1.7|l with respect to t] in powers of "0 as 

and replace rj by this expression in We then obtain a new scalar formal 

equation for ip, under the form fj.) = 0, where 

where 

= — 9^^ + ^^°^ 

and £2 and £3 represent quadratic and cubic terms in tp, respectively. Let us 
write the formal solution found at Theorem 12.31 for ei = 62 = e/2, under the 
form {N > 3) 

- £01 +eV2 + 0(e'), 
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where 

ipi — sinxi cosrx2, 

then we have the identities 

£0^1 = 0, 
£0^-2 +f 2(^-1, V'l) = 0, (H.2) 
-CoV'a + Mi'CiV'i + 2f2(V'i,'02) + f3('0i, "01: V"!) = 0. 

Now, we may observe that the operator H3.7|l we want to invert acts on Scf) = 
6ip — bdrj. Since b is 0{e) and Srj may be expressed formahy hnearly in terms of 
Sip in differentiating formaUy IjH.ll) . we have formally 

d^£{i;, + ni^, A^)) = fi) - n{£, fi) (H.3) 

where £(V', A*) is the linear operator we want to invert, k is the function we 
introduced at Theorem 13.41 and the operator TCiip, ^) is such that formally 

5^ = {i + n{^,v))H. 

and 7^(0, = 0. Since we set 

^j; = Vi^)+0(e^), iV>3 

Ai = +e^/ii(To) + O(e^), fio=fic{To), 



we have 
hence 



£:(^f),Mo+eVi) = 0(e3), 7^ = 0(e3), 



Making now the change of coordinates computed at Lemma 13.81 the new 

(2) 

expressions of operators £o, -Ci, £2{ip6 r), f3(V'ii'0ii ') take the following form 

new£o = £o + e4^^ + + 0{e^), 

new£i = 4"^ +££^^^ +0(£2), 

new£2(V'i,-) = £^°\4'\-)+e£i'\4"\-)+0{e'), 

newfalV-i,^!,-) = £i°\4°\4°\-) + 0(6), 

and the functions ipi , 1/^2 are transformed into 

new Vi = 4°^ + e4i'> + e^i^^ + 0{e^) 
newV'2 = +e42^ +0{e^). 

Moreover, we have thanks to Lemma [ 



k{Y) = l + eKi{Y)+e'^K2{Y) + 0{e^), 
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Hence, identities (|H.2|) lead to 



£o^(°) = 0, 
+ £o^i'^ = 0, (H.4) 

iioV'f +ff (V^''",^^'*") = 0, (H.5) 



£o4°^ + A^i^"'^^"^ + 2£f (./^('"^f ) +ff = 0. (H.6) 

We deduce from these formulae and from ljH.31) , that the linear operator obtained 
after the change of coordinates satisfies 

with 

K,£o +Sj^'^^ + 24°) (V'^") , •) + -CoHi (H.7) 

n,2o + ^1^^'^ + Sj^'^ - + + 2£^'\4,('\-) + 2ff •) + (^W, 

+3£f •) + {£« + 2£f (V^"', •)} + £0^2. 

We now compute the terms under the integral in 1)5. 95f) . First we observe (thanks 

to (|h31), (Hm ) 

hence 

+4£:f (V;(")>f) + £o7^i(4')+24°)) + 

Moreover 

i3(2)^(o) ^ £(2)^(0) + ^^£(0)^(0) + 2£(i)(V;("), V^(°)) + 2£f + 
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and in using again (|H.4|I . (|H.5|I . (|H.6|I . and (|H.7|) we obtain 



i3(2)^(0)_^(l)£^ 1^(1)^(0) 



Hence H5.95|l leads to 



@ 



/ (4' 



Since 



we finally obtain the result of Lemma [5.81 

I Fluid particles dynamics 

The kinematic and dynamic boundary conditions (|1.2I II. 3|) give two equations 
for two unknowns rj and '0- Assume for the moment that we know 77 i.e. the free 
surface E. The question is can we restore ip without solving PDE equations! 
The answer is yes: it suffices to solve a problem of moving a heavy single mass 
point along the free surface, or equivalently to find the corresponding geodesic 
flow on the surface with an appropriate metric. 

Let us begin with the consideration of the motion of a single mass point 
along the surface S = {2:3 = r]{X)}. Assuming that gravity fi acts in the 
—63 = (0, 0, —1) direction we can write the governing equations in the form 



where A is the Lagrange multiplier and n is a normal vector to S. Choosing com- 
ponents of X as generalized coordinates we rewrite equivalently these equation 
in the Lagrange form with the Lagrangian 



where GdX ■ dX is the first fundamental form of the free surface E with 



x + fies^Xn, X3^r]{X), 




More precisely, we have 



^d^L{X,X) - dxUX,X) =0. 



If we define the moments and Hamiltonian by 



y = GX, UiX,y)^-G-'y■y + ^Ir^{X) 
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then the governing equations can be rewritten in the Hamihon form 

X = dyU{X,y), y = -dxn{X,y). (LI) 

Next, suppose that a -generating function S{X) satisfies the Hamilton- Jacobi 
equation 

B.{X,VxS{X)) = h= const , (1.2) 
and the periodicity conditions 

S{X + 27rei) - 27r = S{X + —63) = S{X). (1.3) 

T T 

Suppose also that X{t) is a solution of the equations 

X = G-\X)WxS{X), (1.4) 

then, it is known that {X{t), y{t)), with y{t) = 'S/xS{X{t)), is a solution of the 
Hamiltonian system and the surface y — WxS{X) is an invariant manifold 
of HI.1|) . the flow being defined by (|I.4|I . 

Finally note that due the periodicity conditions, the mapping X {X, WxS{X)) 
defines an embedding of the torus M^/F into M^/F x M^. Therefore, {y = 
VxS{X)} is an invariant torus of lying on the energy surface H = h. In 
coordinates X the Hamiltonian flow on the torus is given by (|I.4|I 

Let us turn to the diamond wave problem. Set 

(p* (x) ^ U.0 ■ X + ip(x), and ij;* (X) ^ Lp* (xi,X2,ri{X)), 

and recall |uo| — 1. In these notations kinematic condition H1.2|) and dynamic 
condition H1.3|l can be rewritten in the equivalent form 

VxV*{xi.^2.xz)^<G{X)-^Vxr{X) for x^^t^{X), (1.5) 

^G-^Vxr -^xr +m=\- (1-6) 

By construction the equation x = V<y9* {x) determines the trajectories of liquid 
particles. From (|I.5|I . such a particle moving along the free surface satisfies 

X = G-\X)Vxr{X). (L7) 

On the other hand, equation Ijl.Bp reads 

H(X,VxV'*) - 1/2. 

Hence S{X) — '4>*{X) is a generating function for the Hamiltonian system Ijl.ip . 

From this we conclude that trajectories of liquid particles X{t), which are 
defined hy X = G~^VxV'*(-^)) along with y{t) = \/xS{X{t)) serve as solutions 
of (|I.1|I and belongs to the invariant torus {y = WxS{X)}. In other words, 
they coincide with projections {X, y) X oi solutions to ljl.l|) belonging to the 
invariant torus {y = Vi/j*{X)} C {H = 1/2}. Moreover, since by (|3.1|l we have 
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V = G ^'S/x4'*7 they also coincide with the integral curves of the vector field 
V. 

Finally note that, by the Maupertuis principle, the projections {X, y) X 
of solutions of belonging to the energy surface H — 1/2, coincide with the 
geodesies on the manifold S endowed with the Jacobi metric 

ds^ = {1/2- ^ri{X))<G{X)dX ■ dX = 2{l/2-JJ)T{X,dX). (1.8) 

Hence the integral curves of the vector field V form the geodesic flow associated 
with the metric H1.8() . 

Corollary I.l Suppose that rj is an arbitrary bi-periodic smooth function so 
that the hamiltonian system (|L1|) has an invariant torus {y = W xS{Xy\ with 
a smooth generating function S{X) satisfying the periodicity conditions p.3p . 
Then the solution if* of the Cauchy problem 

ip*ix)^S{X), dnip*{x) = for X3 = riiX), 

for the Laplace equation, satisfies kinematic and dynamic conditions (|1.2|l . (ll.3|l . 
(The local existence of such a solution follows from the Cauchy-Kowalewski the- 
orem, and the existence and boundedness in the lower half plane is true only for 
the "good" choice of generating function.) 
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